tests/60-pcf-anin-fry.R

In antiphon/Kdirectional: Analysis of anisotropy in point patterns using second order statistics

#' anisotropic inhomogeneous pcf, raw fry point smoothing version (rubbish)

library(devtools)
load_all(".")
library(rstrauss)

if(!exists("pp")){
  set.seed(1)
  d3 <- T
  comp <- 0.6
  #bb <- cbind(c(0,comp) -comp/2, c(0,1/comp)-0.5/comp)
  side <- c(1,2)
  side2<- c(1,3)
  bbm <- cbind(side2,side)
  M <- diag(c(1/comp, comp))
  if(d3)M <- diag(c(sqrt(1/comp),sqrt(1/comp), comp))
  if(d3) bbm<-cbind(bbm, side)
  bb <- bbm%*%solve(M)
  pp0 <- rstrauss(500, .01, R<-0.03+ifelse(d3,0.08,0), perfect=TRUE, bbox=bb, iter=2e4)
  xm <- pp0$x%*%M
  pp <- list(x=xm, bbox=bbm)
}


r = seq(0,.6, l = 100)
l <- rep(nrow(pp$x)/bbox_volume(pp$bbox), nrow(pp$x))

#
pd <- pcf_anin_fry(pp, lambda=l, r = r, border=0)

# compare
p <- pcf_anin(pp, lambda=l, r=r, epsilon = e<-pi/4, stoyan=f<-0.35)
pc <- pcf_anin_cylinder(pp, lambda=l, r=r, epsilon = 0.05, stoyan=f)

par(mfrow=c(2,2))
plot(p, main="sector")
plot(pc, main="cylinder")
plot(pd, main ="fry")

# border correction need to be checked