數據和思路來自:Mixed Models for Missing Data With Repeated Measures
數據:dat
這里Group1為處理,Group2位對照,Time0,1,3,6表示月份時間
1,GenStat重復測量方差分析
圖片.png
結果:
圖片.png
原網頁處理的結果:
圖片.png
方差分析重復測量數據的弱點
- 1,當想定義不同的方差協方差結構是,GLM(一般線性模型)不能實現,只能是LMM(混合線性模型)實現。
- 2,使用方差分析時,我們假定數據滿足球形檢驗(sphericity assumption),如果不滿足我們可以通過一些方法( Greenhouse and Geisser and Hyunh and Feldt)矯正,但是這不是最優選擇。
- 3,如果有缺失值,缺失值的數據在ANOVA重復測量分析時也對結果影響較大。
2,GenStat 中Antedepence Analysis模型
模型:
結果:
3,Antedependence order1
允許subject間有聯系。
圖片.png
結果:AIC 顯示832.06
4,Antedependence order2
結果:AIC為837.07
圖片.png
5,Power model(city-block metric)
結果:AIC842.32
6,Unstructured
結果:AIC833.64
圖片.png
結論
根據AIC可知,Antedependence order1
數據2,按列排:
數據下載:dat2
模型:
結果:
圖片.png