aov不会在R

Question

我对anova摘要结果summary(aov)有一个奇怪的问题。

所以这就是问题所在。我有一个包含6列的数据集。以下是数据集示例：

Panelist    Prod.ID Overall Appearance  Flavor  Texture
1   196 9   9   9   9
1   239 7   9   6   7
1   354 9   8   8   7
1   427 3   8   2   3
1   577 8   9   7   9
1   638 7   9   7   8
1   772 6   4   3   3
1   852 9   8   9   8
2   196 8   8   7   8
2   239 7   7   7   7
2   354 6   5   6   4
2   427 6   7   4   6
2   577 3   6   3   5
2   638 4   4   5   2
2   772 6   2   6   7
2   852 7   6   7   6
3   196 7   9   7   8
3   239 8   9   8   8
3   354 8   8   7   8
3   427 7   8   6   8
3   577 8   9   8   8
3   638 8   9   8   7
3   772 5   8   8   8
3   852 8   9   8   8

无论如何，前两列是因子，其余是响应变量。小组成员和Prod.ID被summary()视为连续变量，因此我将它们转换为as.factor()的因子。

在转换之后，我使用以下模型Overall ~ Panelist * Prod.ID运行anova-test，但作为汇总结果我只得到了：

> summary(aov(Overall ~ Prod.ID * Panelist, data = paneElements))

                   Df Sum Sq Mean Sq
Prod.ID             7  189.6  27.085
Panelist          160 1252.9   7.830
Prod.ID:Panelist 1116 3116.1   2.792

我无法找到导致F测试值和P值消失的任何原因。

非常感谢任何帮助。

非常感谢。

Answer 1

Prod.ID和Panelist的每个组合只有一个观察值（至少在您的样本数据中），因此组的数量等于观察数量。这将导致F检验中的除零，这可能是缺乏报告的F检验和p值的原因。

例如，当我为Prod.ID的{​​{1}}添加额外观察时，我会在输出中报告F和p值。