Добавление столбца к выводу `summary (aov)` в R?

Question

Мне было интересно, можно ли добавить заголовок столбца к выводу, заданному summary(aov)?

Например, после следующего summary(aov) мы можем добавить столбец с названием aa:

aa = c(4.391098, 12.105541,  0.537330,  6.088639,  1.361073)

# add `aa` above to the right-most column of summary output below:

fit <- summary(aov(yield ~ block + N * P + K, data = npk))

Я безуспешно пробовал следующее:

fit[[2]] <- aa

Answer 1

Если вы хотите печатать со значками звезд, вы не можете добавить новый столбец в качестве последнего столбца. Последний столбец должен быть p-значениями:

aa = c(4.391098, 12.105541,  0.537330,  6.088639,  1.361073)

# add `aa` above to the right-most column of summary output below:

fit <- summary(aov(yield ~ block + N * P + K, data = npk))

fit[[1]]$aa <- c(aa, NA) 
fit[[1]] <- fit[[1]][c(1, 2, 3, 6, 4, 5)]

fit
#            Df Sum Sq Mean Sq     aa F value  Pr(>F)   
#block        5  343.3   68.66  4.391   4.391 0.01295 * 
#N            1  189.3  189.28 12.106  12.106 0.00368 **
#P            1    8.4    8.40  0.537   0.537 0.47564   
#K            1   95.2   95.20  6.089   6.089 0.02711 * 
#N:P          1   21.3   21.28  1.361   1.361 0.26284   
#Residuals   14  218.9   15.64                          
#---
#Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Answer 2

Мы можем использовать cbind для обновления извлеченного list элемента

fit[[1]] <- cbind(fit[[1]],aa = c(aa, NA)) 

Но это удалит класс anova, и это будет data.frame