我正在比较两个模型,其中一个使用2级分析,而另一个使用3级分析,当使用anova时,我得到两个模型的相同结果:使用R,lme4
anova(finalModel, finalmodelML)
Data: mydata
Models:
finalModel: healthState ~ days + gender + age + diabetes + stroke + cardiacRehab +
finalModel: drugAdherence + drugCompliance + (days | patientNum)
finalmodelML: healthState ~ days + gender + age + diabetes + stroke + cardiacRehab +
finalmodelML: drugAdherence + drugCompliance + (days | hospital:patientNum)
Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
finalModel 14 -7259.4 -7167 3643.7 -7287.4
finalmodelML 14 -7259.4 -7167 3643.7 -7287.4 0 0 1
我的分析可能有问题吗?
答案 0 :(得分:2)
按PatientNum对随机效果进行分组将与医院分组相同:patientNum,因为patientNum是唯一的。
如果您获得table(patientNum)
和table(interaction(hospital, patientNum))
,则应注意所有单元格大小相同。