fboxplot {rainbow} | R Documentation |
Compute bivariate bagplot, functional bagplot and bivariate HDR boxplot, functional HDR boxplot.
fboxplot(data, plot.type = c("functional", "bivariate"),type = c("bag", "hdr"), alpha = c(0.01, 0.5), h, factor = 2.57, na.rm = TRUE, xlab = data$xname, ylab = data$yname, ...)
data |
An object of class fd or fts . |
plot.type |
Version of boxplot. When plot.type="functional" , a functional plot is provided. When plot.type="bivariate" , a square bivariate plot is provided. |
type |
Type of boxplot. When type="bag" , a bagplot is provided. When type="hdr" , a HDR boxplot is provided. |
alpha |
Coverage probability for the functional HDR boxplot. alpha are the coverage percentages of the outliers and the central region. |
h |
When type="hdr" , bandwidth is selected by using the smoothed cross validation. |
factor |
When type="bag" , the outer region of a bagplot is the convex hull obtained by inflating the inner region by the bagplot factor. |
na.rm |
Remove missing values. |
xlab |
A title for the x axis. |
ylab |
A title for the y axis. |
... |
Other arguments. |
The functional curves are first projected into a finite dimensional subspace. For simiplicity, we choose the subspace as R^2. Based on Tukey (1974)'s halfspace bagplot and Hyndman (1996)'s HDR boxplot, we order each data point in R^2 by data depth and data density. Outliers are those that have either lowest depth or lowest density.
NULL.
Rob J Hyndman, Han Lin Shang
J. W. Tukey (1974) "Mathematics and the picturing of data", Proceedings of the International Congress of Mathematicians, 2, 523-532, Canadian Mathematical Congress, Montreal.
R. Hyndman (1996) "Computing and graphing highest density regions", The American Statistician, 50, 120-126.
R. J. Hyndman and H. L. Shang. (2008) "Rainbow plots, bagplots, and boxplots for functional data", Tech. rep. Department of Econometrics & Business Statistics, Monash University, November 2008, http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/2008/9-08.php.
fboxplot(data = ElNino, plot.type = "functional", type = "bag") fboxplot(data = ElNino, plot.type = "bivariate", type = "bag") fboxplot(data = ElNino, plot.type = "functional", type = "hdr", alpha = c(0.07,0.5)) fboxplot(data = ElNino, plot.type = "bivariate", type = "hdr", alpha = c(0.07,0.5))