Statistical Science

Score, Pseudo-Score and Residual Diagnostics for Spatial Point Process Models

Adrian Baddeley, Ege Rubak, and Jesper Møller

Full-text: Open access


We develop new tools for formal inference and informal model validation in the analysis of spatial point pattern data. The score test is generalized to a “pseudo-score” test derived from Besag’s pseudo-likelihood, and to a class of diagnostics based on point process residuals. The results lend theoretical support to the established practice of using functional summary statistics, such as Ripley’s K-function, when testing for complete spatial randomness; and they provide new tools such as the compensator of the K-function for testing other fitted models. The results also support localization methods such as the scan statistic and smoothed residual plots. Software for computing the diagnostics is provided.

Article information

Statist. Sci., Volume 26, Number 4 (2011), 613-646.

First available in Project Euclid: 28 February 2012

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Compensator functional summary statistics model validation point process residuals pseudo-likelihood


Baddeley, Adrian; Rubak, Ege; Møller, Jesper. Score, Pseudo-Score and Residual Diagnostics for Spatial Point Process Models. Statist. Sci. 26 (2011), no. 4, 613--646. doi:10.1214/11-STS367.

Export citation


  • [1] Alm, S. E. (1998). Approximation and simulation of the distributions of scan statistics for Poisson processes in higher dimensions. Extremes 1 111–126.
  • [2] Atkinson, A. C. (1982). Regression diagnostics, transformations and constructed variables (with discussion). J. Roy. Statist. Soc. Ser. B 44 1–36.
  • [3] Baddeley, A. (1980). A limit theorem for statistics of spatial data. Adv. in Appl. Probab. 12 447–461.
  • [4] Baddeley, A., Møller, J. and Pakes, A. G. (2008). Properties of residuals for spatial point processes. Ann. Inst. Statist. Math. 60 627–649.
  • [5] Baddeley, A. and Turner, R. (2000). Practical maximum pseudolikelihood for spatial point patterns (with discussion). Aust. N. Z. J. Stat. 42 283–322.
  • [6] Baddeley, A. and Turner, R. (2005). Spatstat: An R package for analyzing spatial point patterns. J. Statist. Software 12 1–42.
  • [7] Baddeley, A., Turner, R., Møller, J. and Hazelton, M. (2005). Residual analysis for spatial point processes (with discussion). J. R. Stat. Soc. Ser. B Stat. Methodol. 67 617–666.
  • [8] Baddeley, A. J. (1993). Stereology and survey sampling theory. Bull. Int. Statist. Inst. 50 435–449.
  • [9] Baddeley, A. J. (1999). Spatial sampling and censoring. In Stochastic Geometry (Toulouse, 1996). Monogr. Statist. Appl. Probab. 80 37–78. Chapman & Hall/CRC, Boca Raton, FL.
  • [10] Baddeley, A. J. and Møller, J. (1989). Nearest-neighbour Markov point processes and random sets. Int. Stat. Rev. 57 89–121.
  • [11] Baddeley, A. J. and van Lieshout, M. N. M. (1995). Area-interaction point processes. Ann. Inst. Statist. Math. 47 601–619.
  • [12] Barnard, G. (1963). Discussion of “The spectral analysis of point processes” by M. S. Bartlett. J. R. Stat. Soc. Ser. B Stat. Methodol. 25 294.
  • [13] Berman, M. (1986). Testing for spatial association between a point process and another stochastic process. J. Roy. Statist. Soc. Ser. C 35 54–62.
  • [14] Besag, J. (1978). Some methods of statistical analysis for spatial data. Bull. Int. Statist. Inst. 44 77–92.
  • [15] Chen, C. (1983). Score tests for regression models. J. Amer. Statist. Assoc. 78 158–161.
  • [16] Chetwynd, A. G. and Diggle, P. J. (1998). On estimating the reduced second moment measure of a stationary spatial point process. Aust. N. Z. J. Stat. 40 11–15.
  • [17] Coeurjolly, J. F. and Lavancier, F. (2010). Residuals for stationary marked Gibbs point processes. Preprint. Available at
  • [18] Conniffe, D. (2001). Score tests when a nuisance parameter is unidentified under the null hypothesis. J. Statist. Plann. Inference 97 67–83.
  • [19] Cook, R. D. and Weisberg, S. (1983). Diagnostics for heteroscedasticity in regression. Biometrika 70 1–10.
  • [20] Cordy, C. B. (1993). An extension of the Horvitz–Thompson theorem to point sampling from a continuous universe. Statist. Probab. Lett. 18 353–362.
  • [21] Cox, D. R. (1972). The statistical analysis of dependencies in point processes. In Stochastic Point Processes: Statistical Analysis, Theory, and Applications (Conf., IBM Res. Center, Yorktown Heights, N.Y., 1971) 55–66. Wiley-Interscience, New York.
  • [22] Cox, D. R. and Hinkley, D. V. (1974). Theoretical Statistics. Chapman & Hall, London.
  • [23] Cressie, N. A. C. (1991). Statistics for Spatial Data. Wiley, New York.
  • [24] Daley, D. J. and Vere-Jones, D. (2003). An Introduction to the Theory of Point Processes. Volume I: Elementary Theory and Methods, 2nd ed. Springer, New York.
  • [25] Davies, R. B. (1977). Hypothesis testing when a nuisance parameter is present only under the alternative. Biometrika 64 247–254.
  • [26] Davies, R. B. (1987). Hypothesis testing when a nuisance parameter is present only under the alternative. Biometrika 74 33–43.
  • [27] Diggle, P. J. (1979). On parameter estimation and goodness-of-fit testing for spatial point patterns. Biometrika 35 87–101.
  • [28] Diggle, P. J. (1985). A kernel method for smoothing point process data. J. Roy. Statist. Soc. Ser. C 34 138–147.
  • [29] Diggle, P. J. (2003). Statistical Analysis of Spatial Point Patterns, 2nd ed. Hodder Arnold, London.
  • [30] Georgii, H.-O. (1976). Canonical and grand canonical Gibbs states for continuum systems. Comm. Math. Phys. 48 31–51.
  • [31] Geyer, C. (1999). Likelihood inference for spatial point processes. In Stochastic Geometry (Toulouse, 1996). Monogr. Statist. Appl. Probab. 80 79–140. Chapman & Hall/CRC, Boca Raton, FL.
  • [32] Hanisch, K. H. (1984). Some remarks on estimators of the distribution function of nearest neighbour distance in stationary spatial point processes. Math. Operationsforsch. Statist. Ser. Statist. 15 409–412.
  • [33] Hansen, B. E. (1996). Inference when a nuisance parameter is not identified under the null hypothesis. Econometrica 64 413–430.
  • [34] Heinrich, L. (1988). Asymptotic behaviour of an empirical nearest-neighbour distance function for stationary Poisson cluster processes. Math. Nachr. 136 131–148.
  • [35] Heinrich, L. (1988). Asymptotic Gaussianity of some estimators for reduced factorial moment measures and product densities of stationary Poisson cluster processes. Statistics 19 87–106.
  • [36] Hope, A. C. A. (1968). A simplified Monte Carlo significance test procedure. J. R. Stat. Soc. Ser. B Stat. Methodol. 30 582–598.
  • [37] Huang, F. and Ogata, Y. (1999). Improvements of the maximum pseudo-likelihood estimators in various spatial statistical models. J. Comput. Graph. Statist. 8 510–530.
  • [38] Illian, J., Penttinen, A., Stoyan, H. and Stoyan, D. (2008). Statistical Analysis and Modelling of Spatial Point Patterns. Wiley, Chichester.
  • [39] Jensen, J. L. and Møller, J. (1991). Pseudolikelihood for exponential family models of spatial point processes. Ann. Appl. Probab. 1 445–461.
  • [40] Jolivet, E. (1981). Central limit theorem and convergence of empirical processes for stationary point processes. In Point Processes and Queuing Problems (Colloq., Keszthely, 1978). Colloq. Math. Soc. János Bolyai 24 117–161. North-Holland, Amsterdam.
  • [41] Kallenberg, O. (1978). On conditional intensities of point processes. Z. Wahrsch. Verw. Gebiete 41 205–220.
  • [42] Kallenberg, O. (1984). An informal guide to the theory of conditioning in point processes. Internat. Statist. Rev. 52 151–164.
  • [43] Kelly, F. P. and Ripley, B. D. (1976). A note on Strauss’s model for clustering. Biometrika 63 357–360.
  • [44] Kulldorff, M. (1999). Spatial scan statistics: Models, calculations, and applications. In Scan Statistics and Applications 303–322. Birkhäuser, Boston, MA.
  • [45] Kutoyants, Y. A. (1998). Statistical Inference for Spatial Poisson Processes. Lecture Notes in Statist. 134. Springer, New York.
  • [46] Last, G. and Penrose, M. (2011). Poisson process Fock space representation, chaos expansion and covariance inequalities. Probab. Theory Related Fields 150 663–690.
  • [47] Lawson, A. B. (1993). On the analysis of mortality events around a prespecified fixed point. J. Roy. Statist. Soc. Ser. A 156 363–377.
  • [48] Lotwick, H. W. and Silverman, B. W. (1982). Methods for analysing spatial processes of several types of points. J. Roy. Statist. Soc. Ser. B 44 406–413.
  • [49] Møller, J., Syversveen, A. R. and Waagepetersen, R. P. (1998). Log Gaussian Cox processes. Scand. J. Statist. 25 451–482.
  • [50] Møller, J. and Waagepetersen, R. P. (2004). Statistical Inference and Simulation for Spatial Point Processes. Monogr. Statist. Appl. Probab. 100. Chapman & Hall/CRC, Boca Raton, FL.
  • [51] Møller, J. and Waagepetersen, R. P. (2007). Modern statistics for spatial point processes. Scand. J. Statist. 34 643–684.
  • [52] Nguyen, X.-X. and Zessin, H. (1979). Integral and differential characterizations of the Gibbs process. Math. Nachr. 88 105–115.
  • [53] Numata, M. (1961). Forest vegetation in the vicinity of Choshi—Coastal flora and vegetation at Choshi, Chiba prefecture, IV (in Japanese). Bull. Choshi Mar. Lab. 3 28–48.
  • [54] Numata, M. (1964). Forest vegetation, particularly pine stands in the vicinity of Choshi—Flora and vegetation in Choshi, Chiba prefecture, VI (in Japanese). Bull. Choshi Mar. Lab. 6 27–37.
  • [55] Ogata, Y. and Tanemura, M. (1981). Estimation of interaction potentials of spatial point patterns through the maximum likelihood procedure. Ann. Inst. Statist. Math. 33 315–338.
  • [56] Ogata, Y. and Tanemura, M. (1986). Likelihood estimation of interaction potentials and external fields of inhomogeneous spatial point patterns. In Pacific Statistical Congress (I. S. Francis, B. J. F. Manly and F. C. Lam, eds.) 150–154. Elsevier, Amsterdam.
  • [57] Ohser, J. (1983). On estimators for the reduced second moment measure of point processes. Math. Operationsforsch. Statist. Ser. Statist. 14 63–71.
  • [58] Ohser, J. and Stoyan, D. (1981). On the second-order and orientation analysis of planar stationary point processes. Biometrical J. 23 523–533.
  • [59] Papangelou, F. (1974). The conditional intensity of general point processes and an application to line processes. Z. Wahrsch. Verw. Gebiete 28 207–226.
  • [60] Pregibon, D. (1982). Score tests in GLIM with applications. In GLIM 82: Proceedings of the International Conference on Generalized Linear Models. Lecture Notes in Statist. 14. Springer, New York.
  • [61] Radhakrishna Rao, C. (1948). Large sample tests of statistical hypotheses concerning several parameters with applications to problems of estimation. Proc. Cambridge Philos. Soc. 44 50–57.
  • [62] Rathbun, S. L. and Cressie, N. (1994). Asymptotic properties of estimators for the parameters of spatial inhomogeneous Poisson point processes. Adv. in Appl. Probab. 26 122–154.
  • [63] Ripley, B. D. (1976). The second-order analysis of stationary point processes. J. Appl. Probab. 13 255–266.
  • [64] Ripley, B. D. (1977). Modelling spatial patterns (with discussion). J. Roy. Statist. Soc. Ser. B 39 172–212.
  • [65] Ripley, B. D. (1988). Statistical Inference for Spatial Processes. Cambridge Univ. Press, Cambridge.
  • [66] Ripley, B. D. and Kelly, F. P. (1977). Markov point processes. J. Lond. Math. Soc. (2) 15 188–192.
  • [67] Schladitz, K. and Baddeley, A. J. (2000). A third order point process characteristic. Scand. J. Statist. 27 657–671.
  • [68] Silvapulle, M. J. (1996). A test in the presence of nuisance parameters. J. Amer. Statist. Assoc. 91 1690–1693.
  • [69] Stein, M. L. (1995). An approach to asymptotic inference for spatial point processes. Statist. Sinica 5 221–234.
  • [70] Stillinger, D. K., Stillinger, F. H., Torquato, S., Truskett, T. M. and Debenedetti, P. G. (2000). Triangle distribution and equation of state for classical rigid disks. J. Statist. Phys. 100 49–72.
  • [71] Stoyan, D., Kendall, W. S. and Mecke, J. (1987). Stochastic Geometry and Its Applications. Wiley, Chichester.
  • [72] Stoyan, D. and Stoyan, H. (1995). Fractals, Random Shapes and Point Fields. Wiley, Chichester.
  • [73] Strauss, D. J. (1975). A model for clustering. Biometrika 62 467–475.
  • [74] van Lieshout, M. N. M. (2000). Markov Point Processes and Their Applications. Imperial College Press, London.
  • [75] Wald, A. (1941). Some examples of asymptotically most powerful tests. Ann. Math. Statist. 12 396–408.
  • [76] Waller, L., Turnbull, B., Clark, L. C. and Nasca, P. (1992). Chronic Disease Surveillance and testing of clustering of disease and exposure: Application to leukaemia incidence and TCE-contaminated dumpsites in upstate New York. Environmetrics 3 281–300.
  • [77] Wang, P. C. (1985). Adding a variable in generalized linear models. Technometrics 27 273–276.
  • [78] Widom, B. and Rowlinson, J. S. (1970). New model for the study of liquid–vapor phase transitions. J. Chem. Phys. 52 1670–1684.
  • [79] Wu, L. (2000). A new modified logarithmic Sobolev inequality for Poisson point processes and several applications. Probab. Theory Related Fields 118 427–438.