Я пытаюсь выполнить ANCOVA с повторными измерениями и изучил другие ответы, найденные на SE, но все еще не чувствую, что мой случай соответствует данным других авторов.Другой постер обратился за помощью по этой проблеме, но я хотел бы представить свои собственные данные, которые немного отличаются.
Я хотел бы выполнить 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.