R语言多列数据对应不同mu值执行单样本T检验的实现方法

错误原因 #

你触发报错的核心原因有2个：

• 单样本t检验的mu参数仅支持单次传入单个数值，你直接传入了整个muvals向量，不符合函数参数要求
• 你定义的varlist逻辑错误：c(t1, t2, t3)是将三个向量的数值拼接为一个长向量，并非数据集列名引用；同时$运算符不支持直接接变量名取列，需要用[[]]或[]做索引

实现方案 #

基础版（按列顺序匹配mu值，无需额外定义变量列表） #

直接用base R的mapply同步遍历数据集列和mu值，不需要额外安装依赖包：

# 定义函数
onett <- function(tpoint, muvals) {
  # 同步遍历列数据和对应mu值执行检验
  t_res <- mapply(function(col_data, mu) {
    t.test(col_data, mu = mu)
  }, tpoint, muvals, SIMPLIFY = FALSE)
  # 给结果绑定列名，方便后续查询
  names(t_res) <- colnames(tpoint)
  return(t_res)
}

# 调用示例
# 你的测试数据
t1 <- rnorm(20, 10, 1)
t2 <- rnorm(20, 10, 1)
t3 <- rnorm(20, 10, 1)
test_data <- data.frame(t1, t2, t3)
# 你的mu值列表
muvals <- c(24, 51.8, 21.89)
# 执行检验
test_result <- onett(test_data, muvals)
# 查看某一列的检验结果，比如t1的检验输出
test_result$t1

进阶版（结果结构化导出为数据框） #

如果需要把所有检验结果整理成表格，可在上述基础上增加结果转换逻辑：

result_df <- do.call(rbind, lapply(test_result, function(x) {
  data.frame(
    变量名 = names(x$statistic),
    t值 = x$statistic,
    p值 = x$p.value,
    95%置信区间下限 = x$conf.int[1],
    95%置信区间上限 = x$conf.int[2],
    对比mu值 = x$null.value,
    row.names = NULL
  )
}))

自定义变量匹配版（按指定变量名匹配mu值） #

如果你不需要遍历所有列，只需要对指定列执行检验，可修改函数如下：

onett <- function(tpoint, var_names, muvals) {
  t_res <- mapply(function(var, mu) {
    t.test(tpoint[[var]], mu = mu)
  }, var_names, muvals, SIMPLIFY = FALSE)
  names(t_res) <- var_names
  return(t_res)
}

# 调用示例
# 注意变量名要写为字符串格式
varlist <- c("t1", "t2", "t3")
test_result <- onett(test_data, varlist, muvals)

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