如何将长格式数据转宽格式并添加基于原列的因子行

问题背景 #

我一直没找到这个问题的解决方案，估计是因为这种需求在数据分析规范、最佳实践和数据存储方式里比较少见。我需要把数据转换成指定格式，方便在Primer里做后续分析。小数据集可以手动添加因子行，但大数据集操作太繁琐，希望能实现自动化处理。

我的核心目标是：将长格式数据转换为宽格式，同时在数据框底部添加基于原数据列的因子行。

原始数据 #

我的初始数据框如下：

data <- data.frame(nr = c(1, 2, 3, 4, 5, 6, 7, 8), 
                 year = c(2013, 2013, 2013, 2013, 2022, 2022, 2022, 2022), 
                 depth = c(35, 35, 50, 50, 35, 35, 50, 50), 
                 species = c("A", "B", "A", "D", "C", "B", "D", "A"), 
                 area = c(1.0, 0.5, 3.2, 4.3, 2.0, 5.6, 1.8, 2.3))

期望输出格式 #

最终需要的输出满足以下要求：

• 首行列名由year、depth以及重复编号nr组合而成
• 中间行是各物种对应的area值的宽格式数据（缺失值填充为0）
• 数据框底部添加三行因子信息：分别对应year、depth以及两者的交互项

当前实现进展 #

我已经完成了宽格式转换和列名设置，但还缺少底部的因子行，代码如下：

primer <- data %>%
  pivot_wider(names_from = c(year, depth, nr), values_from = area) %>% 
  mutate(across(.cols = everything(), ~replace_na(.x, 0))) %>% 
  as.data.frame()

现有思路的不足 #

我曾尝试从列名中提取因子并转为向量，再将year、depth等因子行绑定到primer数据框，但这种方法通用性很差——每次更换因子类型或调整因子数量时，都要修改函数逻辑。因此需要一个更通用的解决方案。

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