R执行单因素4水平ANOVA出现contrasts因子报错如何解决

该错误的核心触发逻辑是R的ANOVA拟合要求分组因子至少存在2个有效水平，你的代码逐行计算每个基因的ANOVA p值时，部分基因行满足以下任一条件就会触发报错：

• 预定义的4个分组（is.ctl/II/III/IV）中存在样本量为0的分组：比如dat的列名没有匹配到^NB/^II等正则规则的样本，对应分组的长度为0，生成的分组因子会自动丢弃无样本的水平，最终有效水平数不足2个。
• 单个基因的表达行中多个分组的数值全为NA：aov拟合时会自动删除含NA的行，剔除后仅剩1个分组的有效数据，因子水平不足。
• 单个基因的所有有效表达值无变异（比如全为0、全为同一个常数），也会触发拟合时的水平校验错误。

步骤1：先校验分组样本量 #

运行以下代码确认4个分组都有匹配的样本，如果有任意一个返回值为0，需要先修正grepl的正则匹配规则：

# 输出4个分组的样本数
sum(is.ctl)
sum(II)
sum(III)
sum(IV)

步骤2：修改ANOVA函数增加异常校验 #

对原有函数增加有效水平判断和错误捕获逻辑，不符合计算条件的基因直接返回NA，不会中断整体运行：

aov.all.genes <- function(x,s1,s2,s3,s4) {
  x1 <- as.numeric(x[s1])
  x2 <- as.numeric(x[s2])
  x3 <- as.numeric(x[s3])
  x4 <- as.numeric(x[s4])
  fac <- c(rep("NB",length(x1)), rep("II",length(x2)), rep("III",length(x3)), rep("IV",length(x4)))
  a.dat <- data.frame(factor = as.factor(fac), express = c(x1,x2,x3,x4))
  # 剔除含NA的行
  a.dat <- na.omit(a.dat)
  # 校验有效分组水平数是否≥2
  if(nlevels(droplevels(a.dat$factor)) < 2) {
    return(NA)
  }
  # 捕获其他拟合异常
  p.out <- tryCatch({
    summary(aov(express~factor, a.dat))[[1]][1,5]
  }, error = function(e) {
    return(NA)
  })
  return(p.out)
}

# 运行ANOVA计算
dat.log <- log2(dat)
aov.pv <- apply(dat.log, 1, aov.all.genes,s1=is.ctl,s2=II,s3=III,s4=IV)

运行完成后可自行过滤掉aov.pv中为NA的结果，剩余为有效计算的p值。

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