Use o código a seguir para produzir vários gráficos de caixa, classificados pelo valor médio das variáveis:

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)

(Veja este post para mais detalhe)

Quando altero o código parahorizontal = TRUE, Não consigo alinhar os pontos com os gráficos de caixa. Alguma idéia de como adicionar corretamentepoints para boxplots horizontais?