aov 不会 return R 中的 pvalues

Question

我对方差分析汇总结果有一个 st运行ge 问题 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

无论如何，前两列是因素，其余是响应变量。 Panelist 和 Prod.ID 被 summary() 认为是连续变量，所以我将它们转换为具有 as.factor().

的因子

在那次转换之后，我运行使用以下模型 Overall ~ Panelist * Prod.ID 进行方差分析测试，但作为摘要结果，我只得到了这个：

> 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 值的原因。

例如，当我为 Panelist 的一个水平添加 Prod.ID 196 的额外观察时，我在输出中报告了 F 和 p 值。

aov 不会 return R 中的 pvalues

aov won't return pvalues in R

r

anova

p-value