# Lansing woods data - multitype points with 6 types
 data(lansing)
 <testonly># smaller dataset
   lansing <- lansing[ seq(1,lansing$n,by=45), ]</testonly>
 # compute 6 x 6 array of all cross-type K functions
 a <- alltypes(lansing, "K")
 # extract first three marks only
 b <- a[1:3,1:3]
 plot(b)
 # subset pertaining to hickories
 h <- a[levels(lansing$marks) == "hickory", ]
 plot(h)Run the code above in your browser using DataLab