R语言run_limma函数调用报错：对象contr_1未找到求助

错误原因 #

调用函数时传入的contr_1未加引号，R会将其视为需要预先定义的对象，但你并未创建过名为contr_1的对象，因此触发"对象'contr_1'未找到"的错误。函数逻辑本身没问题，只需传入字符串类型参数，但调用时的语法错误导致R解析失败。

直接解决方法 #

调用函数时，将contr_1用引号包裹，明确指定这是字符串参数：

run_limma(model_1, support_1, fit_1, editing_1, "contr_1", tmp_1)

增强型解决方案（添加参数校验） #

在函数开头加入参数类型检查，避免后续再次出现同类错误：

run_limma <- function (model_x, support_x, fit_x, editing_x, contr_x, tmp_x){
  # 新增参数校验，确保传入的是字符串
  if (!is.character(contr_x)) {
    stop("contr_x必须为字符串类型，请传入带引号的参数，例如\"contr_1\"")
  }
  
  message("starting modeling")
  model_x <- model(support_x, model_x)
  message("starting fitting")
  fit_x <- limma_diff(editing_x, model_x, fit_x)
  message("Making contrasts")
  
  if (contr_x == "contr_1") {
    contr_x <- makeContrasts(diseaseAD - diseaseControl, levels = colnames(coef(fit_x)))
  }
  if (contr_x == "contr_2") {
    contr_x <- makeContrasts(diseaseAD_MCI - diseaseControl, levels = colnames(coef(fit_x)))
  }
  if (contr_x == "contr_3") {
    contr_x <- makeContrasts(diseaseMCI - diseaseControl, levels = colnames(coef(fit_x)))
  }
  if (contr_x == "contr_4") {
    contr_x <- makeContrasts(diseasePD - diseaseControl, levels = colnames(coef(fit_x)))
  }
  
  message("making tmp file")
  tmp_x <- limma_cont(contr_x, fit_x, tmp_x)
  tmp_x
}

代码优化建议 #

用switch语句替代多个if判断，让代码更简洁易读，同时增加无效参数的处理：

run_limma <- function (model_x, support_x, fit_x, editing_x, contr_x, tmp_x){
  if (!is.character(contr_x)) {
    stop("contr_x必须为字符串类型，请传入带引号的参数，例如\"contr_1\"")
  }
  
  message("starting modeling")
  model_x <- model(support_x, model_x)
  message("starting fitting")
  fit_x <- limma_diff(editing_x, model_x, fit_x)
  message("Making contrasts")
  
  # 使用switch替代多分支if
  contr_x <- switch(contr_x,
    "contr_1" = makeContrasts(diseaseAD - diseaseControl, levels = colnames(coef(fit_x))),
    "contr_2" = makeContrasts(diseaseAD_MCI - diseaseControl, levels = colnames(coef(fit_x))),
    "contr_3" = makeContrasts(diseaseMCI - diseaseControl, levels = colnames(coef(fit_x))),
    "contr_4" = makeContrasts(diseasePD - diseaseControl, levels = colnames(coef(fit_x))),
    # 处理未匹配的参数值
    stop("无效的contr_x参数，可选值为：contr_1、contr_2、contr_3、contr_4")
  )
  
  message("making tmp file")
  tmp_x <- limma_cont(contr_x, fit_x, tmp_x)
  tmp_x
}

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