我正在比较两个模型,其中一个使用2级分析,而另一个使用3级分析,使用方差分析时,两个模型的结果相同:使用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
我的分析会不会有问题?
按 patientNum 对随机效应进行分组与按 hospital:patientNum 分组相同,因为 patientNum 是唯一的。
如果您得到table(patientNum)
和table(interaction(hospital, patientNum))
那么您应该注意到所有单元格的大小都相同。