tests/15-Kcross.R

# Test the edge correction in K
library(devtools)
load_all(".")
library(spatstat)

x0 <- x <- rmpoispp(rep(100, 20)*2)
r <- seq(0, .3, l=20)

#Rprof()
t0 <- system.time(  k <- K_cross_all_box(x, r)  )
t1 <- system.time(  k1 <- K_cross_all_box(x, r, correction=F)  )
#summaryRprof()

plot(k[1,1,])
points(k1[1,1,])
antiphon/Kcross documentation built on Oct. 21, 2024, 5:45 p.m.