2013-10-15 34 views
-1

我比較兩種模式其中一個使用2級分析,而其他用途3級的分析和我得到相同的結果,兩種型號使用ANOVA時:使用R,lme4使用anova時,相同的模型適合結果?

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

我的R文件可能有問題嗎?是否有可能得到1的pr(> chisq) –

+0

我只是過於謹慎......這可能是anova的好結果嗎? –

+0

它看起來像你的第二個模型是無視'醫院'。如果這是真的,我猜想finalModel和finalmodelNL是相同的,因此相同的擬合。沒有關於數據的更多信息,很難多說。 – Chris

回答

2

按照patientNum對您的隨機效果進行分組將與由醫院分組相同:patientNum,因爲patientNum是唯一的。

如果你得到table(patientNum)table(interaction(hospital, patientNum))那麼你應該注意到所有的單元格都是相同的大小。

+0

感謝您的回答 –

相關問題