R中使用expss添加的标签写入Stata文件后丢失问题求助

• 标签存储格式不匹配：expss::apply_labels添加的变量标签存储在每个变量的独立label属性中，但foreign::write.dta要求标签存储在整个数据框的全局var.labels属性中（是长度和变量数相等的字符向量），所以会出现单独查看变量标签存在，但写入dta时识别不到的情况。
• 变量名缩写警告的影响：foreign::write.dta默认兼容Stata 8及更早版本，变量名最长仅支持32个字符，你的变量LS_medianpovgap60_disp_wa长度超出限制，会被自动截断缩写，进一步导致标签和变量对应错位。

直接替换foreign包为haven包写入dta文件，haven是tidyverse官方维护的Stata数据交互工具，原生兼容expss的标签格式，支持高版本Stata的长变量名，不需要额外调整属性：

# 未安装haven先执行 install.packages("haven")
library(expss)
library(dplyr)
library(haven)

# 原生成数据、添加标签代码不变
df <- data.frame(country = rep(c("NL", "DE", "FR", "AT"), 2),
                 year = rep(c(2012,2014), 4),
                 LS_medianpovgap60_disp_wa = c(0.448257605781815, 0.468249874784546, 0.473270740126805, 0.483814288478694, 0.486781335455043, 0.49246341926957, 0.51121872756711, 0.556027028656306))

df <- apply_labels(df,
                   country = "Country",
                   year = "Year",
                   LS_medianpovgap60_disp_wa = "Median shortfall from the poverty thresholds using 60% of the median income, disposable income only households with working age (LIS and SILC average)")

# 指定Stata版本为14及以上，支持长变量名
write_dta(df, "df_labelled.dta", version = 14)

写入后不管是用read_dta回读R，还是直接在Stata中打开，变量标签、值标签都会正常显示，也不会出现变量名缩写的警告。

需要手动补齐全局var.labels属性，同时把过长的变量名缩短到32字符以内：

# 提取所有变量的label赋值给全局var.labels属性
attr(df, "var.labels") <- sapply(df, function(x) attr(x, "label", exact = TRUE))
# 把超长变量名改短到32字符以内
colnames(df)[colnames(df) == "LS_medianpovgap60_disp_wa"] <- "LS_medianpovgap_disp_wa"
# 再执行写入
write.dta(df, "df_labelled.dta")

注意该方案对值标签的支持较差，如果有数值标签需要额外转换为factor格式才能被识别，因此优先推荐使用haven方案。

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