使用方差分析功能拟合ANCOVA

时间:2019-06-18 23:44:28

标签: r anova

我调查了在SO上找到的其他答案,但仍然觉得我的情况与其他海报的数据不符。我想对以下数据执行重复测量-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

此数据中的covariatewealth,具有对象间效应Group和重复效应hemitotal是响应。

#Fitting the ANCOVA 
brain.out<-aov(total~group+hemi+Error(ID)+wealth, data=df)
summary(brain.out)

我的问题是:(1)这是我要建模的正确语法吗? (2)如何在aov中执行带错误的ANCOVA?我知道Anova :: car提供了ANCOVA,但在aov语句中不允许Error

感谢您的帮助。谢谢。

0 个答案:

没有答案