web-dev-qa-db-fra.com

Dans R, comment trouver l'erreur type de la moyenne?

Existe-t-il une commande permettant de trouver l'erreur type de la moyenne dans R?

rstatistics
74
alex

Il y a le paquet plotrix avec une fonction intégrée pour ceci: std.error

29
Matt Ball

L'erreur type est simplement l'écart type divisé par la racine carrée de la taille de l'échantillon. Ainsi, vous pouvez facilement créer votre propre fonction:

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

Il est probablement plus efficace d'utiliser var ... puisque vous écrivez deux fois dans votre code, une fois pour obtenir le sd (le code pour sd est en r et est révélé en tapant simplement "sd") ...

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

Une version de la réponse de John ci-dessus qui supprime les satanés NA:

stderr <- function(x, na.rm=FALSE) {
  if (na.rm) x <- na.omit(x)
  sqrt(var(x)/length(x))
}
55
arandomlypickedname

Le paquet sciplot a la fonction intégrée se (x)

4
user2045554

plus généralement, pour les erreurs standard sur tout autre paramètre, vous pouvez utiliser le package de démarrage pour des simulations bootstrap (ou les écrire vous-même))

0
owen