Hello,

I know this has been covered on here before, but as a complete novice, I
need a little more guidance.  I would like to produce boxplots with the
whiskers extending to the 10 and 90th percentiles.  I found this code:

myboxplot.stats <-  function (x, coef = NULL, do.conf = TRUE, do.out =
TRUE)
{
  nna <- !is.na(x)
  n <- sum(nna)
  stats <- quantile(x, c(.1,.25,.5,.75,.9), 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])
}

posted by Mr. Jim Bowers, and additional posts discussing how to make
it work.  The issue I am having is that all the posts say to edit
boxplot.default, but I have no idea how to actually do that.  I've
tried fix(bowplot.default)

and fixInNamespace(...), but what do I actually change?  I tried
including an argument "stats=myboxplot.stats" but that did not change
anything.  Once it is changed, can I just use the same "boxplot(...)"
code as normal or

do I need to use myboxplot?  I know this should be simple, but it is
generating a lot of frustration.  Any help would be greatly
appreciated.

Thanks, James

	[[alternative HTML version deleted]]

