Diagrama de caja con intervalos de confianza para la mediana

Anteriormente vimos cómo crear diagramas de caja en R. En esta ocasión vamos a crear una variante, incluyendo cuñas o muescas(notches) que muestran el intervalo de confianza alrededor de la mediana.

boxplot(iris$Sepal.Width,
        notch = TRUE,
        col = 'palegreen',
        xlab = "Sepal Width",
        horizontal = TRUE)

Las cuñas alrededor de la mediana R las calcula como +/-1,58 RIC/raíz cuadrada(n), donde RIC es el rango intercuartílico y n el número de datos.

Valores

Para obtener los valores del intervalo de confianza de las cuñas:

boxplot.stats(iris$Sepal.Width)$conf

[1] 2.935497 3.064503

Etiquetas de datos

# Horizontal
data <- iris$Sepal.Width
boxplot(data,
        notch = TRUE,
        col = 'palegreen', 
        horizontal = TRUE, 
        axes = FALSE, 
        staplewex =1)
values <- c(round(boxplot.stats(data)$conf, 1), boxplot.stats(data)$stats)
text(x = values, labels = values, y = 1.25)

Si no queremos etiquetar el intervalo de confianza de las cuñas.

values <- c(round(boxplot.stats(data)$stats, 1), boxplot.stats(data)$stats)

# Vertical
data <- iris$Sepal.Width
boxplot(data,
        notch = TRUE,
        col = 'palegreen', 
        axes = FALSE, 
        staplewex =1)
values <- c(round(boxplot.stats(data)$conf, 1), boxplot.stats(data)$stats)
text(y = values, labels = values, x = 1.25)

Entradas relacionadas

• Introducción al diagrama de caja (box plot) en R
• Diagramas de caja con ggplot2 en R. Una introducción.

Referencias

• Función boxplot
• Box Plot Statistics
• Stackoverflow