我正在尝试进行重复测量ANCOVA,并调查了SE上找到的其他答案,但仍然不觉得我的情况与其他发布者的数据相符。另一个poster已寻求帮助解决此问题,但我想介绍一下自己的数据。
我想对以下数据执行RM-ANCOVA:
df<-read.csv(header=T, text="
ID, Group, Hemi, total, wealth
100, 0, 0, 39, 1.99
100, 0, 1, 44, 1.99
101, 0, 0, 38, 2.01
101, 0, 1, 23, 2.01
200, 1, 0, 33, 80.1
200, 1, 1, 23, 80.1
202, 1, 0, 18, 12.1
202, 1, 1, 28, 12.1")
#hemi=0 for left, hemi=1 for right
df$Group<-as.factor(df$Group)
df$Hemi<-as.factor(df$Hemi)
此数据中的协变量为wealth
,对象间效应为Group
,重复效应为hemi
。 total
是响应。
library('car')
ancova.brain<-aov(total~Hemi+Group+Group*Hemi+Error(ID)+wealth, data=df)
Anova(ancova.brain,type="III")
返回错误:
Error in vcov.default(mod, complete = FALSE) :
there is no vcov() method for models of class aovlist, listof
有什么办法可以避免此错误?我认为问题在于Error
语句中不能使用aov
语句。