使用car包Type III ANOVA分析非平衡设计时如何提取分母自由度

以下是两种可行的提取方法，均基于你当前使用的普通线性模型（lm）场景：

• 方法1：从Anova返回结果中直接提取
  你可以先将Anova的输出保存为对象，再索引对应字段取值：

# 保存Anova输出结果到对象
anova_res <- Anova(mod, type = "III")
# 提取所有固定效应共用的分母自由度（即残差自由度）
den_df <- anova_res$Df[nrow(anova_res)]

如果需要查看返回结果的完整可提取字段，可以运行str(anova_res)确认结构。

• 方法2：直接从lm拟合对象提取
  普通线性模型的固定效应F检验分母自由度等价于模型残差自由度，可以直接从你已经拟合好的mod对象中调取，无需重复运行Anova：

den_df <- mod$df.residual

如果需要批量提取单个效应的全部检验结果，可以将Anova输出转为数据框后索引对应行：

anova_df <- as.data.frame(anova_res)
# 示例：提取Partner.type主效应的检验参数
effect_row <- anova_df["Partner.type", ]
num_df <- effect_row$Df # 分子自由度
f_val <- effect_row$`F value` # F值
p_val <- effect_row$`Pr(>F)` # p值

注：上述方案仅适用于当前普通固定效应线性模型（lm）场景，如果后续使用混合效应模型（如lmer），分母自由度会根据不同校正方法（Satterthwaite/Kenward-Roger）发生变化，提取逻辑存在差异。

内容的提问来源于stack exchange，提问作者user310106