如何为数据框批量应用prop.test并添加P值作为新行

先帮你把原始数据框还原出来，方便后续操作：

# 创建示例数据框
TableAvE_Island <- data.frame(
  Type = c("Type1", "Type2"),
  Other = c(10, 19),
  Island = c(4, 45),
  N_Shelf = c(1, 1),
  N_Shore = c(0, 9),
  S_Shore = c(3, 11),
  Sum = c(18, 85)
)

接下来，我们可以用apply()函数批量处理需要计算的列（也就是除了Type和Sum之外的列），步骤如下：

1. 定义一个自定义函数，用来对单列计算prop.test的p值：

get_prop_pvalue <- function(col) {
  # col是当前处理的列，Sum是每行的总和列
  test_result <- prop.test(x = col, n = TableAvE_Island$Sum)
  # 提取p值并返回
  return(test_result$p.value)
}

2. 用apply()遍历目标列，批量计算p值：

# 选择需要处理的列：排除Type和Sum
target_cols <- TableAvE_Island[, !names(TableAvE_Island) %in% c("Type", "Sum")]
# 按列应用函数，得到所有p值的向量
p_values <- apply(target_cols, 2, get_prop_pvalue)

3. 把p值转换成数据框的行，合并到原数据框中：

# 把p值向量转成数据框，设置行名为"p-value"，并补全Type和Sum列（Sum列这里设为NA即可）
p_row <- data.frame(
  Type = "p-value",
  t(p_values),  # 转置向量成列，匹配原数据框结构
  Sum = NA,
  stringsAsFactors = FALSE
)
# 合并原数据框和新行
final_table <- rbind(TableAvE_Island, p_row)

运行完之后，final_table就是包含p值新行的结果啦。

关键细节解释： #

• apply()的第二个参数2表示按列处理（传入1则是按行），第三个参数是我们定义的函数，每一列都会自动传入这个函数完成计算。
• 自定义函数里直接调用prop.test()，然后从返回的结果对象中提取p.value属性，这就是我们需要的核心值。
• 最后用rbind()合并时，要确保新行的列名和原数据框完全一致，所以用t(p_values)把向量转成列结构，再补全Type和Sum列。

如果需要调整置信区间或者其他prop.test的参数，直接在自定义函数里的prop.test()中添加参数就行，比如prop.test(x = col, n = TableAvE_Island$Sum, conf.level = 0.99)。

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