aov не вернет pvalues ​​в 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().

После этого обращения я побежал ANOVA-тест с последующей моделью 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-Test, что может быть причиной отсутствия сообщенных F-тестов и p-значений.

Например, когда я добавляю дополнительное наблюдение для Prod.ID 196 только на один уровень Panelist, я получаю значения F и p, указанные в выводе.