В R, как найти стандартную ошибку среднего значения?


есть ли команда, чтобы найти стандартную ошибку среднего значения в R?

7 64

7 ответов:

здесь plotrix пакет имеет встроенную функцию для этого: std.ошибка

Стандартная ошибка-это просто стандартное отклонение, деленное на квадратный корень из размера выборки. Так что вы можете легко сделать свою собственную функцию:

> std <- function(x) sd(x)/sqrt(length(x))
> std(c(1,2,3,4))
[1] 0.6454972

вероятно, более эффективно использовать var... поскольку вы на самом деле sqrt дважды в своем коде, один раз, чтобы получить sd (код для sd находится в r и раскрывается, просто набрав "sd")...

se <- function(x) sqrt(var(x)/length(x))

версия ответа Джона выше, которая удаляет надоедливые NA:

stderr <- function(x) sqrt(var(x,na.rm=TRUE)/length(na.omit(x)))

пакет sciplot имеет встроенную функцию se (x)

В общем случае, для стандартных ошибок по любому другому параметру вы можете использовать загрузочный пакет для моделирования bootstrap (или написать их самостоятельно)

y <- mean(x, na.rm=TRUE)

sd(y) для стандартного отклонения var(y) для дисперсии.

оба вывода используют n-1 в знаменателе, поэтому они основаны на образцах данных.