Hierarchical modeling for extreme values observed over space and time

We propose a hierarchical modeling approach for explaining a collection of spatially referenced time series of extreme values. We assume that the observations follow generalized extreme value (GEV) distributions whose locations and scales are jointly spatially dependent where the dependence is captured using multivariate Markov random field models specified through coregionalization. In addition, there is temporal dependence in the locations. There are various ways to provide appropriate specifications; we consider four choices. The models can be fitted using a Markov Chain Monte Carlo (MCMC) algorithm to enable inference for parameters and to provide spatio–temporal predictions. We fit the models to a set of gridded interpolated precipitation data collected over a 50-year period for the Cape Floristic Region in South Africa, summarizing results for what appears to be the best choice of model.

Using data augmentation via the Gibbs Sampler to incorporate missing covariate structure in linear models for ecological assessments

Missing covariate values in linear regression models can be an important problem facing environmental researchers. Existing missing value treatment methods such as Multiple Imputation (MI), the EM algorithm and Data Augmentation (DA) have the assumption that both observed and unobserved data come from the same distribution, most commonly a multivariate normal or a conditionally multivariate normal family. These methods do try to incorporate the missing data mechanism and rely on the assumption of Missing At Random (MAR). We present a DA method which does not rely on the MAR assumption and can model missing data mechanisms and covariate structure. This method utilizes the Gibbs Sampler as a tool for incorporating these structures and mechanisms. We apply this method to an ecological data set that relates fish condition to environmental variables. Notice that the presented DA method detects relationships that are not detected when other missing data methods are employed.

Prospective spatial prediction of infectious disease: experience of New York State (USA) with West Nile Virus and proposed directions for improved surveillance

Infectious disease surveillance has become an international top priority due to the perceived risk of bioterrorism. This is driving the improvement of real-time geo-spatial surveillance systems for monitoring disease indicators, which is expected to have many benefits beyond detecting a bioterror event. West Nile Virus surveillance in New York State (USA) is highlighted as a working system that uses dead American Crows (Corvus brachyrhynchos) to prospectively indicate viral activity prior to human onset. A cross-disciplinary review is then presented to argue that this system, and infectious disease surveillance in general, can be improved by complementing spatial cluster detection of an outcome variable with predictive “risk mapping” that incorporates spatiotemporal data on the environment, climate and human population through the flexible class of generalized linear mixed models.

Flexible hierarchical mark-recapture modeling for open populations using WinBUGS

Hierarchical mark-recapture models offer three advantages over classical mark-recapture models: (i) they allow expression of complicated models in terms of simple components; (ii) they provide a convenient way of modeling missing data and latent variables in a way that allows expression of relationships involving latent variables in the model; (iii) they provide a convenient way of introducing parsimony into models involving many nuisance parameters. Expressing models using the complete data likelihood we show how many of the standard mark-recapture models for open populations can be readily fitted using the software WinBUGS. We include examples that illustrate fitting the Cormack–Jolly–Seber model, multi-state and multi-event models, models including auxiliary data, and models including density dependence.

A spatial clustering perspective on autocorrelation and regionalization

We revisit one of the classical problems in geography and cartography where multiple observations on a lattice (N) need to be grouped into many fewer regions (G), especially when this number of desired regions is unknown a priori. Since an optimization through all possible aggregations is not feasible, a hierarchical classification scheme is proposed with an objective function sensitive to spatial pattern. The objective function to be minimized during the assignment of observations to regions (classification) consists of two terms: the first characterizes accuracy and the second, model complexity. For the latter, we introduce a spatial measure that characterizes the number of homogeneous patches rather than the usual number of classes. A simulation study shows that such a classification procedure is less sensitive to random and spatially correlated error (noise) than non-spatial classification. We also show that for conditional autoregressive error (noise) fields the optimal partitioning is the one that has the highest within-units generalized Moran coefficient. The classifier is implemented in ArcView to demonstrate both a socio-economic and an environmental application to illustrate some potential applications.

A note on data augmentation and iterative simulation for survival rate estimation

In this paper we examine the use of data augmentation techniques for simplifying iterative simulation in the context of both Bayesian and classical statistical inference for survival rate estimation. We examine two distinct model families common in population ecology to illustrate our ideas, ring-recovery models and capture–recapture models, and we present the computational advantage of this approach. We discuss also the fact that problems associated with identifiability in the classical framework can be overcome using data augmentation, but highlight the dangers in doing so under both inferential paradigms.

A note on non-parametric and semi-parametric modeling of wildfire hazard in Los Angeles County, California

This paper explores the use of, and problems that arise in, kernel smoothing and parametric estimation of the relationships between wildfire incidence and various meteorological variables. Such relationships may be treated as components in separable point process models for wildfire activity. The resulting models can be used for comparative purposes in order to assess the predictive performance of the Burning Index.

Fiat boundaries: some implications for interpretation,decision-support,and multi-temporal analysis

Polygon-based thematic maps can be composed of boundaries that exist by definition—i.e., bona fide boundaries—or those that exist relative to a specific interpretation of a spatial phenomenon—i.e., fiat boundaries. The construction of maps composed of fiat boundaries is usually based on a subjective interpretive methodology that is affected by the data used to construct the map and the minimum mapping unit employed. That fiat boundaries are not the same as bona fide boundaries affects their use in computer-based spatial decision support tools. This is discussed both in terms of an analysis conducted at one specific moment, and in respect to increasingly common multi-temporal analysis.

Rarefaction and nonrandom spatial dispersion patterns

Rarefaction estimates how many species are expected in a random sample of individuals from a larger collection and allows meaningful comparisons among collections of different sizes. It assumes random spatial dispersion. However, two common dispersion patterns, within-species clumping and segregation among species, can cause rarefaction to overestimate the species richness of a smaller continuous area. We use field studies and computer simulations to determine (1) how robust rarefaction is to nonrandom spatial dispersion and (2) whether simple measures of spatial autocorrelation can predict the bias in rarefaction estimates. Rarefaction does not estimate species richness accurately for many communities, especially at small sample sizes. Measures of spatial autocorrelation of the more abundant species do not reliably predict amount of bias. Survey sites should be standardized to equal-sized areas before sampling. When sites are of equal area but differ in number of individuals sampled, rarefaction can standardize collections. When communities are sampled from different-sized areas, the mean and confidence intervals of species accumulation curves allow more meaningful comparisons among sites. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Seals at sea: modelling seal distribution in the German bight based on aerial survey data

The Wadden Sea is an important habitat for harbour seals and grey seals. They regularly haul-out on sandbanks and islands along the coast. Comparably little is known about the time seals spend at sea and how they use the remainder of the North Sea. Yet, human activity in offshore waters is increasing and information on seal distribution in the North Sea is crucial for conservation and management. Aerial line transect surveys were conducted in the German bight from 2002 to 2007 to investigate the distribution and abundance of marine mammals. Distance sampling methodology was combined with density surface modelling for a spatially explicit analysis of seal distribution in the German North Sea. Depth and distance to coast were found to be relevant predictor variables for seal density. Density surface modelling allowed for a depiction of seal distribution in the study area as well as an abundance estimate. This is the first study to use aerial survey data to develop a density surface model (DSM) for a spatially explicit distribution estimate of seals at sea.

---

Heteroscedastic semiparametric models for domestic water consumption aggregated data

Heteroscedastic additive and multiplicative models are proposed to disaggregate household data on water consumption from Athens and provide individual consumption estimates. The models adjust for heteroscedasticity assuming that variances relate to covariates. Household characteristics that can influence consumption are also included into models in order to allow for a clearer measurement of individual characteristics effects. Estimation is accomplished through a penalized least squares approach. The method is applied to a sample of real data related to domestic water consumption in Athens. The results show a greater consumption of water for males while the single-female households are these that use the lowest quantities of water. The consumption curves by age and gender are constructed presenting differences between the two sexes.

The fire frequency analysis branch of the pyrostatistics tree: sampling decisions and censoring in fire interval data

Statistical characterization of past fire regimes is important for both the ecology and management of fire-prone ecosystems. Survival analysis—or fire frequency analysis as it is often called in the fire literature—has increasingly been used over the last few decades to examine fire interval distributions. These distributions can be generated from a variety of sources (e.g., tree rings and stand age patterns), and analysis typically involves fitting the Weibull model. Given the widespread use of fire frequency analysis and the increasing availability of mapped fire history data, our goal has been to review and to examine some of the issues faced in applying these methods in a spatially explicit context. In particular, through a case study on the massive Cedar Fire in 2003 in southern California, we examine sensitivities of parameter estimates to the spatial resolution of sampling, point- and area-based methods for assigning sample values, current age surfaces versus historical intervals in generating distributions, and the inclusion of censored (i.e., incomplete) observations. Weibull parameter estimates were found to be roughly consistent with previous fire frequency analyses for shrublands (i.e., median age at burning of ~30–50 years and relatively low age dependency). Results indicate, however, that the inclusion or omission of censored observations can have a substantial effect on parameter estimates, far more than other decisions about specifics of sampling.

Spatial structure effects on the detection of patches boundaries using local operators

Landscapes exhibit various degrees of spatial heterogeneity according to the differential intensity and interactions among processes and disturbances that they are subjected to. The management of these spatially dynamical landscapes requires that we can accurately map them and monitor the evolution of their spatial arrangement through time. Such a mapping requires first the delineation of various spatial features present in the landscape such as patches and their boundaries. However, there are several environmental (spatial variability) as well as technical (spatial resolution) factors that impair our ability to accurately delineate patches and their boundaries as polygons. Here, we investigate how the spatial structure and spatial resolution of the data affect the accuracy of detecting patches and their boundaries over simulated landscapes and real data. Simulated landscapes consisted of two patches with parameterized spatial properties (patches’ level of spatial autocorrelation, mean value and variance) separated by a boundary of known location. Real data allowed the investigation of a more complex landscape where there is a known transition between two forest domains with unknown spatial properties. Boundary locations are defined using the lattice-wombling edge detector at various aggregation levels and the degree of patch homogeneity is determined using Getis-Ord’s G*. Results show that boundary detection using a local edge detector is greatly affected by the spatial conditions of the data, namely variance, abruptness of the spatial gradient between two patches and patches’ level of spatial autocorrelation. They also suggest that data aggregation is not a panacea for bringing out the ecological process creating the patches and that indicators derived from local measures of spatial association can be complementary tools for analysing spatial structures affecting boundary delineation.

Improved strategies for coverage estimation by using replicated line-intercept sampling

Coverage, i.e., the area covered by the target attribute in the study region, is a key parameter in many surveys. Coverage estimation is usually performed by adopting a replicated protocol based on line-intercept sampling coupled with a suitable linear homogeneous estimator. Since coverage is a parameter which may be interestingly represented as the integral of a suitable function, improved Monte Carlo strategies for implementing the replicated protocol are introduced in order to achieve estimators with small variance rates. In addition, new specific theoretical results on Monte Carlo integration methods are given to deal with the integrand functions arising in the special coverage estimation setting.

Spatial pattern and persistence of historical fire boundaries in southern interior British Columbia

Recent ecosystem and fire management research aims to quantify, model and understand historical fire disturbances focusing on their frequency, size and distribution. Less attention, however, has been paid to fire boundaries and their location on the landscape. Our study presents a spatial method to quantify the location, pattern and persistence of historical fire boundaries using tree ring fire scar data in the lower Stein watershed (British Columbia). Data from Pinus ponderosa and Pseudotsuga menziesii collected in 35 one-hectare plots over a 412-hectare study area were analyzed for the period between 1879 and 1947 using local spatial statistics and boundary detection techniques. Results of the analysis using local spatial statistic Moran’s I showed significant clustering of boundaries near topographic breaks. To determine and test whether fire boundaries between plots were persistent, we used boundary detection methods and developed a spatially restricted randomization test. The results revealed that out of 86 possible boundary links, 8 were significantly persistent (P < 0.025) while another 8 were significantly rare (P < 0.025). These spatial methods can help determine the historical spatial configuration of persistent boundaries and can be used to maintain natural forest dynamics.

Behavior of yellowfin (<Emphasis Type="Italic">Thunnus albacares</Emphasis>) and bigeye (<Emphasis Type="Italic">T. obesus</Emphasis>) tuna in a network of fish aggregating devices (FADs)

The influence of multiple anchored fish aggregating devices (FADs) on the spatial behavior of yellowfin (Thunnus albacares) and bigeye tuna (T. obesus) was investigated by equipping all thirteen FADs surrounding the island of Oahu (HI, USA) with automated sonic receivers (“listening stations”) and intra-peritoneally implanting individually coded acoustic transmitters in 45 yellowfin and 12 bigeye tuna. Thus, the FAD network became a multi-element passive observatory of the residence and movement characteristics of tuna within the array. Yellowfin tuna were detected within the FAD array for up to 150 days, while bigeye tuna were only observed up to a maximum of 10 days after tagging. Only eight yellowfin tuna (out of 45) and one bigeye tuna (out of 12) visited FADs other than their FAD of release. Those nine fish tended to visit nearest neighboring FADs and, in general, spent more time at their FAD of release than at the others. Fish visiting the same FAD several times or visiting other FADs tended to stay longer in the FAD network. A majority of tagged fish exhibited some synchronicity when departing the FADs but not all tagged fish departed a FAD at the same time: small groups of tagged fish left together while others remained. We hypothesize that tuna (at an individual or collective level) consider local conditions around any given FAD to be representative of the environment on a larger scale (e.g., the entire island) and when those conditions become unfavorable the tuna move to a completely different area. Thus, while the anchored FADs surrounding the island of Oahu might concentrate fish and make them more vulnerable to fishing, at a meso-scale they might not entrain fish longer than if there were no (or very few) FADs in the area. At the existing FAD density, the ‘island effect’ is more likely to be responsible for the general presence of fish around the island than the FADs. We recommend further investigation of this hypothesis.

---

Hierarchical spatial point process analysis for a plant community with high biodiversity

A complex multivariate spatial point pattern of a plant community with high biodiversity is modelled using a hierarchical multivariate point process model. In the model, interactions between plants with different post-fire regeneration strategies are of key interest. We consider initially a maximum likelihood approach to inference where problems arise due to unknown interaction radii for the plants. We next demonstrate that a Bayesian approach provides a flexible framework for incorporating prior information concerning the interaction radii. From an ecological perspective, we are able both to confirm existing knowledge on species’ interactions and to generate new biological questions and hypotheses on species’ interactions.

Maximum likelihood fitting of the Poisson lognormal distribution

In this paper some properties and analytic expressions regarding the Poisson lognormal distribution such as moments, maximum likelihood function and related derivatives are discussed. The author provides a sharp approximation of the integrals related to the Poisson lognormal probabilities and analyzes the choice of the initial values in the fitting procedure. Based on these he describes a new procedure for carrying out the maximum likelihood fitting of the truncated Poisson lognormal distribution. The method and results are illustrated on real data. The computer program for calculations is freely available.