用R做重复测量的方差分析时遇到问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这是个三因素的重复测量方差分析，组间因素是二氧化碳浓度，有两个处理（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这一项的结果

不知道是哪里出了问题，希望大神们为我解答！！万分感谢

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：方差分析 重复测量 treatment Residuals Residual 重复测量方差分析 三因素 R语言

顶贴专用