## Not run:
# # For plotting
# data(sicalis)
# vis.sicalis <- vismodel(sicalis, visual = 'avg.uv')
# tcs.sicalis <- colspace(vis.sicalis, space = 'tcs')
# plot(tcs.sicalis, size = 0.005)
# rgl.postscript('testplot.pdf',fmt='pdf')
# rgl.snapshot('testplot.png')
#
# # For adding points
# patch <- rep(c('C', 'T', 'B'), 7)
# tcs.crown <- subset(tcs.sicalis, 'C')
# tcs.breast <- subset(tcs.sicalis, 'B')
# plot(tcs.crown, col ='blue')
# points(tcs.breast, col ='red')
#
# # For plotting convex hull
# plot(tcs.sicalis, col = 'blue', size = 0.005)
# vol(tcs.sicalis)
# ## End(Not run)
Run the code above in your browser using DataLab