我想在多次插补小鼠后测试我的随机效应术语的意义。我有两个嵌套模型,我试图与池比较进行比较,但是这会产生一些错误。
fm1 <-with(dti.mice1,glmer(治疗〜(1 |医院)+年龄))
fm2 <-与(dti.mice1,glm(治疗〜年龄))
pool.compare(fm1,fm2)
错误:模型'fit1'不大于'fit0'
答案 0 :(得分:0)
问题似乎是,fm1和fm2无法识别为嵌套模型 (其中fm1是更大的型号)。
从小鼠文档中:
fit1
由with.mids()产生的'mira'类对象。fit0
with.mids()产生的'mira'类对象。 fit0中的模型是fit1的嵌套fit0。
因此,完整/较大的模型必须适合1。您的示例就是这种情况。 可能是因为您无法同时使用glm和glmer。如果两个模型都使用glm,会发生什么情况?这行得通吗?