Utilizo el siguiente código para generar múltiples diagramas de caja, clasificados por el valor medio de las variables:

zx <- replicate (5, rnorm(50))
zx_means <- (colMeans(zx, na.rm = TRUE))
colnames (zx) <- seq_len (ncol (zx))
boxplot(zx [, order (zx_means)], horizontal = FALSE, outline = FALSE)
points(zx_means [ order (zx_means)], pch = 22, col = "darkgrey", lwd = 7)

(Ver esta publicación para más detalles)

Cuando cambio el código ahorizontal = TRUE, No puedo hacer que los puntos se alineen con los diagramas de caja. Alguna idea sobre cómo agregar correctamentepoints a diagramas de caja horizontales?