如何用t检验比较对照组与各药物组细胞活力？循环代码失效求助

首先咱们先把你的示例数据集整理成R能正确识别的结构——看起来你给出的原始数据格式有点混乱，先把它转换成标准的数据框：

# 构建你的示例数据集，行名为分组名称
dataset <- data.frame(
  sample_01 = c(100, 50, 90, 40),
  sample02 = c(120, 49, 87, 43),
  sample03 = c(110, 48, 80, 42),
  row.names = c("untreated", "drug1", "drug2", "drug3")
)

接下来咱们看看你原来的循环代码出了哪些问题：

• 结果被覆盖：每次循环里result <- t.test(x,y)都会把之前的检验结果覆盖，最后你只能得到最后一组的结果，前面的都丢了。
• 包含了不必要的组：循环从第1行开始，而第1行就是untreated组，这会导致你拿未处理组和自己做t检验，完全没必要。
• 数据类型隐患：如果直接用dataset["untreated",]提取行，得到的是一个数据框而不是数值向量，可能会让t.test函数报错或者处理异常。

下面是修正后的代码，完美实现你要的功能：

# 先提取未处理组的数值，转换成向量确保t.test能正确处理
untreated_vals <- as.numeric(dataset["untreated", ])

# 创建一个空列表来存储每个药物的检验结果，避免结果被覆盖
t_test_results <- list()

# 循环遍历所有药物组（排除untreated行）
for (drug_name in setdiff(rownames(dataset), "untreated")) {
  # 提取当前药物组的数值向量
  drug_vals <- as.numeric(dataset[drug_name, ])
  # 执行独立样本t检验（根据你的实验设计，这里假设两组独立，你可以按需调整参数，比如var.equal等）
  test_result <- t.test(untreated_vals, drug_vals)
  # 把结果存入列表，用药物名作为键，方便后续查找
  t_test_results[[drug_name]] <- test_result
}

现在你可以轻松查看每个药物的检验结果了，比如查看drug1的结果：

t_test_results[["drug1"]]

或者一次性打印所有结果：

for (name in names(t_test_results)) {
  cat("=== 药物", name, "的t检验结果 ===\n")
  print(t_test_results[[name]])
  cat("\n")
}

如果你想让代码更简洁，符合R的函数式编程风格，可以用lapply替代for循环：

t_test_results_lapply <- lapply(setdiff(rownames(dataset), "untreated"), function(drug) {
  t.test(as.numeric(dataset["untreated", ]), as.numeric(dataset[drug, ]))
})
# 给列表命名，方便识别
names(t_test_results_lapply) <- setdiff(rownames(dataset), "untreated")

这样就能完美完成你要的未处理组和每个药物组的t检验对比啦！

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