Función para error estándar de la media
plotrix::std.error()
Cuantas veces he calculado el error estándar de la media a “mano”, sobre todo, cuando quiero dibujar las barras de error en gráficos de barras. Pues bien, el paquete plotrix contiene una función directa para calcular el error estándar.
plotrix::std.error()
Puede resultar una tontería, pero añadir funciones a tu rutina diaría aunque solo sean para cubrir tareas tan faciles como calcular el error estándar, me ahorra tiempo y posibles errores.
Ejemplo
data_plot =
dplyr::starwars %>%
group_by(gender) %>%
summarise(media_altura = mean(height,na.rm=T),
se_altura = plotrix::std.error(height,na.rm = T)) %>%
drop_na()
Ejemplo con gráfico
Lo dicho, uso especialmente la función cuando quiero añadir las barras de error en gráficos de barras. Dejo un ejemplo rápido abajo. Aprovecho para introducir fugazmente el paquete bbplot para aplicar un tema a los objetos de ggplot muy parecido al estilo que sigue la BBC
ggpubr::ggbarplot(data = data_plot,
x= "gender",
y= "media_altura",
fill = "gender")+
geom_errorbar(data =data_plot ,
aes(ymin= media_altura-se_altura, ymax=media_altura+se_altura),
size = 0.5,
width = 0.5)+
labs(title = "Altura media por género en el universo StarWars",
subtitle = "Un ejemplo rápido para visualizar las barras de error estándar de la media",
caption = "Hecho por Adrián Muñoz el 10Ago del 22")+
ggplot2::coord_cartesian(ylim = c(100,180))+
scale_fill_grey(start = 1, end = 0.5,)+
bbplot::bbc_style()+
theme(legend.position = "none")
Adrián Muñoz García
Data Scientist and Statistics Professor
My research interests are Behavioral Economics, Cognition and Methods