Residual analysis for spatial point processes (with discussion)

Summary.  We define residuals for point process models fitted to spatial point pattern data, and we propose diagnostic plots based on them. The residuals apply to any point process model that has a conditional intensity; the model may exhibit spatial heterogeneity, interpoint interaction and dependence on spatial covariates. Some existing ad hoc methods for model checking (quadrat counts, scan statistic, kernel smoothed intensity and Berman's diagnostic) are recovered as special cases. Diagnostic tools are developed systematically, by using an analogy between our spatial residuals and the usual residuals for (non-spatial) generalized linear models. The conditional intensity λ plays the role of the mean response. This makes it possible to adapt existing knowledge about model validation for generalized linear models to the spatial point process context, giving recommendations for diagnostic plots. A plot of smoothed residuals against spatial location, or against a spatial covariate, is effective in diagnosing spatial trend or co-variate effects. Q – Q -plots of the residuals are effective in diagnosing interpoint interaction.     

Two-step estimation for inhomogeneous spatial point processes

Summary.  The paper is concerned with parameter estimation for inhomogeneous spatial point processes with a regression model for the intensity function and tractable second-order properties ( K -function). Regression parameters are estimated by using a Poisson likelihood score estimating function and in the second step minimum contrast estimation is applied for the residual clustering parameters. Asymptotic normality of parameter estimates is established under certain mixing conditions and we exemplify how the results may be applied in ecological studies of rainforests.     

Modelling columnarity of pyramidal cells in the human cerebral cortex

For modelling the location of pyramidal cells in the human cerebral cortex, we suggest a hierarchical point process in that exhibits anisotropy in the form of cylinders extending along the z-axis. The model consists first of a generalised shot noise Cox process for the xy-coordinates, providing cylindrical clusters, and next of a Markov random field model for the z-coordinates conditioned on the xy-coordinates, providing either repulsion, aggregation or both within specified areas of interaction. Several cases of these hierarchical point processes are fitted to two pyramidal cell data sets, and of these a final model allowing for both repulsion and attraction between the points seem adequate. We discuss how the final model relates to the so-called minicolumn hypothesis in neuroscience.     

Using Weight Functions in Spatial Point Pattern Analysis with Application to Plant Ecology Data

This paper presents procedures for percentile estimation in the three-parameter inverse Gaussian (IG3) and the two-parameter inverse Gaussian (IG2) distributions. All procedures require first the estimation of distribution parameters and second the computation of the desired quantile at the estimated parameters. Parameter estimation is accomplished by maximum likelihood (ML) or a mixed moments (MXM) method. A Newton-Rahpson (NR) procedure is used for inverting the CDF. Simulation and asymptotic results are given for the resulting estimators. Two sets of real data are used to illustrate the procedures.     

Globally intensity-reweighted estimators for K- and pair correlation functions

We introduce new estimators of the inhomogeneous K-function and the pair correlation function of a spatial point process as well as the cross K-function and the cross pair correlation function of a bivariate spatial point process under the assumption of second-order intensity-reweighted stationarity. These estimators rely on a ‘global’ normalisation factor which depends on an aggregation of the intensity function, while the existing estimators depend ‘locally’ on the intensity function at the individual observed points. The advantages of our new global estimators over the existing local estimators are demonstrated by theoretical considerations and a simulation study.     

Structured Spatio-Temporal Shot-Noise Cox Point Process Models, with a View to Modelling Forest Fires

Abstract.  Spatio-temporal Cox point process models with a multiplicative structure for the driving random intensity, incorporating covariate information into temporal and spatial components, and with a residual term modelled by a shot-noise process, are considered. Such models are flexible and tractable for statistical analysis, using spatio-temporal versions of intensity and inhomogeneous K -functions, quick estimation procedures based on composite likelihoods and minimum contrast estimation, and easy simulation techniques. These advantages are demonstrated in connection with the analysis of a relatively large data set consisting of 2796 days and 5834 spatial locations of fires. The model is compared with a spatio-temporal log-Gaussian Cox point process model, and likelihood-based methods are discussed to some extent.     

A spatial analysis of health and pharmaceutical firm survival

The presence of knowledge spillovers and shared human capital is at the heart of the Marhall–Arrow–Romer externalities hypothesis. Most of the earlier empirical contributions on knowledge externalities; however, considered data aggregated at a regional level so that conclusions are based on the arbitrary definition of jurisdictional spatial units: this is the essence of the so-called modifiable areal unit problem. A second limitation of these studies is constituted by the fact that, somewhat surprisingly, while concentrating on the effects of agglomeration on firm creation and growth, the literature has, conversely, largely ignored its effects on firm survival. The present paper aims at contributing to the existing literature by answering to some of the open methodological questions reconciling the literature of Cox proportional hazards model with that on point pattern and thus capturing the true nature of spatial information. We also present some empirical results based on Italian firm demography data collected and managed by the Italian National Institute of Statistics (ISTAT).