ฉันกำลังตรวจสอบงานบางอย่างและได้พบกับสิ่งต่อไปนี้ซึ่งดูเหมือนว่าผิดสำหรับฉัน รุ่นสองแบบผสมถูกติดตั้ง (ใน R) โดยใช้ lmer แบบจำลองนั้นไม่ซ้อนกันและถูกเปรียบเทียบโดยการทดสอบอัตราส่วนความน่าจะเป็น ในระยะสั้นนี่คือตัวอย่างที่ทำซ้ำได้ของสิ่งที่ฉันมี:
set.seed(105)
Resp = rnorm(100)
A = factor(rep(1:5,each=20))
B = factor(rep(1:2,times=50))
C = rep(1:4, times=25)
m1 = lmer(Resp ~ A + (1|C), REML = TRUE)
m2 = lmer(Resp ~ B + (1|C), REML = TRUE)
anova(m1,m2)
เท่าที่ฉันเห็นสามารถlmer
ใช้เพื่อคำนวณความน่าจะเป็นและบันทึกanova
การทดสอบความแตกต่างระหว่างแบบจำลองที่ใช้ไคสแควร์กับองศาอิสระทั่วไป ดูเหมือนจะไม่ถูกต้องสำหรับฉัน ถ้ามันถูกต้องไม่มีใครทราบถึงการอ้างอิงใด ๆ ที่แสดงความชอบธรรมนี้หรือไม่? ฉันตระหนักถึงวิธีการที่ใช้แบบจำลอง (Paper โดย Lewis et al., 2011) และวิธีการที่พัฒนาโดย Vuong (1989) แต่ฉันไม่คิดว่านี่เป็นสิ่งที่ผลิตขึ้นที่นี่ ฉันไม่คิดว่าการใช้anova
คำสั่งนั้นถูกต้อง
anova()
ฟังก์ชั่นใน R ไม่ได้เปรียบเทียบทั้งสองรุ่นที่ติดตั้งภายใต้ REML; มัน refits พวกเขาโดยใช้ ML แล้วทำการทดสอบ ดูซึ่งมีเส้นlme4:::anova.merMod
mods <- lapply(mods, refitML)
(แต่คุณยังคงถูกต้องที่anova()
ไม่สามารถใช้ในการเปรียบเทียบทั้งสองรุ่นเนื่องจากไม่ซ้อนกัน)