Comparison of Two Methods for Caluculating the Partition Functions of Various Spatial Statistical Models

Likelihood computation in spatial statistics requires accurate and efficient calculation of the normalizing constant (i.e. partition function) of the Gibbs distribution of the model. Two available methods to calculate the normalizing constant by Markov chain Monte Carlo methods are compared by simulation experiments for an Ising model, a Gaussian Markov field model and a pairwise interaction point field model.

Distribution of paracoccidioidomycosis: determination of ecologic correlates through spatial analyses

Paracoccidioidomycosis (PCM) is endemic in Latin America and in countries like Brazil it carries a high mortality rate. The fungus' habitat has not been precisely determined. The present study aims to identify ecologic correlates based on PCM distribution in a hyper-endemic area in southeastern Brazil. The Geographic Information System (GIS) and spatial statistics were used to associate environmental attributes, human population density and, PCM distribution. By means of the Pearson r correlation coefficient, the highest statistically significant associations with prevalence density were the percent area (by county) of: basaltic rocks (r = 0.63; P < 0.0001), Podzolic soils (r = - 0.48; P < 0.001), Latosol soils (r = 0.40; P < 0.01), mean annual precipitation between 1500 and 1600 mm (r = 0.46; P < 0.001) and, mean precipitation during the wet season between 940 and 1040 mm (r = - 0.44; P < 0.01). Soil texture and precipitation analyzed together reached r = 0.61 (P < 0.000002) for fine-textured soils with annual precipitation above 1400 mm. Environmental correlates indicate that moisture availability plays an important role in PCM distribution.

Spatial distribution of canine Visceral Leishmaniasis in Ilha Solteira, São Paulo, Brazil

A Leishmaniose Visceral (LV) é causada por protozoários do gênero Leishmania e transmitida por flebotomíneos do gênero Lutzomyia, os quais vêm adaptando-se ao ambiente peridomiciliar, onde o cão é sua principal fonte de alimento, aumentando assim o risco de casos em humanos. Neste trabalho, foram utilizadas técnicas de geoprocessamento e de estastística espacial como contribuição à compreensão da dinâmica epidemiológica da LV na área urbana de Ilha Solteira-SP.

Spatial pattern detection modeling of thrips (Thrips tabaci) on onion fields

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Influence diagnostics in Gaussian spatial linear models

Spatial linear models have been applied in numerous fields such as agriculture, geoscience and environmental sciences, among many others. Spatial dependence structure modelling, using a geostatistical approach, is an indispensable tool to estimate the parameters that define this structure. However, this estimation may be greatly affected by the presence of atypical observations in the sampled data. The purpose of this paper is to use diagnostic techniques to assess the sensitivity of the maximum-likelihood estimators, covariance functions and linear predictor to small perturbations in the data and/or the spatial linear model assumptions. The methodology is illustrated with two real data sets. The results allowed us to conclude that the presence of atypical values in the sample data have a strong influence on thematic maps, changing the spatial dependence structure.

Computational Techniques for Spatial Logistic Regression with Large Datasets

In epidemiological work, outcomes are frequently non-normal, sample sizes may be large, and effects are often small. To relate health outcomes to geographic risk factors, fast and powerful methods for fitting spatial models, particularly for non-normal data, are required. We focus on binary outcomes, with the risk surface a smooth function of space. We compare penalized likelihood models, including the penalized quasi-likelihood (PQL) approach, and Bayesian models based on fit, speed, and ease of implementation. A Bayesian model using a spectral basis representation of the spatial surface provides the best tradeoff of sensitivity and specificity in simulations, detecting real spatial features while limiting overfitting and being more efficient computationally than other Bayesian approaches. One of the contributions of this work is further development of this underused representation. The spectral basis model outperforms the penalized likelihood methods, which are prone to overfitting, but is slower to fit and not as easily implemented. Conclusions based on a real dataset of cancer cases in Taiwan are similar albeit less conclusive with respect to comparing the approaches. The success of the spectral basis with binary data and similar results with count data suggest that it may be generally useful in spatial models and more complicated hierarchical models.

Chapter 3 Industrial Location in Cambodia: Industrial Policies and Economic Statistics

Cambodia has experienced high economic growth in the last decade. Because most of its industries were destroyed during the Pol Pot regime and civil war, in the last 20 years the country has been working hard to liberalize its economy to attract foreign investors With its efforts to join the regional and international community and with changes in the international trade environment, Cambodia started to grow its economy in the late 1990s. Now, in the early 21st century, the Cambodian economy seems to be prepared to take off. We can observe a kind of industrial agglomeration occurring, even though still at a small scale. In this paper, first, I will review the history of Cambodia’s economic development since the late 1980s. Second, I will examine the economic policies, laws, rules, and other environmental factors which have influenced industrial development and industrial location in Cambodia. Third, I will introduce industrial location in the late 2000s. Lastly, I will introduce some statistical data for the future analysis of industrial location in Cambodia.

Chapter 4 A Review of the Availability of and Studies on Spatial Economic Data in Laos

Spatial data are being increasingly used in a wide range of disciplines, a fact that is clearly reflected in the recent trend to add spatial dimensions to the conventional social sciences. Economics is by no means an exception. On one hand, spatial data are indispensable to many branches of economics such as economic geography, new economic geography, or spatial economics. On the other hand, macroeconomic data are becoming available at more and more micro levels, so that academics and analysts take it for granted that they are available not only for an entire country, but also for more detailed levels (e.g. state, province, and even city). The term ‘spatial economics data’ as used in this report refers to any economic data that has spatial information attached. This spatial information can be the coordinates of a location at best or a less precise place name as is used to describe administrative units. Obviously, the latter cannot be used without a map of corresponding administrative units. Maps are therefore indispensible to the analysis of spatial economic data without absolute coordinates. The aim of this report is to review the availability of spatial economic data that pertains specifically to Laos and academic studies conducted on such data up to the present. In regards to the availability of spatial economic data, efforts have been made to identify not only data that has been made available as geographic information systems (GIS) data, but also those with sufficient place labels attached. The rest of the report is organized as follows. Section 2 reviews the maps available for Laos, both in hard copy and editable electronic formats. Section 3 summarizes the spatial economic data available for Laos at the present time, and Section 4 reviews and categorizes the many economic studies utilizing these spatial data. Section 5 give examples of some of the spatial industrial data collected for this research. Section 6 provides a summary of the findings and gives some indication of the direction of the final report due for completion in fiscal 2010.

Competition may explain the fine-scale spatial patterns and genetic structure of two co-occurring plant congeners

1. The spatial distribution of individual plants within a population and the population’s genetic structure are determined by several factors, like dispersal, reproduction mode or biotic interactions. The role of interspecific interactions in shaping the spatial genetic structure of plant populations remains largely unknown. 2. Species with a common evolutionary history are known to interact more closely with each other than unrelated species due to the greater number of traits they share. We hypothesize that plant interactions may shape the fine genetic structure of closely related congeners. 3. We used spatial statistics (georeferenced design) and molecular techniques (ISSR markers) to understand how two closely related congeners, Thymus vulgaris (widespread species) and T. loscosii (narrow endemic) interact at the local scale. Specific cover, number of individuals of both study species and several community attributes were measured in a 10 × 10 m plot. 4. Both species showed similar levels of genetic variation, but differed in their spatial genetic structure. Thymus vulgaris showed spatial aggregation but no spatial genetic structure, while T. loscosii showed spatial genetic structure (positive genetic autocorrelation) at short distances. The spatial pattern of T. vulgaris’ cover showed significant dissociation with that of T. loscosii. The same was true between the spatial patterns of the cover of T. vulgaris and the abundance of T. loscosii and between the abundance of each species. Most importantly, we found a correlation between the genetic structure of T. loscosii and the abundance of T. vulgaris: T. loscosii plants were genetically more similar when they were surrounded by a similar number of T. vulgaris plants. 5. Synthesis. Our results reveal spatially complex genetic structures of both congeners at small spatial scales. The negative association among the spatial patterns of the two species and the genetic structure found for T. loscosii in relation to the abundance of T. vulgaris indicate that competition between the two species may account for the presence of adapted ecotypes of T. loscosii to the abundance of a competing congeneric species. This suggests that the presence and abundance of close congeners can influence the genetic spatial structure of plant species at fine scales.

Ross River virus disease clusters and spatial relationship with mosquito biting exposure in Redland Shire, southern Queensland, Australia

The spatial heterogeneity in the risk of Ross River virus (family Togaviridae, genus Alphavirus, RRV) disease, the most common mosquito-borne disease in Australia, was examined in Redland Shire in southern Queensland, Australia. Disease cases, complaints from residents of intense mosquito biting exposure, and human population data were mapped using a geographic information system. Surface maps of RRV disease age-sex standardized morbidity ratios and mosquito biting complaint morbidity ratios were created. To determine whether there was significant spatial variation in disease and complaint patterns, a spatial scan analysis method was used to test whether the number of cases and complaints was distributed according to underlying population at risk. Several noncontiguous areas in proximity to productive saline water habitats of Aedes vigilax (Skuse), a recognized vector of RRV, had higher than expected numbers of RRV disease cases and complaints. Disease rates in human populations in areas which had high numbers of adult Ae. vigilax in carbon dioxide- and octenol-baited light traps were up to 2.9 times those in areas that rarely had high numbers of mosquitoes. It was estimated that targeted control of adult Ae. vigilax in these high-risk areas could potentially reduce the RRV disease incidence by an average of 13.6%. Spatial correlation was found between RRV disease risk and complaints from residents of mosquito biting. Based on historical patterns of RRV transmission throughout Redland Shire and estimated future human population growth in areas with higher than average RRV disease incidence, it was estimated that RRV incidence rates will increase by 8% between 2001 and 2021. The use of arbitrary administrative areas that ranged in size from 4.6 to 318.3 km2, has the potential to mask any small scale heterogeneity in disease patterns. With the availability of georeferenced data sets and high-resolution imagery, it is becoming more feasible to undertake spatial analyses at relatively small scales.

Higher-order co-occurrences for exploratory point pattern analysis and decision tree clustering on spatial data

Analyzing geographical patterns by collocating events, objects or their attributes has a long history in surveillance and monitoring, and is particularly applied in environmental contexts, such as ecology or epidemiology. The identification of patterns or structures at some scales can be addressed using spatial statistics, particularly marked point processes methodologies. Classification and regression trees are also related to this goal of finding "patterns" by deducing the hierarchy of influence of variables on a dependent outcome. Such variable selection methods have been applied to spatial data, but, often without explicitly acknowledging the spatial dependence. Many methods routinely used in exploratory point pattern analysis are2nd-order statistics, used in a univariate context, though there is also a wide literature on modelling methods for multivariate point pattern processes. This paper proposes an exploratory approach for multivariate spatial data using higher-order statistics built from co-occurrences of events or marks given by the point processes. A spatial entropy measure, derived from these multinomial distributions of co-occurrences at a given order, constitutes the basis of the proposed exploratory methods. © 2010 Elsevier Ltd.

Higher-order co-occurrences for exploratory point pattern analysis and decision tree clustering on spatial data

Analyzing geographical patterns by collocating events, objects or their attributes has a long history in surveillance and monitoring, and is particularly applied in environmental contexts, such as ecology or epidemiology. The identification of patterns or structures at some scales can be addressed using spatial statistics, particularly marked point processes methodologies. Classification and regression trees are also related to this goal of finding "patterns" by deducing the hierarchy of influence of variables on a dependent outcome. Such variable selection methods have been applied to spatial data, but, often without explicitly acknowledging the spatial dependence. Many methods routinely used in exploratory point pattern analysis are2nd-order statistics, used in a univariate context, though there is also a wide literature on modelling methods for multivariate point pattern processes. This paper proposes an exploratory approach for multivariate spatial data using higher-order statistics built from co-occurrences of events or marks given by the point processes. A spatial entropy measure, derived from these multinomial distributions of co-occurrences at a given order, constitutes the basis of the proposed exploratory methods. © 2010 Elsevier Ltd.

GIS technologies and spatial organization patterns among hunter-gatherers: an archaeological case study from Southern Patagonia (southern South America)

This paper discusses some aspects of hunter-gatherer spatial organization in southern South Patagonia, in later times to 10,000 cal yr BP. Various methods of spatial analysis, elaborated with a Geographic Information System (GIS) were applied to the distributional pattern of archaeological sites with radiocarbon dates. The shift in the distributional pattern of chronological information was assessed in conjunction with other lines of evidence within a biogeographic framework. Accordingly, the varying degrees of occupation and integration of coastal and interior spaces in human spatial organization are explained in association with the adaptive strategies hunter-gatherers have used over time. Both are part of the same human response to changes in risk and uncertainty variability in the region in terms of resource availability and environmental dynamics.

Hierarchical models for 2D presence/absence data having ambiguous zeroes: With a biogeographical case study on dingo behaviour

This dissertation is primarily an applied statistical modelling investigation, motivated by a case study comprising real data and real questions. Theoretical questions on modelling and computation of normalization constants arose from pursuit of these data analytic questions. The essence of the thesis can be described as follows. Consider binary data observed on a two-dimensional lattice. A common problem with such data is the ambiguity of zeroes recorded. These may represent zero response given some threshold (presence) or that the threshold has not been triggered (absence). Suppose that the researcher wishes to estimate the effects of covariates on the binary responses, whilst taking into account underlying spatial variation, which is itself of some interest. This situation arises in many contexts and the dingo, cypress and toad case studies described in the motivation chapter are examples of this. Two main approaches to modelling and inference are investigated in this thesis. The first is frequentist and based on generalized linear models, with spatial variation modelled by using a block structure or by smoothing the residuals spatially. The EM algorithm can be used to obtain point estimates, coupled with bootstrapping or asymptotic MLE estimates for standard errors. The second approach is Bayesian and based on a three- or four-tier hierarchical model, comprising a logistic regression with covariates for the data layer, a binary Markov Random field (MRF) for the underlying spatial process, and suitable priors for parameters in these main models. The three-parameter autologistic model is a particular MRF of interest. Markov chain Monte Carlo (MCMC) methods comprising hybrid Metropolis/Gibbs samplers is suitable for computation in this situation. Model performance can be gauged by MCMC diagnostics. Model choice can be assessed by incorporating another tier in the modelling hierarchy. This requires evaluation of a normalization constant, a notoriously difficult problem. Difficulty with estimating the normalization constant for the MRF can be overcome by using a path integral approach, although this is a highly computationally intensive method. Different methods of estimating ratios of normalization constants (N Cs) are investigated, including importance sampling Monte Carlo (ISMC), dependent Monte Carlo based on MCMC simulations (MCMC), and reverse logistic regression (RLR). I develop an idea present though not fully developed in the literature, and propose the Integrated mean canonical statistic (IMCS) method for estimating log NC ratios for binary MRFs. The IMCS method falls within the framework of the newly identified path sampling methods of Gelman & Meng (1998) and outperforms ISMC, MCMC and RLR. It also does not rely on simplifying assumptions, such as ignoring spatio-temporal dependence in the process. A thorough investigation is made of the application of IMCS to the three-parameter Autologistic model. This work introduces background computations required for the full implementation of the four-tier model in Chapter 7. Two different extensions of the three-tier model to a four-tier version are investigated. The first extension incorporates temporal dependence in the underlying spatio-temporal process. The second extensions allows the successes and failures in the data layer to depend on time. The MCMC computational method is extended to incorporate the extra layer. A major contribution of the thesis is the development of a fully Bayesian approach to inference for these hierarchical models for the first time. Note: The author of this thesis has agreed to make it open access but invites people downloading the thesis to send her an email via the 'Contact Author' function.