在模型摘要的模型表的输出中添加自由度(df或df.error)

Question

我正在用模型摘要来拟合模型列表和生成结果表。下面是一个可重复的例子：

library(palmerpenguins)
library(modelsummary)

m_test <- list(
  lm(bill_length_mm ~ flipper_length_mm, data = penguins),
  lm(bill_length_mm ~ flipper_length_mm + body_mass_g, data = penguins))

modelsummary(m_test)

默认的modelsummary表输出估计的系数和带有标准错误的截取，以及观察的数量、AIC、BIC、F、RMSE、R2和R2 Adj.。但是，我也希望这个表为每个模型打印自由度。我试过modelsummary(..., metrics="all",...)和modelsummary(..., metrics=c("R2","AIC","df.error",...) ( DF的各种名称如"df.residual")，但这些都没有用。

Answer 1

解决这一问题的方法之一是此帮助文件用于模型摘要。下面的代码修改模型摘要列表对象，以便打印的内容包括df.error。诀窍是找出所需的变量是df.error，并且它驻留在tidy中，但我能让它打印的唯一方法是将它分配给glance。也许有一个不那么麻烦的解决方案，但这对我有效。

对于一个单一的模型：

m_simple <- modelsummary(list(
  lm(bill_length_mm ~ flipper_length_mm, data = penguins)),
  output = "modelsummary_list"
)

m_simple$glance$df.error <- m_simple$tidy$df.error[1] # IMPORTANT: you should visually inspect this to be sure it's the metric of df you want! 
modelsummary(m_simple)

关于一份模型清单：

m_prep <- modelsummary(list(
  lm(bill_length_mm ~ flipper_length_mm, data = penguins),
  lm(bill_length_mm ~ flipper_length_mm + body_mass_g, data = penguins)
),output = "modelsummary_list" )
  
m_out <- as.list(m_prep)
for(i in 1:length(m_out)){
  m_out[[i]]$glance$df.error <- m_out[[i]]$tidy$df.error[1]
}

modelsummary(m_out) 

注意，如果不亲自检查模型，df.error列表中的第一个元素就是您想要的元素，这很容易打印错误的度量。