【发布时间】：2014-06-22 20:51:58
【问题描述】：
当响应是向量值时，我应该在 R 中为 ANOVA 模型使用什么代码......即 假设我有 20 个人的纵向数据，每个人都有 10 个时间点的测量值......现在我有一个具有 3 个水平的因子 X，比如 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
【问题讨论】：
-
您遇到的具体问题是什么？ Stack Overflow 不是代码编写服务，因此您需要展示您目前编写的代码或准确说明您遇到问题的部分。请先阅读帮助：stackoverflow.com/help/on-topic
-
在上述样本数据中，Y的每一行都是对特定个体的观察
标签： anova