RM-ANCOVA,但在vcov.default(mod,complete = F)...消息中出错

时间:2019-07-08 12:48:47

标签: r covariance anova

我正在尝试进行重复测量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,重复效应为hemitotal是响应。

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语句。

0 个答案:

没有答案