如何从 R 中的 aov 摘要中提取子集?

How to subset from aov summary in R?

也许这是一个简单的问题,但我想知道如何将 DFF.values 子集用于出现在 aov 摘要中的术语?

例如,使用基础 R 内置数据集 npk,如何提取残差和其他 DFF.values 出现在以下模型的摘要中:

fit <- summary(aov(yield ~ block + N * P + K, data = npk)) # example is fully reproducible

P.S. 我正在寻找基础 R 解决方案。

fit 输出是 length 1 的 list(通过检查 str(fit))。我们用 [[ 提取它,然后 $[[ 提取组件

fit[[1]]$Df
#[1]  5  1  1  1  1 14 #where 14 is the Residuals df
fit[[1]]$`F value`
#[1]  4.391098 12.105541  0.537330  6.088639  1.361073        NA