Pertanyaan yang diberi tag «r»

Saya baru dalam bidang statistik dan saya mencoba memahami perbedaan antara ANOVA dan regresi linier. Saya menggunakan R untuk mengeksplorasi ini. Saya membaca berbagai artikel tentang mengapa ANOVA dan regresi berbeda tetapi masih sama dan bagaimana dapat divisualisasikan dll. Saya pikir saya cantik di sana tetapi satu bit masih hilang. …

9 r  regression  anova 

Saya ingin menilai keakuratan tes normalitas pada ukuran sampel yang berbeda dalam R (Saya menyadari bahwa tes normalitas dapat menyesatkan ). Sebagai contoh, untuk melihat tes Shapiro-Wilk, saya sedang melakukan simulasi berikut (dan juga merencanakan hasilnya) dan berharap bahwa ketika ukuran sampel meningkatkan kemungkinan menolak nol menurun: n <- 1000 …

9 r  simulation  power-analysis  normality-assumption 

Saya mencoba memperkirakan efek 2 obat ( drug1, drug2) pada kemungkinan pasien jatuh ( event). Para pasien dapat terjatuh lebih dari satu kali dan dapat memakai atau melepas obat kapan saja. Pertanyaan saya adalah bagaimana data harus disusun berkaitan dengan periode waktu (hari), khususnya apakah perlu ada tumpang tindih antara …

9 r  survival  cox-model 

Saya sedang mengerjakan estimasi vektor auto-regression (VAR) dan impulse response function (IRFs) berdasarkan data panel dengan 33 orang lebih dari 77 kuartal. Bagaimana seharusnya jenis situasi ini dianalisis? Algoritma apa yang ada untuk tujuan ini? Saya lebih suka melakukan analisis ini dalam R, jadi jika ada yang akrab dengan kode …

9 r  econometrics  panel-data  var  impulse-response 

Saya memiliki model regresi logistik (cocok melalui glmnet di R dengan regularisasi jaring elastis), dan saya ingin memaksimalkan perbedaan antara positif benar dan positif palsu. Untuk melakukan ini, prosedur berikut muncul dalam pikiran: Sesuai dengan model regresi logistik standar Dengan menggunakan ambang prediksi 0,5, identifikasi semua prediksi positif Tetapkan bobot …

9 r  regression  logistic  classification  glmnet 

Saya tertarik menemukan prosedur untuk mensimulasikan data yang konsisten dengan model mediasi yang ditentukan. Menurut kerangka model persamaan struktural linear umum untuk menguji model mediasi pertama kali diuraikan oleh Barron dan Kenny (1986) dan dijelaskan di tempat lain seperti Judd, Yzerbyt, & Muller (2013) , model mediasi untuk hasil YYY …

9 r  regression  simulation  random-generation  mediation 

Dalam Dixon, Coles ( 1997 ), mereka telah menggunakan estimasi kemungkinan maksimum untuk dua model Poisson independen yang dimodifikasi dalam (4.3) untuk memodelkan skor dalam sepak bola. Saya mencoba menggunakan R untuk "mereproduksi" alfa dan beta serta parameter efek rumah (hal 274, Tabel 4) tanpa menggunakan paket apa pun (menggunakan …

9 r  modeling  maximum-likelihood  games 

Saya menganalisis set data tentang komunitas intertidal. Data berupa persen tutupan (dari rumput laut, teritip, kerang, dll) dalam kuadrat. Saya terbiasa memikirkan analisis korespondensi (CA) dalam hal jumlah spesies , dan analisis komponen utama (PCA) sebagai sesuatu yang lebih berguna untuk tren lingkungan (bukan spesies) linier. Saya belum benar-benar beruntung …

9 r  pca  multivariate-analysis  correspondence-analysis 

Saya memiliki data berikut yang terletak di sini . Saya mencoba menghitung interval kepercayaan 95% pada kemurnian rata-rata ketika persentase hidrokarbon adalah 1,0. Di R, saya masukkan yang berikut ini. > predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95) fit lwr upr 1 89.66431 87.51017 91.81845 Namun, bagaimana saya bisa mendapatkan hasil ini sendiri? …

9 r  regression  confidence-interval  prediction-interval 

Saya punya contoh data yang diatur sebagai berikut: Volume <- seq(1,20,0.1) var1 <- 100 x2 <- 1000000 x3 <- 30 x4 = sqrt(x2/pi) H = x3 - Volume r = (x4*H)/(H + Volume) Power = (var1*x2)/(100*(pi*Volume/3)*(x4*x4 + x4*r + r*r)) Power <- jitter(Power, factor = 1, amount = 0.1) plot(Volume,Power) …

9 r  hypothesis-testing  nonlinear-regression  goodness-of-fit  heteroscedasticity 

Ini mungkin pertanyaan yang konyol, tetapi ketika membuat model dengan tanda sisipan dan menggunakan sesuatu seperti LOOCVatau (bahkan lebih tepatnya) LGOCV, apa manfaat membagi data ke dalam kereta dan set tes jika ini pada dasarnya merupakan langkah validasi silang ngomong-ngomong Saya membaca beberapa pertanyaan terkait dan mereka menyarankan bahwa beberapa …

9 r  machine-learning  cross-validation  caret 

Saya punya model seperti ini: require(nlme) set.seed(123) n <- 100 k <- 5 cat <- as.factor(rep(1:k, n)) cat_i <- 1:k # intercept per kategorie x <- rep(1:n, each = k) sigma <- 0.2 alpha <- 0.001 y <- cat_i[cat] + alpha * x + rnorm(n*k, 0, sigma) plot(x, y) m1 …

9 r  mixed-model  aic 

Saya melakukan meta-analisis ukuran efek d dalam R menggunakan paket metafor. d mewakili perbedaan dalam skor memori antara pasien dan sehat. Namun beberapa penelitian melaporkan hanya subskala ukuran minat d (misalnya beberapa skor memori yang berbeda atau skor dari tiga blok pengujian memori yang terpisah). Silakan lihat kumpulan data dummy …

9 r  meta-analysis  effect-size  meta-regression 

Saya ingin memahami penggunaan simulasi Monte Carlo dalam chisq.test()fungsi di R. Saya memiliki variabel kualitatif yang memiliki 128 level / kelas. Ukuran sampel saya adalah 26 (saya tidak dapat mengambil sampel lebih banyak "individu"). Jadi jelas, saya akan memiliki beberapa level dengan 0 "individu". Tetapi kenyataannya adalah bahwa saya hanya …

9 r  chi-squared  monte-carlo 

Saya memiliki randomForestmodel klasifikasi yang baik yang ingin saya gunakan dalam aplikasi yang memprediksi kelas kasus baru. Kasing baru memiliki nilai yang pasti hilang. Prediksi tidak akan berfungsi seperti itu untuk NAS. Bagaimana saya harus melakukan ini? data(iris) # create first the new case with missing values na.row<-45 na.col<-c(3,5) case.na<-iris[na.row,] …

9 r  random-forest  prediction  missing-data