## Not run: ------------------------------------
# data(sicalis)
# bysic <- gsub("^ind[0-9].",'', names(sicalis)[-1])
# aggplot(sicalis, bysic)
# aggplot(sicalis, bysic, shade = spec2rgb(sicalis), lcol = 1)
# aggplot(sicalis, bysic, lcol = 1, FUN.error = function(x) sd(x)/sqrt(length(x)))
## ---------------------------------------------
Run the code above in your browser using DataLab