如何将数据框列表转换为SIBER包extractPosteriorMeans可用格式

我正在使用R包SIBER中的extractPosteriorMeans函数提取贝叶斯后验分布的均值，该函数要求输入由siberMVN生成的特定格式列表（示例为double_list）。但由于数据集限制，我只能按组-群落组合分开处理，得到数据框列表df_list，当前格式无法适配extractPosteriorMeans。处于分析流程后期，不想重新处理数据集，希望将df_list转换为与double_list一致的格式。

原始代码 #

library(SIBER)

# 从数据框创建自定义列表
df <- structure(list(Sigma2.1.1. = c(0.313422199359486, 0.367808733163473, 
0.287379338451926, 0.266462306294853), Sigma2.2.1. = c(0.543763352968645, 
0.729034607855539, 0.188017384643087, 0.10580503578473), Sigma2.1.2. = c(0.543763352968645, 
0.729034607855539, 0.188017384643087, 0.10580503578473), Sigma2.2.2. = c(1.90928916482304, 
2.63224559339719, 0.318783548730578, 0.25108375378826), mu.1. = c(-21.8944454084601, 
-21.6319111120951, -21.3223533707681, -21.2709073118823), mu.2. = c(4.72633431332481, 
4.80720104049834, 4.47663667743955, 4.61790878072818), group_comm_names = c("2019_claws.ADULT", 
"2019_claws.ADULT", "2019_claws.JU", "2019_claws.JU")), row.names = c(45001L, 
45002L, 54001L, 54002L), class = "data.frame")

df_list <- split(df[setdiff(names(df), "group_comm_names")], f = df$group_comm_names, drop = TRUE)
# 此列表格式无法适配extractPosteriorMeans

# 由siberMVN生成的正确格式示例（目标格式）
double_list <- list(`2019_feathers.ADULT` = structure(c(0.321658761453638, 0.195543225238878, 
0.583801586825703, 0.26254908714109, 0.583801586825703, 0.26254908714109, 
1.89508038050268, 1.68514983009331, -21.8041135046668, -21.6033433747941, 
4.70782882682883, 4.76493678037111), dim = c(2L, 6L), dimnames = list(
    NULL, c("Sigma2[1,1]", "Sigma2[2,1]", "Sigma2[1,2]", "Sigma2[2,2]", 
    "mu[1]", "mu[2]"))), `2021_feathers.ADULT` = structure(c(0.423467729835709, 
0.450662348987855, 0.145852582124701, 0.142899300399445, 0.145852582124701, 
0.142899300399445, 0.246099001638212, 0.34807362924927, -21.0033901314136, 
-20.980878158095, 4.60002257083648, 4.85210622120396), dim = c(2L, 
6L), dimnames = list(NULL, c("Sigma2[1,1]", "Sigma2[2,1]", "Sigma2[1,2]", 
"Sigma2[2,2]", "mu[1]", "mu[2]"))))

解决方案：转换df_list为目标格式 #

要将df_list转换成double_list的格式，需完成以下操作：

• 对每个数据框转置（目标格式行对应后验样本，列对应参数）
• 重命名列名，匹配double_list的参数命名规则
• 将转置后的矩阵转换为带维度名称的结构

执行代码如下：

# 定义转换函数
convert_to_siber_format <- function(df) {
  # 转置数据框为矩阵，调整行列对应关系
  mat <- t(as.matrix(df))
  # 替换列名，匹配目标格式的参数命名
  colnames(mat) <- c("Sigma2[1,1]", "Sigma2[2,1]", "Sigma2[1,2]", "Sigma2[2,2]", "mu[1]", "mu[2]")
  return(mat)
}

# 应用转换函数到整个列表
formatted_list <- lapply(df_list, convert_to_siber_format)

# 验证格式是否匹配目标结构
str(formatted_list)
str(double_list)

验证与使用 #

转换完成后，formatted_list的结构与double_list完全一致，可直接传入extractPosteriorMeans函数：

posterior_means <- extractPosteriorMeans(formatted_list)

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