Build multiple data series in a single package with R

Question

Let's look at an example for a tablet , which is included in the aplpack library in R. A baggage bag is a two-dimensional generalization of a box and, therefore, gives an idea of ​​the distribution of data points in both axes.

Package example:

Code for example:

# example of Rousseeuw et al., see R-package rpart cardata <- structure(as.integer( c(2560,2345,1845,2260,2440, 2285, 2275, 2350, 2295, 1900, 2390, 2075, 2330, 3320, 2885, 3310, 2695, 2170, 2710, 2775, 2840, 2485, 2670, 2640, 2655, 3065, 2750, 2920, 2780, 2745, 3110, 2920, 2645, 2575, 2935, 2920, 2985, 3265, 2880, 2975, 3450, 3145, 3190, 3610, 2885, 3480, 3200, 2765, 3220, 3480, 3325, 3855, 3850, 3195, 3735, 3665, 3735, 3415, 3185, 3690, 97, 114, 81, 91, 113, 97, 97, 98, 109, 73, 97, 89, 109, 305, 153, 302, 133, 97, 125, 146, 107, 109, 121, 151, 133, 181, 141, 132, 133, 122, 181, 146, 151, 116, 135, 122, 141, 163, 151, 153, 202, 180, 182, 232, 143, 180, 180, 151, 189, 180, 231, 305, 302, 151, 202, 182, 181, 143, 146, 146)), .Dim = as.integer(c(60, 2)), .Dimnames = list(NULL, c("Weight", "Disp."))) bagplot(cardata,factor=3,show.baghull=TRUE, show.loophull=TRUE,precision=1,dkmethod=2) title("car data Chambers/Hastie 1992") # points of y=x*x bagplot(x=1:30,y=(1:30)^2,verbose=FALSE,dkmethod=2) 

The aplpack bag only seems to support the creation of a “bag” for one data series. Even more interesting would be to build two (or three) series of data in one package, where a visual comparison of the "bags" of a data series gives an idea of ​​the differences in the distribution of data in a data series. Does anyone know if (and if so, how) it can be done in R?

Answer 1

If we modify some aplpack::bagplot , we can create a new geometry for ggplot2. Then we can compare the groups in the dataset using the usual ggplot2 methods. Here is one example:

 library(ggplot2) ggplot(iris, aes(Sepal.Length, Sepal.Width, colour = Species, fill = Species)) + geom_bag() + theme_minimal() 

and we can show the points with the bag:

 ggplot(iris, aes(Sepal.Length, Sepal.Width, colour = Species, fill = Species)) + geom_bag() + geom_point() + theme_minimal() 

Here is the code for the geom_bag function and the modified aplpack::bagplot : https://gist.github.com/benmarwick/00772ccea2dd0b0f1745