Pertanyaan yang diberi tag «r»

Saya telah melihat ada berbagai pertanyaan mengenai interpretasi dan pembangunan gams, yang tampaknya menggambarkan kesulitan bagi non-ahli statistik untuk berurusan dengan itu. Sayangnya, dari tidak ada thread atau tutorial yang saya baca, saya bisa mendapatkan pemahaman yang jelas tentang bagaimana membangun model yang bermakna. Saat ini saya sedang mempelajari pengaruh …

9 r  model-selection  gam 

Saya menggunakan adonis()fungsi dalam veganpaket untuk menentukan 1) jika spesies inang yang muncul secara bersamaan bervariasi dalam komunitas mikroba mereka di beberapa situs, dan 2) jika situs berbeda. Saya telah memeriksa semua posting di CV dan SO, dan tidak ada jawaban yang jelas tentang bagaimana menentukan signifikansi beberapa faktor menggunakan …

9 r  permutation-test  sums-of-squares  blocking 

Saya telah ditanya pertanyaan tentang model linier yang dibuat dengan R's lm: "Apakah regresi menggunakan kuadrat terkecil linier atau non-linear?" Saya mencari sedikit dan [berpikir bahwa saya] memahami perbedaan antara keduanya, tetapi tidak dapat menemukan bukti penggunaan R dari linear kuadrat terkecil di lm(yang saya pikir itu menggunakan). Saya menyisir …

9 r  linear-model  lm 

The arimaxfungsi dalam TSApaket adalah untuk pengetahuan saya satu-satunya Rpaket yang akan cocok dengan fungsi transfer untuk model intervensi. Tidak memiliki fungsi prediksi meskipun yang kadang-kadang diperlukan. Apakah langkah-langkah berikut untuk mengatasi masalah ini, memanfaatkan forecastpaket yang sangat baik ? Apakah interval prediksi benar? Dalam contoh saya, kesalahan std adalah …

9 r  time-series  forecasting  arima  intervention-analysis 

Bagaimana cara menghitung mean, varians, median, standar deviasi dan modus dari distribusi? Jika saya secara acak menghasilkan angka yang membentuk distribusi normal, saya telah menetapkan mean sebagai m=24.2standar deviasi seperti sd=2.2: > dist = rnorm(n=1000, m=24.2, sd=2.2) Maka saya bisa melakukan hal berikut: Berarti: > mean(dist) [1] 24.17485 Perbedaan: > …

9 r  distributions  mean 

Misalkan dalam analisis regresi dalam R, saya memiliki jenis variabel variabel independen dengan 3 level dalam dataset kereta saya. Tetapi dalam data uji mengatur bahwa variabel faktor yang sama memiliki 5 level. Oleh karena itu saya tidak dapat memprediksi nilai respons untuk dataset uji. Apa yang harus dilakukan dalam kasus …

9 r  regression  categorical-data  cross-validation  rare-events 

Data mentah saya terdiri dari seri waktu 60 hari dengan tren menurun. Data mingguan sehingga frekuensi diatur ke 7. Saya menghitung selisih data yang terlihat seperti ini Ketika saya menjalankan plot ACF dan PACF pada perbedaannya, saya sepertinya mendapatkan hasil yang kontradiktif? ACF menunjukkan dampak positif dari term lagged pertama …

9 r  time-series  autocorrelation 

Saat memilih jumlah simpul yang tepat untuk GAM, orang mungkin ingin mempertimbangkan jumlah data dan peningkatan pada sumbu x. Bagaimana jika kita memiliki 100 peningkatan pada sumbu x dengan 1000 titik data pada setiap kenaikan. Info di sini mengatakan: Jika tidak disediakan maka simpul spline ditempatkan secara merata di seluruh …

9 r  gam  splines 

Saya memiliki tiga variabel ekonomi makro (ICS - sentimen konsumen, ER - tingkat pekerjaan, DGO - pesanan barang tahan lama) dan telah menjalankan tes kausalitas Granger dalam R pada mereka. Saya tidak benar-benar tahu bagaimana menafsirkan hasil tes Granger. Adakah yang bisa membantu saya memahami hasilnya? Saya tahu bahwa kami …

9 r  granger-causality 

Saya memiliki beberapa kesulitan untuk mendapatkan hasil yang setara antara aovmodel tindakan berulang antara-dalam dan lmermodel campuran. Data dan skrip saya terlihat sebagai berikut data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id FITNESS TEST PULSE 1 1 pilates CYCLING 91 2 2 pilates CYCLING 82 3 3 pilates CYCLING 65 4 4 pilates CYCLING …

9 r  anova  mixed-model  repeated-measures  lme4-nlme 

Saya ingin menghitung tingkat kejadian untuk disajikan di sepanjang rasio bahaya untuk menyajikan ukuran risiko relatif dan absolut. Saya melihat dalam penelitian lain bahwa tingkat kejadian seperti itu dapat dihitung menggunakan model poisson dengan waktu tindak lanjut dalam model sebagai offset. Jadi saya mencobanya di R sebagai berikut: library(survival) # …

9 r  survival  cox-model  poisson-regression  incidence-rate-ratio 

library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1 0.92498 0.18330 5.046 4.51e-07 *** period2 -0.06698 0.22845 -0.293 …

9 r  binomial  logit  lme4-nlme 

Dari apa yang saya mengerti, argumen cp ke rpartfungsi membantu pra-pangkas pohon dengan cara yang sama dengan argumen menitplit atau minbucket. Yang tidak saya mengerti adalah bagaimana nilai CP dihitung. Sebagai contoh df<-data.frame(x=c(1,2,3,3,3,4), y=as.factor(c(TRUE, TRUE, FALSE, TRUE, FALSE, FALSE)), method="class") mytree<-rpart(y ~ x, data = df, minbucket = 1, minsplit=1) …

9 r  cart  rpart 

Saya berjuang dengan poin proyeksi dalam analisis diskriminan linier (LDA). Banyak buku tentang metode statistik multivariat menggambarkan gagasan LDA dengan gambar di bawah ini. Deskripsi masalahnya adalah sebagai berikut. Pertama, kita perlu menggambar batas keputusan, menambahkan garis tegak lurus dan kemudian memplot proyeksi titik data di atasnya. Saya ingin tahu …

9 r  machine-learning  classification  discriminant-analysis 

Berikut datanya: > tires <- data.frame(Wear = c(17, 14, 12, 13, 14, 14, 12, 11, 13, 13, 10, 11, 13, 8, 9, 9), Brand = rep(LETTERS[1:4], 4), Car = as.character(as.roman(rep(1:4, each = 4)))) > tires Wear Brand Car 1 17 A I 2 14 B I 3 12 C I …

9 r  anova  dataset  degrees-of-freedom