Я пытался запустить симуляцию, чтобы определить, каково стандартное отклонение суммы 16 стандартных нормальных переменных

Question

2 * (1 - pnorm(z_score, mean = 0, sd = 1))
2 * (1 - pnorm(diff, mean = 0, sd = se))

Я попробовал вышеуказанные команды, и я новичок в R, пожалуйста, помогите.

Answer 1

Хорошо, если я понимаю, что вы хотите сделать хорошо, следующий код будет работать для вас, при условии, что стандартные нормали независимы.

X <- matrix(NA, nrow = 10000, ncol = 16)
for(i in 1:16) X[,i] <- rnorm(10000)
Y <- rowSums(X)

rbind(mean = mean(Y), var = var(Y), sd = sd(Y))
            [,1]
mean  0.03440269
var  16.17512070
sd    4.02183052

Стандартное отклонение будет равно 4 для суммы нормальных (независимых) переменных, поскольку дисперсия является суммой дисперсий.И, наконец, стандартное отклонение - это квадратный корень из дисперсии.