如何让glht函数及summary输出其计算p值时使用的自由度？

有两种实用方法可以获取或显示glht计算p值时用到的自由度：

1. 直接提取summary对象中的自由度 #

先运行summary()得到结果对象，再提取其中的df字段：

# 示例代码
library(multcomp)
# 拟合线性模型
fit <- lm(mpg ~ factor(cyl) + hp, data = mtcars)
# 生成glht对象
glht_obj <- glht(fit, linfct = mcp(cyl = "Tukey"))
# 生成summary对象
sum_glht <- summary(glht_obj)
# 提取自由度
sum_glht$df

运行后会直接返回计算p值时使用的自由度（比如示例中是28）。

2. 自定义打印函数让summary直接显示自由度 #

如果想让summary输出时自动打印自由度，可以写一个简单的自定义打印函数：

# 自定义打印函数
print.summary.glht_with_df <- function(x, ...) {
  # 调用原打印函数输出常规结果
  print.summary.glht(x, ...)
  # 追加打印自由度信息
  cat("\n计算p值所用自由度：", x$df, "\n")
}

# 使用自定义函数打印summary结果
sum_glht <- summary(glht_obj)
print.summary.glht_with_df(sum_glht)

这样输出的结果会在常规summary内容后，额外显示自由度信息。

补充说明 #

• glht默认会自动从拟合的模型（比如lm模型）中提取残差自由度用于p值计算，不需要手动指定；如果需要手动指定自由度，可以在summary()中添加df参数，比如summary(glht_obj, df = 25)。
• 不同类型的模型（比如广义线性模型），glht使用的自由度逻辑会有差异，此时sum_glht$df也会返回对应场景下的自由度值。

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