## 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