Base R实现：展开聚合操作后嵌套在数据框列中的多返回值矩阵

我完全懂你碰到的这个麻烦——用aggregate()做聚合时，如果传入的函数（比如这里的quantile()）返回多个值，结果会把这些值打包成一个矩阵嵌套在数据框的单个列里。虽然打印出来看起来像是多列，但实际数据框只有g和x两列；要是直接用do.call(cbind, x.a)，又会把分组的因子列强制转成数值，这显然不是我们想要的效果。

下面给你一个纯Base R的解决方案，既能展开嵌套的矩阵列，又能保留分组列的原始类型：

步骤1：准备示例数据（和你的MWE一致） #

n = 10
df = data.frame(
  x = runif(n),
  g = sample(c('A','B'), n, rep = TRUE)
)
# 聚合得到带嵌套矩阵的结果
x.a = aggregate(x ~ g, df, quantile, probs = c(.025, .5, .975))

步骤2：展开嵌套矩阵并合并分组列 #

# 提取分组列，保留原有的因子类型
group_data <- x.a["g"]

# 将嵌套的矩阵列转换为数据框，自动继承分位数的列名
expanded_x <- as.data.frame(x.a$x)
# 如果想给列名加上x前缀（和你期望的一致），可以手动修改
colnames(expanded_x) <- paste0("x.", colnames(expanded_x))

# 合并分组列和展开后的数值列
final_result <- cbind(group_data, expanded_x)

验证结果 #

运行后查看结果和列名：

print(final_result)
print(names(final_result))

输出的列名会是"g" "x.2.5%" "x.50%" "x.97.5%"，而且final_result$g依然是因子类型，不会被转成数值。

原理说明 #

这个方法的关键是单独处理分组列和嵌套矩阵列：先把分组列单独提取出来（保持数据框格式，避免因子被强制转换），再把嵌套的矩阵转成独立的数据框，最后用cbind()合并两者。这样既保留了分组列的类型，又把矩阵里的多值拆分成了单独的列。

备注：内容来源于stack exchange，提问作者jessexknight