Saya menggunakan lme4 dalam R agar sesuai dengan model campuran
lmer(value~status+(1|experiment)))
di mana nilai kontinu, status dan percobaan adalah faktor, dan saya dapatkan
Linear mixed model fit by REML
Formula: value ~ status + (1 | experiment)
AIC BIC logLik deviance REMLdev
29.1 46.98 -9.548 5.911 19.1
Random effects:
Groups Name Variance Std.Dev.
experiment (Intercept) 0.065526 0.25598
Residual 0.053029 0.23028
Number of obs: 264, groups: experiment, 10
Fixed effects:
Estimate Std. Error t value
(Intercept) 2.78004 0.08448 32.91
statusD 0.20493 0.03389 6.05
statusR 0.88690 0.03583 24.76
Correlation of Fixed Effects:
(Intr) statsD
statusD -0.204
statusR -0.193 0.476
Bagaimana saya bisa tahu bahwa pengaruh status itu signifikan? R hanya melaporkan nilai- dan bukan nilai- .p
Jawaban:
Ada banyak informasi tentang topik ini di FAQ GLMM . Namun, dalam kasus khusus Anda, saya sarankan menggunakan
karena Anda tidak memerlukan barang apa pun yang
lmer
menawarkan (kecepatan lebih tinggi, penanganan efek acak silang, GLMMs ...).lme
harus memberikan persis sama koefisien dan varians estimasi tetapi juga akan df menghitung dan p-nilai untuk Anda (yang melakukan masuk akal dalam "klasik" desain seperti Anda tampaknya memiliki). Anda mungkin juga ingin mempertimbangkan istilah acak~status|experiment
(memungkinkan variasi efek status di seluruh blok, atau setara dengan interaksi status-per-eksperimen). Poster di atas juga benar bahwat
statistik Anda sangat besar sehingga nilai-p Anda pasti akan <0,05, tetapi saya bisa membayangkan Anda ingin nilai-p "nyata".sumber
lmer
bisa dengan mudah melaporkan jenis-nilai p yang sama tetapi tidak untuk alasan yang valid. Saya kira itu adalah komentar bahwa ada nilai-p "nyata" di sini yang mengganggu saya. Anda dapat berargumen bahwa Anda dapat menemukan satu kemungkinan cutoff, dan bahwa cutoff wajar apa pun dilewati. Tapi Anda tidak bisa berdebat ada nilai p nyata.summary(m1)
bukan (saya menggunakan ini dengan paket nlme)Anda bisa menggunakan paket lmerTest . Anda cukup menginstal / memuatnya dan model lmer akan diperpanjang. Jadi misalnya
akan memberi Anda hasil dengan nilai-p. Jika nilai-p adalah indikasi yang tepat agak sedikit diperdebatkan, tetapi jika Anda ingin memilikinya, ini adalah cara untuk mendapatkannya.
sumber
Jika Anda dapat menangani meninggalkan nilai-p ( dan Anda harus ), Anda dapat menghitung rasio kemungkinan yang akan mewakili bobot bukti untuk efek status melalui:
sumber
Masalahnya adalah bahwa perhitungan nilai-p untuk model-model ini tidak sepele, lihat diskusi di sini sehingga penulis
lme4
paket ini sengaja memilih untuk tidak memasukkan nilai-p dalam output. Anda mungkin menemukan metode penghitungan ini, tetapi mereka tidak selalu benar.sumber
Pertimbangkan apa yang Anda tanyakan. Jika Anda hanya ingin tahu apakah nilai p keseluruhan untuk efek status melewati semacam nilai cutoff sewenang-wenang, seperti 0,05, maka itu mudah. Pertama, Anda ingin mengetahui efek keseluruhannya. Anda bisa mendapatkannya dari
anova
.Sekarang Anda memiliki nilai F. Anda bisa mengambilnya dan mencarinya di beberapa tabel F. Pilih saja denom serendah mungkin. derajat kebebasan. Batasnya akan ada sekitar 20. F Anda mungkin lebih besar dari itu tetapi saya bisa salah. Bahkan jika tidak, lihat jumlah derajat kebebasan dari perhitungan ANOVA konvensional di sini menggunakan jumlah percobaan yang Anda miliki. Menempel nilai itu di Anda turun menjadi sekitar 5 untuk cutoff. Sekarang Anda dengan mudah lulus dalam studi Anda. Df 'benar' untuk model Anda akan menjadi sesuatu yang lebih tinggi dari itu karena Anda memodelkan setiap titik data yang bertentangan dengan nilai agregat yang akan dimodelkan oleh ANOVA.
Jika Anda benar-benar menginginkan nilai-p yang tepat, tidak ada yang namanya kecuali Anda bersedia membuat pernyataan teoretis tentang hal itu. Jika Anda membaca Pinheiro & Bates (2001, dan mungkin beberapa buku lagi tentang masalah ini ... lihat tautan lain dalam jawaban ini) dan Anda keluar dengan argumen untuk df tertentu maka Anda dapat menggunakannya. Tapi Anda sebenarnya tidak mencari nilai p yang tepat. Saya menyebutkan ini karena Anda seharusnya tidak melaporkan nilai-p yang tepat, hanya saja cutoff Anda dilewati.
Anda harus benar-benar mempertimbangkan jawaban Mike Lawrence karena seluruh gagasan hanya bertahan dengan titik lulus untuk nilai-p sebagai informasi akhir dan paling penting untuk diekstrak dari data Anda pada umumnya salah arah (tetapi mungkin tidak ada dalam kasus Anda karena kami tidak Saya tidak punya cukup informasi untuk diketahui). Mike menggunakan versi hewan peliharaan dari perhitungan LR yang menarik, tetapi mungkin sulit untuk menemukan banyak dokumentasi di dalamnya. Jika Anda melihat pemilihan model dan interpretasi menggunakan AIC Anda mungkin menyukainya.
sumber
Sunting: Metode ini tidak lagi didukung di versi lme4 yang lebih baru. Gunakan paket lmerTest seperti yang disarankan dalam jawaban ini oleh pbx101 .
Ada posting di daftar R oleh penulis lme4 tentang mengapa nilai-p tidak ditampilkan. Dia menyarankan menggunakan sampel MCMC sebagai gantinya, yang Anda lakukan menggunakan pvals.fnc dari paket languageR:
Lihat http://www2.hawaii.edu/~kdrager/MixedEffectsModels.pdf untuk contoh dan detail.
sumber
Apakah Anda tertarik mengetahui apakah efek gabungan dari
status
berpengaruh signifikanvalue
? Jika demikian, Anda dapat menggunakanAnova
fungsi dalamcar
paket (jangan bingung dengananova
fungsi di pangkalanR
).Lihatlah
?Anova
setelah memuatcar
paket.sumber
car::Anova()
menghindari masalah lengket seputar perhitungan nilai-p yang ditautkan Michelle?anova
Perintah reguler akan memberi Anda nilai F.Fungsi
pvals.fnc
ini tidak lagi didukung oleh lme4. Menggunakan paket paket lmerTest, dimungkinkan untuk menggunakan metode lain untuk menghitung nilai-p, seperti perkiraan Kenward-Rogersumber
Hanya dengan memuat paket afex akan mencetak nilai-p dalam output fungsi lmer dari paket lme4 (Anda tidak perlu menggunakan afex; muat saja):
Ini akan secara otomatis menambahkan kolom nilai p ke output lmer (model Anda) untuk efek tetap.
sumber