if (FALSE) {
# example data, see manual page for details
data(sp5)
# slice top 25 cm
# 24-25cm is the last slice
s <- dice(sp5, 0:24 ~ .)
# check some of the data
par(mar=c(0,0,0,0))
plotSPC(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