这是个三因素的重复测量方差分析,组间因素是二氧化碳浓度,有两个处理(Treatment);另外两个因素是年份(Year)和平均海平面高度(MSL),观测变量是 生物量(SM),总共有10个生长箱来经行这些实验(Chamber),我认为Year和MSL都是重复测量因素,因此我的代码是:
model<-aov(SM~Treatment*MSL*Year+Error(Chamber/(MSL+Year)),data=DATA)
summary(model)
但提示错误:
Warning message:
In aov(SM ~ Treatment * MSL * Year + Error(Chamber/(MSL + Year)), :
Error() model is singular
并且结果是这样的:
Error: Chamber
Df Sum Sq Mean Sq F value Pr(>F)
Treatment 1 58009 58009 0.245 0.634
Residuals 8 1891035 236379
Error: Chamber:MSL
Df Sum Sq Mean Sq F value Pr(>F)
MSL 5 448928 89786 15.556 1.71e-08 ***
Treatment:MSL 5 30260 6052 1.049 0.403
Residuals 40 230865 5772
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
没有其他的交互项,并且我尝试了不做重复测量的普通三因素方差分析:model<-aov(SM~Treatment*MSL*Year),data=DATA)结果任然没有含有Year这一项的结果
不知道是哪里出了问题,希望大神们为我解答!!万分感谢