Какие коды следует использовать в R для модели ANOVA, когда ответ имеет векторное значение... т.е. Предположим, у меня есть лонгитюдные данные для 20 человек, каждый из которых имеет измерения в 10 моментов времени... теперь у меня есть фактор X, имеющий 3 уровня, скажем, 0,1,2... Мне нужно проверить, значительно ли уровни отличаются друг от друга. .. Я должен проверить средний вектор (вектор, поскольку каждый человек содержит 10 моментов времени)..т.е. если средний вектор для уровня 0, средний вектор для уровня 1 и средний вектор для уровня 2 существенно различаются...
Мои образцы данных:
Y
[1,] 9.759608 15.02230 17.70331
[2,] 9.596711 15.50542 18.49343
[3,] 11.298570 17.44781 19.48276
[4,] 8.519376 13.73086 17.05881
[5,] 10.232851 15.85302 19.87476
[6,] 10.888219 16.05568 20.12624
[7,] 9.688724 15.50494 18.82778
[8,] 10.309219 16.78230 18.80428
[9,] 9.620743 15.84582 19.32465
[10,] 10.418802 16.18098 17.94019
>treatment=c(0,1,1,2,0,2,1,1,0,1)
>treatment=factor(treatment)
> result=aov(Y~treatment)
Error in model.frame.default(formula = Y ~ treatment, drop.unused.levels = TRUE) :
object is not a matrix