## Not run:
# data(sicalis)
# vis.sicalis <- vismodel(sicalis, visual = 'avg.uv')
# tcs.sicalis <- colspace(vis.sicalis, space = 'tcs')
# 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