Есть ли какая-нибудь команда, чтобы найти стандартную ошибку среднего в R?
Стандартная ошибка - это просто стандартное отклонение, деленное на квадратный корень размера выборки. Таким образом, вы можете легко сделать свою собственную функцию:
> 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, na.rm=FALSE) { if (na.rm) x <- na.omit(x) sqrt(var(x)/length(x)) }
Существует пакет plotrix со встроенной функцией для этого: std.error
Пакет Sciplot имеет встроенную функцию se (x)
y <- mean(x, na.rm=TRUE)
sd(y) для стандартного отклонения var(y) для отклонения.
sd(y)
var(y)
Оба деривата используют n-1 в знаменателе, поэтому они основаны на данных выборки.
n-1
в более общем смысле, для стандартных ошибок по любому другому параметру, вы можете использовать загрузочный пакет для моделирования начальной загрузки (или написать их самостоятельно)