производительность rho2hat и ppm в spatstat

Как упоминалось в комментарии @ adrian-baddeley, предлагаемая вами стратегия просто измеряет разницу между двумя результатами. Вы можете рассчитать ожидаемое количество точек в разных регионах на основе оценочной интенсивности следующим образом:

library(spatstat)
#> Loading required package: nlme
#> Loading required package: rpart
#> 
#> spatstat 1.51-0.035       (nickname: 'Spatfefe') 
#> For an introduction to spatstat, type 'beginner'
elev <- bei.extra$elev
grad <- bei.extra$grad
non <- rho2hat(bei, elev, grad)
pred <- predict(non)
grid <- quadrats(pred, nx = 4, ny = 2)
by(pred, grid, integral)
#> 1:
#> [1] 522.0247
#> 
#> 2:
#> [1] 503.6255
#> 
#> 3:
#> [1] 344.8552
#> 
#> 4:
#> [1] 522.834
#> 
#> 5:
#> [1] 454.5958
#> 
#> 6:
#> [1] 470.2484
#> 
#> 7:
#> [1] 556.1989
#> 
#> 8:
#> [1] 415.8889