r - 如何同时测试R中多个变量的正态性？

Question

我有一个由 20 个观察值和 35 个变量组成的数据框。

一个变量的正态性检验将是

shapiro.test(mydata$var1)

我想同时测试所有变量的正态性。我怎样才能在 R 中做到这一点？

Answer 1

这在很大程度上取决于您正在寻找什么样的输出以及您希望做什么来控制家庭错误率。这是使用 Bonferroni 校正的一种解决方案

# example data
t <- as.data.frame(matrix(rnorm(700), 20,35))

shapiro_test_df <- function(df, bonf= TRUE, alpha= 0.05) {
  l <- lapply(df, shapiro.test)
  s <- do.call("c", lapply(l, "[[", 1))
  p <- do.call("c", lapply(l, "[[", 2))
  if (bonf == TRUE) {
    sig <- ifelse(p > alpha / length(l), "H0", "Ha")
  } else {
    sig <- ifelse(p > alpha, "H0", "Ha")
  }
  return(list(statistic= s,
              p.value= p,
              significance= sig,
              method= ifelse(bonf == TRUE, "Shapiro-Wilks test with Bonferroni Correction",
                             "Shapiro-Wilks test without Bonferroni Correction")))
}

shapiro_test_df(t)