如何在R中为数据框所有数值变量生成QQ图评估单变量正态性

需求 #

为ggplot2包中的txhousing数据框所有数值变量生成独立QQ图，仅显示图形，用于评估单变量正态性。

问题分析 #

• 使用apply调用EnvStats::qqPlot时的警告：仅为数据中存在NA的提示，属于正常信息，图形仍能生成；若需消除警告，可显式移除NA。
• 原for循环报错'x' must be a numeric vector：因为txhousing是tibble类型，df[,i]返回的是单列tibble而非数值向量，而qqPlot要求输入向量，改用df[[i]]即可解决。

可行方案 #

方案1：使用EnvStats包（带统计增强的QQ图） #

library(ggplot2)
library(EnvStats)

df <- txhousing

# 遍历所有列，仅处理数值型变量
for (col_name in names(df)) {
  col_data <- df[[col_name]]
  if (is.numeric(col_data)) {
    # 生成QQ图，显式移除NA，添加变量名作为标题
    qqPlot(na.omit(col_data), main = paste("QQ Plot for", col_name))
  }
}

说明：na.omit(col_data)提前移除NA，避免警告；df[[col_name]]确保获取数值向量，符合qqPlot的参数要求。

方案2：使用ggplot2包（自定义风格QQ图） #

如果偏好ggplot2的可视化风格，可结合purrr批量生成：

library(ggplot2)
library(purrr)

df <- txhousing
# 筛选所有数值列
num_columns <- df |> dplyr::select(where(is.numeric))

# 遍历数值列生成并显示QQ图
walk(names(num_columns), function(col) {
  plot <- ggplot(num_columns, aes(sample = .data[[col]])) +
    geom_qq(color = "#2c3e50") +
    geom_qq_line(color = "#e74c3c") +
    labs(title = paste("QQ Plot for", col),
         subtitle = "Assessing Univariate Normality") +
    theme_minimal()
  print(plot)
})

说明：walk函数用于执行绘图动作而不存储结果；.data[[col]]正确引用列名，避免ggplot2的非标准评估问题；可自定义图形颜色、主题等样式。

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