Try require(Hmisc) ?panel.bpplot This implements extended box plots that can show a variety of quantiles. Frank
Paul Raftery wrote: > > Hi all, > > I'm just getting started with R and I would appreciate some help. I'm > having > trouble creating a boxplot with whiskers at the 95th and 5th percentiles > instead of at 1.5 * IQR. I have read the relevant documentation, and > checked > existing mails on this topic. I found a small modification that should > work > : https://stat.ethz.ch/pipermail/r-help/2001-November/016817.html and > tried > to implement it. > > Basically, it says to replace boxplot.stats with: > > myboxplot.stats <- function (x, coef = NULL, do.conf = TRUE, do.out = > TRUE) > { > nna <- !is.na(x) > n <- sum(nna) > stats <- quantile(x, c(.05,.25,.5,.75,.95), na.rm = TRUE) > iqr <- diff(stats[c(2, 4)]) > out <- x < stats[1] | x > stats[5] > conf <- if (do.conf) > stats[3] + c(-1.58, 1.58) * diff(stats[c(2, 4)])/sqrt(n) > list(stats = stats, n = n, conf = conf, out = x[out & nna]) > } > > I entered the new function, and used fix(boxplot.default) to modify > boxplot.default so that it references myboxplot.stats instead of the > original boxplot.stats function. > > If I now type boxplot.default, I can see that the code has been modified > as > expected. However, I get the exact same result as before when I create a > boxplot - it shows the whiskers at 1.5 * IQR. You can test this out by > creating a boxplot from the iris dataset supplied with R using > boxplot(iris$Sepal.Length ~ iris$Species). You see that the boxplot is the > same before and after the fix. Does anybody know why this occurs, and how > I > can get around this issue? > > Thanks, > -- > Regards, > Paul > > > ===================== > Contact Details > ===================== > Paul Raftery, BEng(Hons) (Mech), Fulbright Fellow, PhD > http://www.paulraftery.com/ > <http://www.paulraftery.com/> > Postdoctoral Research Engineer > Informatics Research Unit for Sustainable Engineering (IRUSE) > http://www.iruse.ie/ > > Department of Civil Engineering, > National University of Ireland, Galway, > University Road, > Galway, > Ireland. > > Landline: +353 91 49 3086 > Mobile: +353 85 124 7947 > Skype: praftery > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@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. > ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Box-plot-with-5th-and-95th-percentiles-instead-of-1-5-IQR-problems-implementing-an-existing-solution-tp3456123p3456187.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@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.
