## Not run:
# data(sicalis)
# vis.sicalis <- vismodel(sicalis, visual='avg.uv')
# tcs.sicalis <- tcs(vis.sicalis)
# projplot(tcs.sicalis, pch=16, col=setNames(rep(1:3, 7), rep(c('C','T','B'),7))) ## End(Not run)
Run the code above in your browser using DataLab