Dalam kode berikut ini saya melakukan regresi logistik pada data yang dikelompokkan menggunakan glm dan "dengan tangan" menggunakan mle2. Mengapa fungsi logLik di R memberi saya kemungkinan logLik (fit.glm) = - 2.336 yang berbeda dari yang logLik (fit.ml) = - 5.514 saya dapatkan dengan...
Saya hanya mencoba untuk menghitung ulang dengan dnorm () log-kemungkinan yang disediakan oleh fungsi logLik dari model lm (dalam R). Ini berfungsi (hampir sempurna) untuk jumlah data yang tinggi (mis. N = 1000): > n <- 1000 > x <- 1:n > set.seed(1) > y <- 10 + 2*x +...
Ketika mempelajari pemilihan kovarians, saya pernah membaca contoh berikut. Sehubungan dengan model berikut: Matriks kovariansnya dan matriks kovarians terbalik diberikan sebagai berikut, Saya tidak mengerti mengapa independensi dan diputuskan oleh kovarians terbalik di sini?yxxxyyy Apa...
Pertanyaan / topik ini muncul dalam diskusi dengan seorang kolega dan saya mencari beberapa pendapat tentang ini: Saya memodelkan beberapa data menggunakan regresi logistik efek acak, lebih tepatnya regresi logistik intersep acak. Untuk efek tetap, saya memiliki 9 variabel yang menarik dan...
Sebagai contoh, dalam manual BUGS atau buku yang akan datang oleh Lee dan Wagenmakers ( pdf ) dan di banyak tempat lain jenis notasi digunakan yang bagi saya tampaknya sangat fleksibel karena dapat digunakan untuk menggambarkan secara ringkas sebagian besar model statistik. Contoh dari notasi ini...
Dalam model campuran, kami mengasumsikan efek acak (parameter) adalah variabel acak yang mengikuti distribusi normal. Itu terlihat sangat mirip dengan metode Bayesian, di mana semua parameter dianggap acak. Jadi, apakah model efek acak jenis kasus khusus metode
ESL bab 2.4 tampaknya mengklasifikasikan regresi linier sebagai "berbasis model", karena mengasumsikan , sedangkan tidak ada perkiraan serupa yang dinyatakan untuk k-tetangga terdekat. Tapi bukankah kedua metode membuat asumsi tentang ?f ( x )f( x ) ≈ x ⋅ βf(x)≈x⋅βf(x) \approx x\cdot\betaf( x...
Saya baru dalam analisis bertahan hidup dan saya baru-baru ini belajar bahwa ada berbagai cara untuk melakukannya dengan tujuan tertentu. Saya tertarik pada implementasi aktual dan kesesuaian metode ini. Saya diberi Cox Proportional-Hazards tradisional , model waktu kegagalan yang dipercepat dan...
Saya menganalisis hasil percobaan waktu reaksi dalam R. Saya menjalankan pengukuran ANOVA berulang (1 faktor dalam subjek dengan 2 level dan 1 antara faktor subjek dengan 2 level). Saya menjalankan model campuran linear yang sama dan saya ingin meringkas hasil lmer dalam bentuk tabel ANOVA...
Saya tidak mengerti mengapa mengubah jaringan Bayesian menjadi grafik faktor baik untuk inferensi Bayesian? Pertanyaan saya adalah: Apa manfaat menggunakan grafik faktor dalam penalaran Bayesian? Apa yang akan terjadi jika kita tidak menggunakannya? Setiap contoh nyata akan
Bagaimana struktur varians-kovarians default untuk efek acak dalam glmeratau lmerdalam lme4paket? Bagaimana cara menentukan struktur varians-kovarians lain untuk efek-acak dalam kode? Saya tidak dapat menemukan informasi mengenai ini di
Saya hanya mencoba untuk mengimplementasikan algoritma Baum-Welch yang diskalakan dan saya telah mengalami masalah di mana variabel terbelakang saya, setelah penskalaan, berada di atas nilai 1. Apakah ini normal? Bagaimanapun, probabilitas tidak boleh lebih dari 1. Saya menggunakan faktor skala...
Dalam regresi logistik, rasio odds 2 berarti bahwa kejadian tersebut 2 kali lebih mungkin diberikan peningkatan satu unit dalam prediktor. Dalam regresi Cox, rasio bahaya 2 berarti peristiwa tersebut akan terjadi dua kali lebih sering pada setiap titik waktu mengingat peningkatan satu unit dalam...
Saya menganalisa beberapa data perilaku menggunakan lme4di R, sebagian besar mengikuti tutorial Bodo Winter baik , tapi saya tidak mengerti jika saya menangani interaksi dengan benar. Lebih buruk lagi, tidak ada orang lain yang terlibat dalam penelitian ini menggunakan model campuran, jadi saya...
Saya memiliki sampel 1,449 titik data yang tidak berkorelasi (r-squared 0,006). Ketika menganalisis data, saya menemukan bahwa dengan memecah nilai-nilai variabel independen menjadi kelompok-kelompok positif dan negatif, tampaknya ada perbedaan yang signifikan dalam rata-rata variabel dependen...
Saya memiliki model berikut yang m_plotdilengkapi dengan lme4::lmerefek acak silang untuk peserta ( lfdn) dan item ( content): Random effects: Groups Name Variance Std.Dev. Corr lfdn (Intercept) 172.173 13.121 role1 62.351 7.896 0.03 inference1 24.640 4.964 0.08 -0.30 inference2 52.366...
Sebelum mengajukan pertanyaan ini, saya mencari di situs kami dan menemukan banyak pertanyaan serupa, (seperti di sini , di sini , dan di sini ). Tetapi saya merasa pertanyaan-pertanyaan terkait itu tidak ditanggapi atau didiskusikan dengan baik, sehingga ingin mengajukan pertanyaan ini lagi. Saya...
Baik bayesglm()(dalam paket arm R) dan berbagai fungsi dalam paket MCMCpack ditujukan untuk melakukan estimasi Bayesian dari model linier umum, tetapi saya tidak yakin mereka benar-benar menghitung hal yang sama. Fungsi MCMCpack menggunakan rantai Markov Monte Carlo untuk mendapatkan sampel...
Dalam regresi linear saya menemukan hasil yang menyenangkan jika kita cocok dengan model E[Y]=β1X1+β2X2+c,E[Y]=β1X1+β2X2+c,E[Y] = \beta_1 X_1 + \beta_2 X_2 + c, kemudian, jika kita membuat standar dan memusatkan data , dan
Saya sungguh berharap bahwa saya telah mengutarakan pertanyaan ini sedemikian rupa sehingga dapat dijawab secara pasti - jika tidak, tolong beri tahu saya dan saya akan mencoba lagi! Saya juga harus menebak bahwa saya akan menggunakan R untuk analisis ini. Saya memiliki beberapa ukuran plant...