У нас есть набор данных о 3 методах диеты и похудания. Мы хотим запустить тест Anova, чтобы сравнить ЕД диеты номер 1 и диеты 2, и сравнить диеты 1 и 2 с 3 (используя контрасты, чтобы увидеть разницу в эффекте).
Название диеты наши данные: diet_data
Итак, мы попытались запустить это:
res.aov <- aov((dif_vector) ~ factor(diet_d), data = diet_data)
summary(res.aov)
ok2<- as.factor(diet_d)
c1 <- c(-1,1,0)
c2 <- c(-1,-1,2)
mat <- cbind(c1,c2)
contrasts(ok2) <- mat
levels(ok2)
summary.aov(res.aov,split=list(ok2=list("1 vs.2"=1, "3 vs 1,2"=2)))
и получили сообщение об ошибке:
Error in summary.aov(res.aov, split = list(ok2 = list(`1 vs.2` = 1, `3 vs 1,2` = 2))) :
unknown name ‘ok2’ in the 'split' list
Do you есть идеи, в чем проблема? Мы фактически скопировали этот код у кого-то, кому это удалось, но для нас это не сработает.
Это проблема, что три метода диеты не имеют одинакового размера выборки?