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]] ______________________________________________ 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.