Hi:

Following up on Ben's suggestion re ggplot2, here's a manufactured example:

# Fake data:
mortrate <- round(runif(100), 3)
dd <- data.frame(rate = mortrate, moe = 1.96 * sqrt(mortrate * (1 -
mortrate))/10,
                  hosp = factor(paste('H', 1:100, sep = '')))
dim(dd)
[1] 100   3

# Set up the plot: ymin and ymax are input parameters to geom_pointrange
g <- ggplot(dd, aes(x = hosp, y = rate, ymin = rate - moe, ymax = rate +
moe))

# Version with hospitals in lexicographic order (rather useless):
g + geom_pointrange() +
    geom_hline(aes(yintercept = mean(rate))) +
    scale_x_discrete(expand = c(0.01, 0.01)) +
    coord_flip() + theme_bw() +
    opts(axis.text.y = theme_text(size = 5, hjust = 1))

# Version with hospitals sorted by rate (better):
g + geom_pointrange(aes(x = reorder(hosp, rate, mean))) +
    geom_hline(aes(yintercept = mean(rate))) +
    scale_x_discrete(expand = c(0.01, 0.01)) +
    coord_flip() + theme_bw() +
    opts(axis.text.y = theme_text(size = 5, hjust = 1))

Hope you like it...
Dennis

On Tue, Oct 19, 2010 at 4:59 PM, XINLI LI <lihawaii@gmail.com> wrote:

> Dear R Users:
>
>       I have the individual mortality rate and 95% CI of 100 hospitals,
> how  to do the  plot with the individual hospital in the Yaxis, and the
> mortality rate and 95% CI in the Xais and a overall mean as a reference
> line?
>
> Thanks and regards,
>
> Xin
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]

