Re: [R] Custom Plot - means, SD & 5th-95th% (Plotmeans or Boxplot)

From: <h.wickham_at_gmail.com>
Date: Sun, 17 Feb 2008 21:03:46 -0600

It's fairly simple to set up something like this for ggplot2:

install.packages("ggplot2")
library(ggplot2)

library(ggplot2)

q5 <- function(data) {
 q <- function(p) unname(quantile(data$y, p))  data.frame(min = q(0.05), max = q(0.95)) }

ggplot(diamonds, aes(x = cut, y = price)) + stat_summary(fun="mean_sdl", geom="crossbar", mult=1) + stat_summary(fun="q5", geom="linerange")

(see Hmisc::smean.sdl for more details on the first summary function).

Although from this example you can see that that choice of summary statistics probably isn't the best. You can read more about stat_summary (and ggplot2 in general) at http://had.co.nz/ggplot2/stat_summary.html.

Hadley

On 2/17/08, Stropharia <stevenworthington_at_hotmail.com> wrote:
>
> Any help with this problem would be greatly appreciated:
>
> I need to produce a custom plot i haven't come across in R. Basically, I
> want to show means, 1st standard deviation and 5th and 95th percentiles
> visually, using something resembling a boxplot. Is it possible to completely
> customize a boxplot so that it shows means as the bar (instead of, not as
> well as medians), standard deviations at the hinges (instead of IQR) and 5th
> & 95th percentiles at the brackets? The plotmeans function (ggplots) allows
> means to be plotted, but it seems only with confidence intervals, not 5th
> and 95th percentiles, and also without hinges (for standard deviations).
> I've searched the forums and various books and have drawn a blank on this.
> Thanks.
>
> Steve
>
> ~~~~~~~~~~~~~~~~
> Steve Worthington
> Ph.D. Candidate
> New York Consortium
> in Evolutionary Primatology
> --
> View this message in context:
> http://www.nabble.com/Custom-Plot---means%2C-SD---5th-95th--%28Plotmeans-or-Boxplot%29--tp15537706p15537706.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help_at_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.
>

-- 
http://had.co.nz/

______________________________________________
R-help_at_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.
Received on Mon 18 Feb 2008 - 03:06:45 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 18 Feb 2008 - 04:30:15 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive