performance de rho2hat et ppm dans le spatstat

Question

Je veux comparer la performance de l'estimateur non paramétrique de la fonction d'intensité "rho2hat" avec "ppm".

Ma question comment puis-je effectuer un test pour voir l'estimation fonctionne mieux? Je ne pouvais pas utiliser quadrat.test?

Answer 1

Comme mentionné dans le commentaire par @ adrian-Baddeley votre stratégie suggérée va simplement mesurer la différence entre les deux résultats. La façon dont vous pouvez calculer le nombre attendu de points dans différentes régions d'une intensité estimée est comme suit:

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