使用reshape2熔解数据遇阻：16个英国家统计局合并表堆叠求助

嘿，我懂你现在的困扰——用reshape2处理英国国家统计局那16张企业规模表格，搞定了一个就卡壳了是吧？别慌，咱一步步把所有表格都规整成想要的长格式数据，保证顺畅完成熔解操作。

因为你每张表格对应一类企业规模（大型、中型…），熔解前必须给每张表格加一列用来标记它的规模类型——不然等下所有数据堆在一起，你根本分不清哪行是大型哪行是中型。

比如读入的大型企业表格叫large_firms，直接加列：

large_firms$firm_size <- "大型企业"

同理，中型企业表格medium_firms就加：

medium_firms$firm_size <- "中型企业"

把16张表格都按这个逻辑加好标识列，别漏掉任何一张。

既然每张表格都是「16个行业+6年数据」的结构，必须保证所有表格的列名完全一致——比如行业列都叫industry，年份列都叫2017/2018…（按你实际的年份来）。如果有表格列名不一样，比如有的叫Sector，先统一：

# 把某表格的Sector列改成industry
colnames(medium_firms)[colnames(medium_firms) == "Sector"] <- "industry"

这步千万别偷懒，列名不统一的话，后面合并直接报错。

现在可以把加了标识列、列名统一的表格全合并到一个数据框里。如果表格少，直接用rbind()：

all_firms <- rbind(large_firms, medium_firms, small_firms, ...) # 把16张表格依次列进去

要是表格多嫌麻烦，就把它们放进列表批量处理，效率更高：

# 先把所有表格放进一个列表
firm_tables <- list(large_firms, medium_firms, small_firms, ...)
# 批量给每个表格加规模标识
firm_tables <- mapply(function(df, size_label) {
  df$firm_size <- size_label
  return(df)
}, firm_tables, c("大型企业", "中型企业", "小型企业", ...), SIMPLIFY = FALSE)
# 一键合并所有表格
all_firms <- do.call(rbind, firm_tables)

现在终于可以用reshape2的melt()函数把宽格式的年份列转成规整的长格式了，指定id.vars为我们的标识列和行业列就行：

library(reshape2)
melted_data <- melt(all_firms, 
                    id.vars = c("industry", "firm_size"), # 保留不熔解的列
                    variable.name = "year", # 转后的年份列名
                    value.name = "firm_count") # 转后的企业数量列名

运行完你就能得到一行对应「行业+企业规模+年份+企业数量」的干净数据了。

如果之前操作失败，大概率是这几个原因：

• 没加firm_size标识列，导致熔解后数据混乱，分不清规模
• 不同表格列名不一致，rbind合并时直接报错
• melt的id.vars参数设错，漏了行业或标识列

给你整个极简版的模拟示例，方便你对照：

# 模拟大型企业数据
large <- data.frame(industry = c("制造业", "零售业"), 
                    "2017" = c(100, 80), "2018" = c(105, 85))
large$firm_size <- "大型企业"

# 模拟中型企业数据
medium <- data.frame(industry = c("制造业", "零售业"), 
                     "2017" = c(200, 150), "2018" = c(210, 155))
medium$firm_size <- "中型企业"

# 合并+熔解
all_firms <- rbind(large, medium)
melted_data <- melt(all_firms, id.vars = c("industry", "firm_size"), 
                    variable.name = "year", value.name = "firm_count")

看看输出的melted_data，是不是就是你想要的结构？

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