[R] Impaired boxplot functionality - mean instead of median
maustin at amgen.com
Fri Dec 2 01:32:50 CET 2005
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,
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
> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch
> [mailto:r-help-bounces at 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 at 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
> > 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?
> R-help at stat.math.ethz.ch mailing list
> PLEASE do read the posting guide!
More information about the R-help