Modelling Spatial Intensity for Replicated Inhomogeneous Point Patterns in Brain Imaging
- 14 April 2004
- journal article
- research article
- Published by Oxford University Press (OUP) in Journal of the Royal Statistical Society Series B: Statistical Methodology
- Vol. 66 (2) , 429-446
- https://doi.org/10.1046/j.1369-7412.2003.05285.x
Abstract
Summary: Pharmacological experiments in brain microscopy study patterns of cellular activ- ation in response to psychotropic drugs for connected neuroanatomic regions. A typical ex- perimental design produces replicated point patterns having highly complex spatial variability. Modelling this variability hierarchically can enhance the inference for comparing treatments. We propose a semiparametric formulation that combines the robustness of a nonparametric kernel method with the efficiency of likelihood-based parameter estimation. In the convenient framework of a generalized linear mixed model, we decompose pattern variation by kriging the intensities of a hierarchically heterogeneous spatial point process. This approximation entails discretizing the inhomogeneous Poisson likelihood by Voronoi tiling of augmented point patterns. The resulting intensity-weighted log-linear model accommodates spatial smoothing through a reduced rank penalized linear spline. To correct for anatomic distortion between subjects, we interpolate point locations via an isomorphic mapping so that smoothing occurs relative to common neuroanatomical atlas co-ordinates. We propose a criterion for choosing the degree and spatial locale of smoothing based on truncating the ordered set of smoothing covariates to minimize residual extra-dispersion. Additional spatial covariates, experimental design factors, hierarchical random effects and intensity functions are readily accommodated in the linear predictor, enabling comprehensive analyses of the salient properties underlying replicated point patterns. We illustrate our method through application to data from a novel study of drug effects on neuronal activation patterns in the brain of rats.Funding Information
- National Institutes of Health (AI51951, ES10844, NS37483, MH60450)
This publication has 33 references indexed in Scilit:
- Geoadditive ModelsJournal of the Royal Statistical Society Series C: Applied Statistics, 2003
- Practical Maximum Pseudolikelihood for Spatial Point PatternsAustralian & New Zealand Journal of Statistics, 2000
- A comparison between parametric and non-parametric approaches to the analysis of replicated spatial point patternsAdvances in Applied Probability, 2000
- Log Gaussian Cox ProcessesScandinavian Journal of Statistics, 1998
- Poisson/gamma random field models for spatial statisticsBiometrika, 1998
- Maximum Likelihood Algorithms for Generalized Linear Mixed ModelsJournal of the American Statistical Association, 1997
- Algorithm 751: TRIPACKACM Transactions on Mathematical Software, 1996
- Approximating Point Process Likelihoods with GLIMJournal of the Royal Statistical Society Series C: Applied Statistics, 1992
- Analysis of Variance for Replicated Spatial Point Patterns in Clinical NeuroanatomyJournal of the American Statistical Association, 1991
- Longitudinal data analysis using generalized linear modelsBiometrika, 1986