Proporsi varian yang dijelaskan dalam model efek-campuran

18

Saya tidak tahu apakah ini sudah ditanyakan sebelumnya, tetapi saya tidak menemukan apa-apa tentang itu. Pertanyaan saya adalah apakah ada yang bisa memberikan referensi yang baik untuk belajar bagaimana mendapatkan proporsi varian yang dijelaskan oleh masing-masing faktor tetap dan acak dalam model efek campuran.

Manuel Ramón
sumber
4
Pertanyaan bagus, tapi saya tidak punya (referensi untuk) jawaban yang bagus. Ada lebih dari satu tingkat variasi dalam model campuran, jadi ada lebih dari satu komponen varians untuk dijelaskan, ditambah lagi masih dapat diperdebatkan apakah efek acak benar-benar dapat dikatakan 'menjelaskan' varians. Saya pikir seluruh konsep 'proporsi varian dijelaskan' kurang berguna dalam model campuran.
onestop
Berikut adalah beberapa diskusi lebih lanjut tentang topik ini: stat.ethz.ch/pipermail/r-sig-mixed-models/2010q1/003363.html
user5475
1
Pendekatan Gelmans "Bayesian ANOVA" mungkin juga berguna.
N Brouwer

Jawaban:

12

Saya dapat memberikan beberapa referensi:

Xu, R. (2003). Mengukur variasi yang dijelaskan dalam model efek campuran linier. Statistik dalam Kedokteran , 22 , 3527-3541. DOI: 10.1002 / sim.1572

R2

Hössjer, O. (2008). Pada koefisien determinasi untuk model regresi campuran. Jurnal Perencanaan Statistik dan Inferensi , 138 , 3022-3038. DOI: 10.1016 / j.jspi.2007.11.010

R2

Selamat membaca!

Wolfgang
sumber
5

MuMInR21

#load packages
library(lme4)
library(MuMIn)

#Fit Model
m <- lmer(mpg ~ gear + disp + (1|cyl), data = mtcars)

#Determine R2:
r.squaredGLMM(m) 

       R2m       R2c 
 0.5476160 0.7150239  

Output untuk fungsi r.squaredGLMMmenyediakan:

  • R2m : nilai kuadrat R marjinal yang terkait dengan efek tetap

  • Nilai R2 bersyarat R2c terkait dengan efek tetap ditambah efek acak.

Catatan: komentar pada posting blog tertaut menunjukkan bahwa pendekatan alternatif yang diilhami Nakagawa & Schielzeth yang dikembangkan oleh Jon Lefcheck (menggunakan sem.model.fitsfungsi dalam piecewiseSEMpaket) menghasilkan hasil yang identik. [Jadi Anda punya opsi: p].

  • Saya tidak menguji fungsi terakhir ini, tapi saya menguji r.squaredGLMM()fungsi dalam MuMInpaket dan dapat membuktikan bahwa itu masih berfungsi hari ini (2018).

  • 2


1: Nakagawa, S., dan Schielzeth, H. 2013. Metode umum dan sederhana untuk mendapatkan R2 dari model efek campuran linier umum. Metode dalam Ekologi dan Evolusi 4 (2): 133-142.

2: Johnson, PCD 2014, ekstensi R2GLMM Nakagawa & Schielzeth ke model lereng acak. Metode dalam Ekologi dan Evolusi 5: 44-946.

ahli kehutanan
sumber
1
Terima kasih @theforestecologist atas jawaban Anda. Saya akan melihat paket-paket yang disebutkan.
Manuel Ramón