This function plots a boxplot of the data and the boundaries are based on percentiles.
boxplotperc(x, ..., quant = c(0.02, 0.98), width = NULL, varwidth = FALSE,
notch = FALSE, outline = TRUE, names, plot = TRUE, border = par("fg"),
col = NULL, log = "", pars = list(boxwex = 0.8, staplewex = 0.5, outwex = 0.5),
horizontal = FALSE, add = FALSE, at = NULL)
a vector of length 5, containing the extreme of the lower whisker, the lower "hinge", the median, the upper "hinge" and the extreme of the upper whisker (backtransformed)
the number of non-NA observations in the sample
the lower and upper extremes of the "notch"
the values of any data points which lie beyond the extremes of the whiskers (backtransformed)
the group
the attributes
data
further arguments for creating the list
the underlying percentages
a vector giving the relative widths of the boxes making up the plot
if varwidth is TRUE, the boxes are drawn with widths proportional to the square-roots of the number of observations in the groups.
if notch is TRUE, a notch is drawn in each side of the boxes
if outliers is FALSE, the outliers are not drawn
define the names of the attributes
if plot is TRUE the boxplot is plotted in the current plot
character or numeric (vector) which indicates the color of the box borders
defines the colour
character, indicating if any axis should be drawn in logarithmic scale
some graphical parameters can be specified
logical parameter indicating if the boxplots should be horizontal; FALSE means vertical boxes
if TRUE the boxplot is added to the current plot
numeric vector giving the locations of the boxplots
Peter Filzmoser <P.Filzmoser@tuwien.ac.at> http://cstat.tuwien.ac.at/filz/
The default value for quant is the 2% and 98% quantile and this argument defines the percentiles for the upper and lower whiskers.
C. Reimann, P. Filzmoser, R.G. Garrett, and R. Dutter: Statistical Data Analysis Explained. Applied Environmental Statistics with R. John Wiley and Sons, Chichester, 2008.
boxplotlog
data(chorizon)
Ba=chorizon[,"Ba"]
boxplotperc(Ba,quant=c(0.05,0.95),horizontal=TRUE,xlab="Ba [mg/kg]",cex.lab=1.2,pch=3)
Run the code above in your browser using DataLab