循环函数中忽略NA值进行异常值检测的问题求助

核心错误修复 #

报错的直接原因是quantile()函数默认不处理NA值，必须显式指定na.rm=TRUE。同时你的代码存在语法和逻辑问题，逐一修正如下：

1. 修正异常值检测函数 #

在quantile()中添加na.rm=TRUE，同时处理NA值在逻辑判断中的返回结果（避免出现NA，将其转为FALSE）：

# 创建异常值检测函数
outlier <- function(x) {
  # 计算四分位数，忽略NA值
  Quantile1 <- quantile(x, probs = .25, na.rm = TRUE)
  Quantile3 <- quantile(x, probs = .75, na.rm = TRUE)
  IQR <- Quantile3 - Quantile1
  # 检测异常值，将NA转为FALSE避免索引错误
  result <- x > Quantile3 + (IQR * 1.5) | x < Quantile1 - (IQR * 1.5)
  result[is.na(result)] <- FALSE
  return(result)
}

2. 修正异常值移除函数 #

• 修复赋值符号的空格错误（pheno_bw < -改为pheno_bw <-）
• 替换未定义的pheno_bw1为函数内的变量pheno_bw
• 优化逻辑：基于原始数据标记所有指定列的异常行后统一删除（避免逐列删除导致的迭代偏差）：

# 创建异常值移除函数
remove_outlier <- function(pheno_bw, columns = names(pheno_bw)) {
  # 初始化异常行标记为FALSE
  is_outlier <- rep(FALSE, nrow(pheno_bw))
  # 遍历目标列，标记所有异常行
  for (col in columns) {
    is_outlier <- is_outlier | outlier(pheno_bw[[col]])
  }
  # 移除异常行
  pheno_bw_clean <- pheno_bw[!is_outlier, ]
  # 输出并返回结果
  print("移除异常值后的数据框：")
  print(pheno_bw_clean)
  return(pheno_bw_clean)
}

3. 调用函数 #

# 调用移除函数，指定目标列
cleaned_data <- remove_outlier(pheno_BW_stat, c('BW_M', 'BW_Pa'))

可选逻辑（逐列移除） #

如果你的需求是逐列移除该列的异常值（每处理一列就删除对应异常行，后续列基于已删除的数据处理），保留原循环逻辑并修正错误即可：

remove_outlier <- function(pheno_bw, columns = names(pheno_bw)) {
  for (col in columns) {
    # 修正赋值符号的空格
    pheno_bw <- pheno_bw[!outlier(pheno_bw[[col]]), ]
  }
  print("移除异常值后的数据框：")
  print(pheno_bw)
  return(pheno_bw)
}

验证说明 #

修正后的代码中，quantile()会自动忽略NA值，异常值检测逻辑也会处理NA，不会再抛出"missing values ​​and NaN are not allowed if 'na.rm' is FALSE"的错误。

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