## 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)))
# ## End(Not run)
Run the code above in your browser using DataLab