重复测量数据的分析方法
2017-03-29 本文已影响0人
育种数据分析之放飞自我
数据和思路来自: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间有联系。
结果:AIC 显示832.06
4,Antedependence order2
结果:AIC为837.07
图片.png5,Power model(city-block metric)
结果:AIC842.32
6,Unstructured
结果:AIC833.64
图片.png结论
根据AIC可知,Antedependence order1
数据2,按列排:
数据下载:dat2
模型:
结果:
图片.png