Pertanyaan yang diberi tag «r»

Harap berikan kode R yang memungkinkan seseorang untuk melakukan ANOVA antar-subyek dengan -3, -1, 1, 3 kontras. Saya mengerti ada perdebatan mengenai tipe Sum of Squares (SS) yang sesuai untuk analisis semacam itu. Namun, sebagai tipe standar SS yang digunakan dalam SAS dan SPSS (Tipe III) dianggap standar di daerah …

26 r  anova  contrasts  sums-of-squares 

Saya baru mengenal R, memerintahkan regresi logistik, dan polr. Bagian "Contoh" di bagian bawah halaman bantuan untuk polr (yang cocok dengan model regresi logistik atau probit ke respons faktor yang dipesan) menunjukkan options(contrasts = c("contr.treatment", "contr.poly")) house.plr <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data = …

26 r  logistic 

Dalam sebuah pertanyaan di tempat lain di situs ini, beberapa jawaban menyebutkan bahwa AIC setara dengan validasi silang leave-one-out (LOO) dan bahwa BIC setara dengan validasi silang K-fold. Apakah ada cara untuk secara empiris menunjukkan ini dalam R sehingga teknik yang terlibat dalam LOO dan K-fold diperjelas dan didemonstrasikan setara …

26 r  aic  cross-validation  bic 

Saya memiliki matriks korelasi pengembalian keamanan yang determinannya nol. (Ini agak mengejutkan karena matriks korelasi sampel dan matriks kovarian yang sesuai secara teoritis pasti positif.) Hipotesis saya adalah bahwa setidaknya satu sekuritas secara linear tergantung pada sekuritas lainnya. Apakah ada fungsi dalam R yang secara berurutan menguji setiap kolom sebuah …

26 r  correlation  pca  linear-model  svd 

Saya menjalankan model regresi LOESS dalam R, dan saya ingin membandingkan output dari 12 model yang berbeda dengan ukuran sampel yang bervariasi. Saya dapat menjelaskan model yang sebenarnya lebih terinci jika itu membantu dengan menjawab pertanyaan. Berikut adalah ukuran sampel: Fastballs vs RHH 2008-09: 2002 Fastballs vs LHH 2008-09: 2209 …

26 r  regression  loess 

Apakah mungkin untuk menemukan nilai-p dalam korelasi pearson dalam R? Untuk menemukan korelasi pearson, saya biasanya melakukan ini col1 = c(1,2,3,4) col2 = c(1,4,3,5) cor(col1,col2) # [1] 0.8315218 Tetapi bagaimana saya bisa menemukan nilai p ini?

26 r  correlation  p-value  pearson-r 

Saya telah menjelajahi banyak situs bantuan dan masih bingung tentang cara menentukan istilah bersarang yang lebih rumit dalam model campuran juga. Saya juga bingung karena penggunaan :dan /dan |dalam menentukan interaksi dan bersarang dengan faktor acak yang digunakan lmer()dalam lme4paket di R. Untuk keperluan pertanyaan ini, anggap saya telah secara …

26 r  mixed-model  lme4-nlme 

Ini mirip dengan pertanyaan metode pengambilan sampel Caret , meskipun itu benar-benar tidak pernah menjawab bagian dari pertanyaan ini dengan cara yang disepakati. fungsi kereta caret menawarkan cvdan repeatedcv. Apa perbedaan dalam mengatakan melakukan: MyTrainControl=trainControl( method = "cv", number=5, repeats=5 ) vs. MyTrainControl=trainControl( method = "repeatedcv", number=5, repeats=5 ) Saya …

26 r  machine-learning  caret 

Terinspirasi oleh ceramah Peter Donnelly di TED , di mana ia membahas berapa lama untuk pola tertentu muncul dalam serangkaian lemparan koin, saya membuat skrip berikut dalam R. Diberikan dua pola 'hth' dan 'htt', itu menghitung berapa lama (yaitu berapa banyak lemparan koin) rata-rata sebelum Anda menekan salah satu pola …

26 r  probability  stochastic-processes 

Saya sebelumnya telah menggunakan pro perkiraan untuk meramalkan seri waktu univariat, tetapi saya mengalihkan alur kerja saya ke R. Paket ramalan untuk R berisi banyak fungsi yang bermanfaat, tetapi satu hal yang tidak dilakukannya adalah segala jenis transformasi data sebelum menjalankan otomatis .arima (). Dalam beberapa kasus ramalan pro memutuskan …

26 r  time-series  data-transformation  forecasting  arima 

Pertimbangkan model rintangan yang memprediksi data jumlah ydari prediktor normal x: set.seed(1839) # simulate poisson with many zeros x <- rnorm(100) e <- rnorm(100) y <- rpois(100, exp(-1.5 + x + e)) # how many zeroes? table(y == 0) FALSE TRUE 31 69 Dalam hal ini, saya memiliki data jumlah …

25 r  count-data  zero-inflation 

Kutipan blok di bawah ini, dari para pemimpin di bidang pemodelan efek campuran, mengklaim bahwa mengoordinasikan pergeseran dalam model dengan korelasi nol antara efek acak (model 'ZCP') mengubah prediksi model. Tapi, bisakah seseorang menguraikan atau lebih jauh membenarkan klaim mereka? Pernyataan yang dipertanyakan berasal dari makalah Bates et al 2015 …

25 r  mixed-model  lme4-nlme 

Saya ingin menunjukkan bagaimana nilai-nilai variabel tertentu (~ 15) berubah seiring waktu, tetapi saya juga ingin menunjukkan bagaimana variabel berbeda satu sama lain di setiap tahun. Jadi saya membuat plot ini: Tetapi bahkan ketika mengubah skema warna atau menambahkan jenis garis / bentuk yang berbeda ini terlihat berantakan. Apakah ada …

25 r  data-visualization 

The lmerTestpaket menyediakan anova()fungsi untuk model linear campuran dengan opsional pendekatan Satterthwaite (default) atau Kenward-Roger dari derajat kebebasan (df). Apa perbedaan antara kedua pendekatan ini? Kapan memilih yang mana?

25 r  anova  mixed-model  lme4-nlme  degrees-of-freedom 

Saat ini saya menjalankan beberapa model linear efek campuran. Saya menggunakan paket "lme4" di R. Model saya mengambil bentuk: model <- lmer(response ~ predictor1 + predictor2 + (1 | random effect)) Sebelum menjalankan model saya, saya memeriksa kemungkinan multikolinieritas antara prediktor. Saya melakukan ini dengan: Buat kerangka data prediksi dummy_df …

25 r  correlation  mixed-model  lme4-nlme  multicollinearity