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