library(spatstat)
# Longleaf Pine data: real-valued marks represent tree diameter
data(longleaf)
# Subset of this large pattern
swcorner <- owin(c(0,100),c(0,100))
sub <- longleaf[ , swcorner]
# mark correlation function
mc <- markcorr(sub)
plot(mc$r, mc$trans,
xlab="r", ylab="mcf(r)",
type="l", ylim=c(0,2), xlim=c(0,40))
abline(1,0,lty=2)
# or:
conspire(mc, cbind(trans,theo) ~ r, xlim=c(0,40),ylim=c(0,2))
# simulated data with independent marks
X <- rpoispp(100)
X$marks <- runif(X$n)
Xc <- markcorr(X)
conspire(Xc, cbind(iso,trans,theo) ~ r, xlim=c(0,0.4),ylim=c(0,2))
# multitype pattern
data(ganglia)
gmc <- markcorr(ganglia, function(m1,m2){m1 == m2})
conspire(gmc, cbind(iso,trans,theo)~ r, ylim=c(0,2), xlim=c(0,0.4))
# shows strong negative association between types at close distances
<testonly>niets <- markcorr(ganglia, function(m1,m2){m1 == m2}, method="loess")
niets <- markcorr(X, correction="isotropic", method="smrep")</testonly>
Run the code above in your browser using DataCamp Workspace