if (FALSE) {
# example data, see manual page for details
data(sp5)
# slice top 25 cm
s <- slice(sp5, 1:25 ~ .)
# check some of the data
par(mar=c(0,0,0,0))
plot(sample(s, 25), divide.hz=FALSE, name='', print.id=FALSE, width=0.5)
# colors
previewColors(unique(s$soil_color))
# compute marginal quantiles and L1 median
cq <- colorQuantiles(s$soil_color)
# simple graphical display of results
plotColorQuantiles(cq)
}
Run the code above in your browser using DataLab