[R] Impaired boxplot functionality - mean instead of median

Austin, Matt
2005-12-02T00:32:50Z

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
>