Second-order analysis of inhomogeneous spatial point processes using case-control data

Methods for the statistical analysis of stationary spatial point process data are now well established, methods for nonstationary processes less so. One of many sources of nonstationary point process data is a case-control study in environmental epidemiology. In that context, the data consist of a realization of each of two spatial point processes representing the locations, within a specified geographical region, of individual cases of a disease and of controls drawn at random from the population at risk. In this article, we extend work by Baddeley, Moller, and Waagepetersen (2000, Statistica Neerlandica 54, 329-350) concerning estimation of the second-order properties of a nonstationary spatial point process. First, we show how case-control data can be used to overcome the problems encountered when using the same data to estimate both a spatially varying intensity and second-order properties. Second, we propose a semiparametric method for adjusting the estimate of intensity so as to take account of explanatory variables attached to the cases and controls. Our primary focus is estimation, but we also propose a new test for spatial clustering that we show to be competitive with existing tests. We describe an application to an ecological study in which juvenile and surviving adult trees assume the roles of controls and cases.