方差分析的自由度不正确

Question

我的方差分析结果如下，

mnn=aov(gy~gen+rep, data=my)
summary(mnn)
            Df Sum Sq Mean Sq F value   Pr(>F)    
gen          6 1923.1   320.5  31.206 5.42e-07 ***
rep          1    7.1     7.1   0.695    0.419    
Residuals   13  133.5    10.3                     

为什么当 rep 的级别数是 3在数据.

这是我的示例数据集，

   gen rep gy  sy
1    a   1 40  95
2    b   1 50 120
3    c   1 55 120
4    d   1 60 140
5    e   1 40 110
6    f   1 50 125
7    g   1 65 145
8    a   2 35 100
9    b   2 50 125
10   c   2 59 130
11   d   2 65 150
12   e   2 40 110
13   f   2 55 130
14   g   2 60 145
15   a   3 40 100
16   b   3 50 120
17   c   3 50 130
18   d   3 65 145
19   e   3 40 115
20   f   3 55 130
21   g   3 70 155

Answer 1

您可能打算将 rep 视为 分类 预测变量。因为数据集中 rep 的值是数字 (1,2,3)，R 将假设 rep 是一个连续预测变量（因此您正在做一种形式的 ANCOVA 而不是两个方差分析）。在这种情况下，您必须明确指定变量是分类变量，方法是在数据集中使用 factor() 对其进行转换：

my$rep <- factor(my$rep)

或在您的 aov() 通话期间：

aov(gy~gen+factor(rep), data=my)