# Load the 'sicalis' dataset
data(sicalis)
# Generate a visual model
vm1 <- vismodel(sicalis)
# Make a tetracolorspace
tcs1 <- tcs(vm1)
# Subset all 'crown' patches (C in file names)
head(subset(sicalis, "B"))
subset(vm1, "B")
subset(tcs1, "B")[, 1:5]
Run the code above in your browser using DataLab