使用plm::pdata.frame后labelled变量标签在msummary中丢失的解决方法

方法1：提取原变量标签并手动映射给msummary #

plm::pdata.frame()转换后，labelled包设置的标签通常还保存在变量的label属性中，你可以提取这些标签，再通过modelsummary的参数指定显示内容：

1. 提取变量标签：

# 假设你的pdata.frame对象名为p_df
var_labels <- sapply(p_df, function(x) attr(x, "label"))
# 过滤无标签的变量
var_labels <- var_labels[!is.na(var_labels)]

2. 在msummary()中用coef_rename参数映射标签：

modelsummary::msummary(your_fitted_model, coef_rename = var_labels)

如果模型包含滞后项（如lag(var)），需要单独给滞后项补充标签：

var_labels["lag(var)"] <- "滞后的[原变量标签]"

方法2：转换后重新设置变量标签 #

部分场景下转换会丢失标签属性，直接重新用labelled::var_label()赋值即可：

labelled::var_label(p_df$var1) <- "变量1的自定义标签"
labelled::var_label(p_df$var2) <- "变量2的自定义标签"

完成后重新拟合模型再调用msummary()，标签一般会被正常识别。

方法3：用modelsummary钩子函数自定义标签提取逻辑 #

如果前两种方法失效，可以通过钩子函数强制匹配标签：

# 定义自定义提取函数
get_estimates_custom <- function(model, ...) {
  est <- modelsummary::get_estimates(model, ...)
  # 从模型的原始数据中提取标签
  df <- model$model
  labels <- sapply(df, function(x) attr(x, "label"))
  # 替换变量名为对应标签
  est$term <- ifelse(est$term %in% names(labels), labels[est$term], est$term)
  return(est)
}

# 调用msummary时指定钩子
modelsummary::msummary(your_fitted_model, get_estimates = get_estimates_custom)

注意：如果模型包含交互项、多项式项等复杂形式，需要调整钩子函数中的匹配逻辑，确保标签能正确对应。

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