aov不會返回p值在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 

無論如何，前兩列是因素，其餘是響應變量。 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值在輸出中報告。