[R] Impaired boxplot functionality - mean instead of median

[Austin, Matt]

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 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 
> 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 at 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
>