小鼠插补后比较glmer与glm

时间:2019-06-09 07:42:05

标签: statistics imputation

我想在多次插补小鼠后测试我的随机效应术语的意义。我有两个嵌套模型,我试图与池比较进行比较,但是这会产生一些错误。

fm1 <-with(dti.mice1,glmer(治疗〜(1 |医院)+年龄))

fm2 <-与(dti.mice1,glm(治疗〜年龄))

pool.compare(fm1,fm2)

错误:模型'fit1'不大于'fit0'

1 个答案:

答案 0 :(得分:0)

问题似乎是,fm1和fm2无法识别为嵌套模型 (其中fm1是更大的型号)。

小鼠文档中:

  

fit1
  由with.mids()产生的'mira'类对象。

     

fit0
  with.mids()产生的'mira'类对象。 fit0中的模型是fit1的嵌套fit0。

因此,完整/较大的模型必须适合1。您的示例就是这种情况。 可能是因为您无法同时使用glm和glmer。如果两个模型都使用glm,会发生什么情况?这行得通吗?