How to draw a chart with sorted horizontal error bars (sorted barcharts with error marks)?

lattice, plot, r

Solution

See R-Help: Adding error bars to lattice plots

prepanel.ci <- function(x, y, lx, ux, subscripts, ...) {
    x <- as.numeric(x)
    lx <- as.numeric(lx[subscripts])
    ux <- as.numeric(ux[subscripts])
    list(xlim = range(0, x, ux, lx, finite = TRUE))
}


panel.ci <- function(x, y, lx, ux, subscripts, ...) {
    x <- as.numeric(x)
    y <- as.numeric(y)
    lx <- as.numeric(lx[subscripts])
    ux <- as.numeric(ux[subscripts])
    panel.barchart(x, y, ...)
    panel.arrows(lx, y, ux, y, col = 'black',
                 length = 0.25, unit = "native",
                 angle = 90, code = 3)
}

p <- barchart(reorder(var, mean) ~ mean, data=plot.data,
              lx=plot.data$mean-plot.data$error,
              ux=plot.data$mean+plot.data$error,
              panel=panel.ci,
              prepanel=prepanel.ci)
print(p)

Problem

I would like to plot means and standard errors as a horizontal barchart, and I want the mean sorted. I've found the way to plot horizontal sorted barcharts using lattice, but I do not know how to add error marks. The following are my data and the R code I came up with. ``` data <- structure(c(0.67, 0.67, 0.76, 0.66, 0.71, 0.6, 0.52, 0.6, 0.71, 0.76, 0.76, 0.71, 0.6, 0.61, 0.9, 0.5, 0.58, 0.84, 0.68, 0.88, 0.89, 0.96, 1, 0.95, 1, 1, 0.98, 0.78, 0.98, 1, 1, 0.99, 1, 1, 0.95, 0.92, 1, 0.91, 1, 0.87, 0.91, 0.72, 0.73, 0.55, 0.82, 0.87, 0.64, 0.75, 0.75, 1, 0.81, 0.79, 1, 0.74, 0.57, 0.84, 1, 0.95, 0.78, 0.95), .Dim = c(20L, 3L), .Dimnames = list( c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20"), c("A", "B", "C"))) means <- apply(data, 2, mean) errors <- apply(data, 2, sd) plot.data <- data.frame(colnames(data), means, errors) colnames(plot.data) <- c("var", "mean", "error") library("lattice") plot.new() barchart(reorder(var, mean) ~ mean, plot.data, xlim = c(0, 1)) ``` Is there any way to add error marks to this chart? If not, any suggestion on how to plot the chart I want in R? Thank you in advance!

Original source