[R] Чудеса науки, техники и теста Шапиро-Уилка

test <- function() {
  r <- c();
  for (i in 1:100) {
   r <- append(r, shapiro.test(rnorm(5000))$statistic)
  }
  summary(r);
}

может p.value криво аппроксимируется. потому что значение статистики нормальное вроде

p.s. всё забыл уже :(

p.value: an approximate p-value for the test. This is said in Royston
(1995) to be adequate for 'p.value < 0.1'.

help(shapiro.test)

Оцените улучшенную реализацию этой функции:

test2 <- function() {
  summary(sapply(1:100, function(x) {shapiro.test(rnorm(5000))$p.value}))
}

Второй вариант лучше первого?

Мне кажется первый понятнее

> Мне кажется первый понятнее

Вторая версия типа в функциональном стиле написана... типа haskell и все такое...