Re: [R] Impaired boxplot functionality - mean instead of median

From: Austin, Matt <maustin_at_amgen.com>
Date: Fri 02 Dec 2005 - 11:32:50 EST

Check your syntax on the bwplot call.

fa <- data.frame(doz=sample(500:2000, size=500), fabp2=rep(1:20, 25))

bwplot(factor(fabp2) ~ doz, data=fa, panel=panel.bpplot)

fa.sum <- summarize( fa$doz, list( fabp2 = fa$fabp2), smean.sd, stat.name="doz")

Dotplot( factor(fabp2) ~ Cbind(doz, doz - SD, doz + SD) , data=fa.sum)

You can ignore the warning, I'm sure Dr. Harrell has already fixed that issue.

--Matt

> -----Original Message-----
> From: r-help-bounces@stat.math.ethz.ch
> [mailto:r-help-bounces@stat.math.ethz.ch]On Behalf Of Evgeniy Kachalin
> Sent: Thursday, December 01, 2005 2:43 PM
> To: Frank E Harrell Jr
> Cc: r-help@stat.math.ethz.ch
> Subject: Re: [R] Impaired boxplot functionality - mean
> instead of median
>
>
> Frank E Harrell Jr пишет:
> > Evgeniy Kachalin wrote:
> >
> >> Marc Schwartz (via MN) пишет:
> >>
> >>>> Marc Schwartz (via MN) пишет:
>
> >>
> >
> > library(Hmisc)
> > library(lattice)
> > ?panel.bpplot
> >
> > bwplot(...., panel=panel.bpplot)
> >
> > By default, panel.bpplot shows the mean (dot) and median
> (line) plus
> > several quantiles. To bother Martin in a friendly way, I
> think that
> > means can be useful additions - not that they are so useful by
> > themselves, but that when they differ a lot from the median,
> > non-statisticians gain further information about asymmetry.
> Also, even
> > though the simple box plot is elegant, I sometimes think it
> has a high
> > ink to information ratio. I have gained a lot from seeing outer
> > quantiles on the plot, and I don't like to show outer
> points for fear of
> > someone labeling them outliers. For describing raw data
> distributions,
> > I never find standard deviations useful, however.
> >
>
> => fa
> doz fabp2
> 1 900 2
> 4 1500 2
> 6 1000 2
> 8 750 3
> 10 750 1
> 11 1750 2
> 12 500 3
> ....
> ....
> ....
> ....
>
>
> > bwplot(doz~factor(fabp2),data=fa,panel=panel.bpplot)
> Error in sort(x, partial = unique(c(lo, hi))) :
> unsupported options for partial sorting
>
>
> That's NOT simple way.
>
> I need just one change.
> Is there any good way?
> $-(
>
> --
> Evgeniy
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Dec 02 11:39:00 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:41:26 EST