The biogeography of invasive alien plants in California: an application of GIS and spatial regression analysis

The spatial distribution of invasive alien plants has been poorly documented in California. However, with the increased availability of GIS software and spatially explicit data, the distribution of invasive alien plants can be explored. Using bioregions as defined in Hickman (1993 ), I compared the distribution of invasive alien plants (n = 78) and noninvasive alien plants (n = 1097). The distribution of both categories of alien plants was similar with the exception of a higher concentration of invasive alien plants in the North Coast bioregion. Spatial autocorrelation analysis using Moran's I indicated significant spatial dependence for both invasive and noninvasive alien plant species. I used both ordinary least squares (OLS) and spatial autoregressive (SAR) models to assess the relationship between alien plant species distribution and native plant species richness, road density, population density, elevation, area of sample unit, and precipitation. The OLS model for invasive alien plants included two significant effects; native plant species richness and elevation. The SAR model for invasive alien plants included three significant effects; elevation, road density, and native plant species richness. The SAR model for noninvasive alien plants resulted in the same significant effects as invasive alien plants. Both invasive and noninvasive alien plants are found in regions with low elevation, high road density, and high native‐plant species richness. This is in congruity with previous spatial pattern studies of alien plant species. However, the similarity in effects for both categories of alien plants alludes to the importance of autecological attributes, such as pollination system, dispersal system and differing responses to disturbance in the distribution of invasive plant species. In addition, this study emphasizes the critical importance of testing for spatial autocorrelation in spatial pattern studies and using SAR models when appropriate.     

Testing for spatial correlation in nonstationary binary data, with application to aberrant crypt foci in colon carcinogenesis

In an experiment to understand colon carcinogenesis, all animals were exposed to a carcinogen, with half the animals also being exposed to radiation. Spatially, we measured the existence of what are referred to as aberrant crypt foci (ACF), namely, morphologically changed colonic crypts that are known to be precursors of colon cancer development. The biological question of interest is whether the locations of these ACFs are spatially correlated: if so, this indicates that damage to the colon due to carcinogens and radiation is localized. Statistically, the data take the form of binary outcomes (corresponding to the existence of an ACF) on a regular grid. We develop score-type methods based upon the Matern and conditionally autoregressive (CAR) correlation models to test for the spatial correlation in such data, while allowing for nonstationarity. Because of a technical peculiarity of the score-type test, we also develop robust versions of the method. The methods are compared to a generalization of Moran's test for continuous outcomes, and are shown via simulation to have the potential for increased power. When applied to our data, the methods indicate the existence of spatial correlation, and hence indicate localization of damage.     

Spatio-temporal modeling of localized brain activity

Functional neuroimaging, including positron emission tomography (PET) and functional magnetic resonance imaging (fMRI), plays an important role in identifying specific brain regions associated with experimental stimuli or psychiatric disorders such as schizophrenia. PET and fMRI produce massive data sets that contain both temporal correlations from repeated scans and complex spatial correlations. Several methods exist for handling temporal correlations, some of which rely on transforming the response data to induce either a known or an independence covariance structure. Despite the presence of spatial correlations between the volume elements (voxels) comprising a brain scan, conventional methods perform voxel-by-voxel analyses of measured brain activity. We propose a two-stage spatio-temporal model for the estimation and testing of localized activity. Our second-stage model specifies a spatial auto-regression, capturing correlations within neural processing clusters defined by a data-driven cluster analysis. We use maximum likelihood methods to estimate parameters from our spatial autoregressive model. Our model protects against type-I errors, enables the detection of both localized and regional activations (including volume of interest effects), provides information on functional connectivity in the brain, and establishes a framework to produce spatially smoothed maps of distributed brain activity for each individual. We illustrate the application of our model using PET data from a study of working memory in individuals with schizophrenia.     

Parametric spatial cure rate models for interval-censored time-to-relapse data

Several recent papers (e.g., Chen, Ibrahim, and Sinha, 1999, Journal of the American Statistical Association 94, 909-919; Ibrahim, Chen, and Sinha, 2001a, Biometrics 57, 383-388) have described statistical methods for use with time-to-event data featuring a surviving fraction (i.e., a proportion of the population that never experiences the event). Such cure rate models and their multivariate generalizations are quite useful in studies of multiple diseases to which an individual may never succumb, or from which an individual may reasonably be expected to recover following treatment (e.g., various types of cancer). In this article we extend these models to allow for spatial correlation (estimable via zip code identifiers for the subjects) as well as interval censoring. Our approach is Bayesian, where posterior summaries are obtained via a hybrid Markov chain Monte Carlo algorithm. We compare across a broad collection of rather high-dimensional hierarchical models using the deviance information criterion, a tool recently developed for just this purpose. We apply our approach to the analysis of a smoking cessation study where the subjects reside in 53 southeastern Minnesota zip codes. In addition to the usual posterior estimates, our approach yields smoothed zip code level maps of model parameters related to the relapse rates over time and the ultimate proportion of quitters (the cure rates).     

Effects of residual smoothing on the posterior of the fixed effects in disease-mapping models

Disease-mapping models for areal data often have fixed effects to measure the effect of spatially varying covariates and random effects with a conditionally autoregressive (CAR) prior to account for spatial clustering. In such spatial regressions, the objective may be to estimate the fixed effects while accounting for the spatial correlation. But adding the CAR random effects can cause large changes in the posterior mean and variance of fixed effects compared to the nonspatial regression model. This article explores the impact of adding spatial random effects on fixed effect estimates and posterior variance. Diagnostics are proposed to measure posterior variance inflation from collinearity between the fixed effect covariates and the CAR random effects and to measure each region's influence on the change in the fixed effect's estimates by adding the CAR random effects. A new model that alleviates the collinearity between the fixed effect covariates and the CAR random effects is developed and extensions of these methods to point-referenced data models are discussed.     

Modeling longitudinal spatial periodontal data: a spatially adaptive model with tools for specifying priors and checking fit

Summary .   Attachment loss (AL), the distance down a tooth's root that is no longer attached to surrounding bone by periodontal ligament, is a common measure of periodontal disease. In this article, we develop a spatiotemporal model to monitor the progression of AL. Our model is an extension of the conditionally autoregressive (CAR) prior, which spatially smooths estimates toward their neighbors. However, because AL often exhibits a burst of large values in space and time, we develop a nonstationary spatiotemporal CAR model that allows the degree of spatial and temporal smoothing to vary in different regions of the mouth. To do this, we assign each AL measurement site its own set of variance parameters and spatially smooth the variances with spatial priors. We propose a heuristic to measure the complexity of the site-specific variances, and use it to select priors that ensure parameters in the model are well identified. In data from a clinical trial, this model improves the fit compared to the usual dynamic CAR model for 90 of 99 patients' AL measurements.     

A Parameter-Elimination Method for Nonlinear Regression with Linear Parameters and Autocorrelated Errors

In this study, we are interested in the problem of estimating the parameters in a nonlinear regression model when the error terms are correlated. Throughout this work, we restrict ourselves to the special case when the error terms follow a pth order stationary autoregressive model (AR(p)). Following the idea of LAWTON and SYLVESTRE (1971) and GALLANT and GOEBEL (1976), a parameter-elimination method is proposed, which has the advantages that it is not sensitive to the initial values and convergence of the procedure may be more stable because of the reduced dimension of the problem. The parameter-elimination method is compared with the methods by GALLANT and GOEBEL (1976) and GLASBEY (1980) by Monte Carlo Simulation, and the results of applying the first two methods to the real data obtained from the Environmental Protection Administration of the Executive Yuan of the Republic of China are presented.