-1

我正在比较两个模型,其中一个使用 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

我的分析可能有问题吗?

4

1 回答 1

2

按 PatientNum 对随机效应进行分组与按 Hospital:patientNum 分组相同,因为 patientNum 是唯一的。

如果你得到table(patientNum)然后table(interaction(hospital, patientNum))你应该注意到所有的单元格都是相同的大小。

于 2013-10-15T15:56:55.387 回答