studio - text in r
En R, ¿cómo encontrar el error estándar de la media? (7)
¿Hay algún comando para encontrar el error estándar de la media en R?
El error estándar es solo la desviación estándar dividida por la raíz cuadrada del tamaño de muestra. Para que pueda hacer fácilmente su propia función:
> std <- function(x) sd(x)/sqrt(length(x))
> std(c(1,2,3,4))
[1] 0.6454972
El paquete sciplot tiene la función incorporada se (x)
Probablemente sea más eficiente usar var ... ya que en realidad se inserta dos veces en el código, una para obtener el sd (el código para sd está en r y se revela simplemente escribiendo "sd") ...
se <- function(x) sqrt(var(x)/length(x))
Una versión de la respuesta de John que elimina las fastidiosas NA:
stderr <- function(x) sqrt(var(x,na.rm=TRUE)/length(na.omit(x)))
más en general, para errores estándar en cualquier otro parámetro, puede usar el paquete de arranque para simulaciones de arranque (o escríbalos por su cuenta)
y <- mean(x, na.rm=TRUE)
sd(y) para la desviación estándar var(y) para la varianza.
Ambas derivaciones usan n-1 en el denominador, por lo que se basan en datos de muestra.