Estimating the effects of environmental exposures using a weighted mean of monitoring stations

The health effects of environmental hazards are often examined using time series of the association between a daily response variable (e.g., death) and a daily level of exposure (e.g., temperature). Exposures are usually the average from a network of stations. This gives each station equal importance, and negates the opportunity for some stations to be better measures of exposure. We used a Bayesian hierarchical model that weighted stations using random variables between zero and one. We compared the weighted estimates to the standard model using data on health outcomes (deaths and hospital admissions) and exposures (air pollution and temperature) in Brisbane, Australia. The improvements in model fit were relatively small, and the estimated health effects of pollution were similar using either the standard or weighted estimates. Spatial weighted exposures would be probably more worthwhile when there is either greater spatial detail in the health outcome, or a greater spatial variation in exposure.

Marginal reversible jump Markov chain Monte Carlo with application to motor unit number estimation

Motor unit number estimation (MUNE) is a method which aims to provide a quantitative indicator of progression of diseases that lead to loss of motor units, such as motor neurone disease. However the development of a reliable, repeatable and fast real-time MUNE method has proved elusive hitherto. Ridall et al. (2007) implement a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to produce a posterior distribution for the number of motor units using a Bayesian hierarchical model that takes into account biological information about motor unit activation. However we find that the approach can be unreliable for some datasets since it can suffer from poor cross-dimensional mixing. Here we focus on improved inference by marginalising over latent variables to create the likelihood. In particular we explore how this can improve the RJMCMC mixing and investigate alternative approaches that utilise the likelihood (e.g. DIC (Spiegelhalter et al., 2002)). For this model the marginalisation is over latent variables which, for a larger number of motor units, is an intractable summation over all combinations of a set of latent binary variables whose joint sample space increases exponentially with the number of motor units. We provide a tractable and accurate approximation for this quantity and also investigate simulation approaches incorporated into RJMCMC using results of Andrieu and Roberts (2009).

Endotoxins in indoor air and settled dust in primary schools in a subtropical climate

Endotoxins can significantly affect the air quality in school environments. However, there is currently no reliable method for the measurement of endotoxins and there is a lack of reference values for endotoxin concentrations to aid in the interpretation of measurement results in school settings. We benchmarked the “baseline” range of endotoxin concentration in indoor air, together with endotoxin load in floor dust, and evaluated the correlation between endotoxin levels in indoor air and settled dust, as well as the effects of temperature and humidity on these levels in subtropical school settings. Bayesian hierarchical modeling indicated that the concentration in indoor air and the load in floor dust were generally (<95th percentile) < 13 EU/m3 and < 24,570 EU/m2, respectively. Exceeding these levels would indicate abnormal sources of endotoxins in the school environment, and the need for further investigation. Metaregression indicated no relationship between endotoxin concentration and load, which points to the necessity for measuring endotoxin levels in both the air and settled dust. Temperature increases were associated with lower concentrations in indoor air and higher loads in floor dust. Higher levels of humidity may be associated with lower airborne endotoxin concentrations.

School children’s personal exposure to ultrafine particles in the urban environment

There has been considerable scientific interest in personal exposure to ultrafine particles (UFP). In this study, the inhaled particle surface area doses and dose relative intensities in the tracheobronchial and alveolar regions of lungs were calculated using the measured 24-hour UFP time series of school children personal exposures for each recorded activity. Bayesian hierarchical modelling was used to determine mean doses and dose intensities for the various microenvironments. Analysis of measured personal exposures for 137 participating children from 25 schools in the Brisbane Metropolitan Area showed similar trends for all the participating children. Bayesian regression modelling was performed to calculate the daily proportion of children's total doses at different microenvironments. The proportion of alveolar doses in the total daily dose for \emph{home}, \emph{school}, \emph{commuting} and \emph{other} were 55.3\%, 35.3\%, 4.5\% and 5.0\%, respectively, with the \emph{home} microenvironment contributing a majority of children's total daily dose. Children's mean indoor dose was never higher than the outdoor's at any of the schools, indicating there were no persistent indoor particle sources in the classrooms during the measurements. Outdoor activities, eating/cooking at home and commuting were the three activities with the highest dose intensities. Personal exposure was more influenced by the ambient particle levels than immediate traffic.

An additive model for spatio-temporal smoothing of cancer mortality rates

In this paper, a Bayesian hierarchical model is used to anaylze the female breast cancer mortality rates for the State of Missouri from 1969 through 2001. The logit transformations of the mortality rates are assumed to be linear over the time with additive spatial and age effects as intercepts and slopes. Objective priors of the hierarchical model are explored. The Bayesian estimates are quite robustness in terms change of the hyperparamaters. The spatial correlations are appeared in both intercepts and slopes.

Geographical distribution of human Schistosoma japonicum infection in the Philippines : tools to support disease control and further elimination

S. japonicum infection is believed to be endemic in 28 of the 80 provinces of the Philippines and the most recent data on schistosomiasis prevalence have shown considerable variability between provinces. In order to increase the efficient allocation of parasitic disease control resources in the country, we aimed to describe the small scale spatial variation in S. japonicum prevalence across the Philippines, quantify the role of the physical environment in driving the spatial variation of S. japonicum, and develop a predictive risk map of S. japonicum infection. Data on S. japonicum infection from 35,754 individuals across the country were geo-located at the barangay level and included in the analysis. The analysis was then stratified geographically for Luzon, the Visayas and Mindanao. Zero-inflated binomial Bayesian geostatistical models of S. japonicum prevalence were developed and diagnostic uncertainty was incorporated. Results of the analysis show that in the three regions, males and individuals aged ≥ 20 years had significantly higher prevalence of S. japonicum compared with females and children <5 years. The role of the environmental variables differed between regions of the Philippines. S. japonicum infection was widespread in the Visayas whereas it was much more focal in Luzon and Mindanao. This analysis revealed significant spatial variation in prevalence of S. japonicum infection in the Philippines. This suggests that a spatially targeted approach to schistosomiasis interventions, including mass drug administration, is warranted. When financially possible, additional schistosomiasis surveys should be prioritized to areas identified to be at high risk, but which were underrepresented in our dataset.

Airborne culturable fungi in naturally ventilated primary school environments in a subtropical climate

There is currently a lack of reference values for indoor air fungal concentrations to allow for the interpretation of measurement results in subtropical school settings. Analysis of the results of this work established that, in the majority of properly maintained subtropical school buildings, without any major affecting events such as floods or visible mould or moisture contamination, indoor culturable fungi levels were driven by outdoor concentration. The results also allowed us to benchmark the “baseline range” concentrations for total culturable fungi, Penicillium spp., Cladosporium spp. and Aspergillus spp. in such school settings. The measured concentration of total culturable fungi and three individual fungal genera were estimated using Bayesian hierarchical modelling. Pooling of these estimates provided a predictive distribution for concentrations at an unobserved school. The results indicated that “baseline” indoor concentration levels for indoor total fungi, Penicillium spp., Cladosporium spp. and Aspergillus spp. in such school settings were generally ≤ 1450, ≤ 680, ≤ 480 and ≤ 90 cfu/m3, respectively, and elevated levels would indicate mould damage in building structures. The indoor/outdoor ratio for most classrooms had 95% credible intervals containing 1, indicating that fungi concentrations are generally the same indoors and outdoors at each school. Bayesian fixed effects regression modeling showed that increasing both temperature and humidity resulted in higher levels of fungi concentration.

Predictors of mental health in midlife and older Australian women: A multilevel investigation

The purpose of this article is to examine the factors associated with women's mental health. A random sample of 340 Australian women aged 40–55 completed surveys on menopausal and lifestyle factors and mental health at three time points. We used hierarchical models to show that decrements in mental health were associated with a corresponding increase in some midlife symptoms (p < .01), time (p < .01), and poor physical health (p < .01), but the effect was not permanent. In older women, mental health was associated with physical functioning, climacteric symptoms, and time, while individual variations in mental health score were largely explained by lifestyle factors.

Field evaluation of tolerance to Tobacco streak virus in sunflower germplasm, and observations of seasonal disease spread

Strong statistical evidence was found for differences in tolerance to natural infections of Tobacco streak virus (TSV) in sunflower hybrids. Data from 470 plots involving 23 different sunflower hybrids tested in multiple trials over 5 years in Australia were analysed. Using a Bayesian Hierarchical Logistic Regression model for analysis provided: (i) a rigorous method for investigating the relative effects of hybrid, seasonal rainfall and proximity to inoculum source on the incidence of severe TSV disease; (ii) a natural method for estimating the probability distributions of disease incidence in different hybrids under historical rainfall conditions; and (iii) a method for undertaking all pairwise comparisons of disease incidence between hybrids whilst controlling the familywise error rate without any drastic reduction in statistical power. The tolerance identified in field trials was effective against the main TSV strain associated with disease outbreaks, TSV-parthenium. Glasshouse tests indicate this tolerance to also be effective against the other TSV strain found in central Queensland, TSV-crownbeard. The use of tolerant germplasm is critical to minimise the risk of TSV epidemics in sunflower in this region. We found strong statistical evidence that rainfall during the early growing months of March and April had a negative effect on the incidence of severe infection with greatly reduced disease incidence in years that had high rainfall during this period.

Mapping the risk of soil-transmitted helminthic infections in the Philippines

Background In order to increase the efficient allocation of soil-transmitted helminth (STH) disease control resources in the Philippines, we aimed to describe for the first time the spatial variation in the prevalence of A. lumbricoides, T. trichiura and hookworm across the country, quantify the association between the physical environment and spatial variation of STH infection and develop predictive risk maps for each infection. Methodology/Principal Findings Data on STH infection from 35,573 individuals across the country were geolocated at the barangay level and included in the analysis. The analysis was stratified geographically in two major regions: 1) Luzon and the Visayas and 2) Mindanao. Bayesian geostatistical models of STH prevalence were developed, including age and sex of individuals and environmental variables (rainfall, land surface temperature and distance to inland water bodies) as predictors, and diagnostic uncertainty was incorporated. The role of environmental variables was different between regions of the Philippines. This analysis revealed that while A. lumbricoides and T. trichiura infections were widespread and highly endemic, hookworm infections were more circumscribed to smaller foci in the Visayas and Mindanao. Conclusions/Significance This analysis revealed significant spatial variation in STH infection prevalence within provinces of the Philippines. This suggests that a spatially targeted approach to STH interventions, including mass drug administration, is warranted. When financially possible, additional STH surveys should be prioritized to high-risk areas identified by our study in Luzon.

空间关系类别的形成与计算

Spatial categorization and computation were studied systematically and sequentially in the dissertation, in which, category learning paradigm was transferred to spatial domain. Four parts were developed in our study, in which, the former two parts investigated the mechanism of spatial categorization, and the latter two parts explored two applied issues. Spatial category formation was investigated by two aspects in distance and angle respectively in research one. The results showed that subjects can make accurate judgment and rating for corresponding spatial categorical prototypes and boundaries. In the second research, the formation of spatial relationship categorization was examined through two factors of number of categories and different spatial configurations. The results were repeated in this part, and also, different influences of reference of frame were found in the process of spatial categorization. In the third research, the influences of spatial categorization on spatial localization were found, that in subject could make different location judgment with new acquired spatial categories. Spatial position distribution has no effect on spatial localization on the other hand. The result that reference objects have significant influences on applicability of spatial relationship was found in research four. Our Bayesian hierarchical model could be applied in simulation and prediction of perception of “above” spatial relationship. Theoretical significance and applied importance were discussed in detail.

Understanding GPU programming for statistical computation: Studies in massively parallel massive mixtures

This article describes advances in statistical computation for large-scale data analysis in structured Bayesian mixture models via graphics processing unit (GPU) programming. The developments are partly motivated by computational challenges arising in fitting models of increasing heterogeneity to increasingly large datasets. An example context concerns common biological studies using high-throughput technologies generating many, very large datasets and requiring increasingly high-dimensional mixture models with large numbers of mixture components.We outline important strategies and processes for GPU computation in Bayesian simulation and optimization approaches, give examples of the benefits of GPU implementations in terms of processing speed and scale-up in ability to analyze large datasets, and provide a detailed, tutorial-style exposition that will benefit readers interested in developing GPU-based approaches in other statistical models. Novel, GPU-oriented approaches to modifying existing algorithms software design can lead to vast speed-up and, critically, enable statistical analyses that presently will not be performed due to compute time limitations in traditional computational environments. Supplementalmaterials are provided with all source code, example data, and details that will enable readers to implement and explore the GPU approach in this mixture modeling context. © 2010 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America.

La visión artificial, un nuevo aliado para el análisis de imágenes artísticas

En este estudio se evalúa el rendimiento de los métodos de Bag-of-Visualterms (BOV) para la clasificación automática de imágenes digitales de la base de datos del artista Miquel Planas. Estas imágenes intervienen en la ideación y diseño de su producción escultórica. Constituye un interesante desafío dada la dificultad de la categorización de escenas cuando éstas difieren más por los contenidos semánticos que por los objetos que contienen. Hemos empleado un método de reconocimiento basado en Kernels introducido por Lazebnik, Schmid y Ponce en 2006. Los resultados son prometedores, en promedio, la puntuación del rendimiento es aproximadamente del 70%. Los experimentos sugieren que la categorización automática de imágenes basada en métodos de visión artificial puede proporcionar principios objetivos en la catalogación de imágenes y que los resultados obtenidos pueden ser aplicados en diferentes campos de la creación artística.

Do non-native invasive fish support elevated lamprey populations?

Summary
1: Managing populations of predators and their prey to achieve conservation or resource management goals is usually technically challenging and frequently socially controversial. This is true even in the simplest ecosystems but can be made much worse when predator–prey relationships are in?uenced by complex interactions, such as biological invasions, population trends or animal movements.
2: Lough Neagh in Northern Ireland is a European stronghold for pollan Coregonus autumnalis, a coregonine ?sh and for river lampreyLampetra ?uviatilis, which feeds parasitically as an adult. Both species are of high conservation importance. Lampreys are known to consume pollan but detailed knowledge of their interactions is scant. While pollan is well known to be a landlocked species in Ireland, the life cycle of normally anadromous river lamprey in Lough Neagh has been unclear. The Lough is also a highly perturbed ecosystem, supporting several invasive, non-native ?sh species that have the potential to in?uence lamprey–pollan interactions.
3: We applied stable isotope techniques to resolve both the movement patterns of lamprey and trophic interactions in this complex community. Recognizing that stable isotope studies are often hampered by high-levels of variability and uncertainty in the systems of interest, we employed novel Bayesian mixing models, which incorporate variability and uncertainty.
4: Stable isotope analyses identi?ed troutSalmo trutta and non-native breamAbramis brama as the main items in lamprey diet. Pollan only represented a major food source for lamprey between May and July.
5: Stable isotope ratios of carbon in tissues from 71 adult lamprey showed no evidence of marine carbon sources, strongly suggesting that Lough Neagh is host to a highly unusual, nonanadromous freshwater population. This ?nding marks out the Lough’s lamprey population as of particular scienti?c interest and enhances the conservation signi?cance of this feature of the Lough.
6: Synthesis and applications.Our Bayesian isotopic mixing models illustrate an unusual pattern of animal movement, enhancing conservation interest in an already threatened population. We have also revealed a complex relationship between lamprey and their food species that is suggestive of hyperpredation, whereby non-native species may sustain high lamprey populations that may in turn be detrimental to native pollan.Long-term conservation of lamprey and pollan in this system is likely to require management intervention, but in light of this exceptional complexity, no simple management options are currently supported. Conservation plans will require better characterization ofpopulation-level interactions and simulation modelling of interventions. More generally, our study demonstrates the importance of considering a full range of possible trophic interactions, particularly in complex ecosystems, and highlights Bayesian isotopic mixing models as powerful tools in resolving trophic relationships.
Key-words: Bayesian, conservation dilemma, Coregonus autumnalis, hyperpredation, Lampetra ?uviatilis, pollan, potamodromous, River lamprey, stable isotope analysis in R, stable isotope

Incorporating uncertainty in predictive species distribution modelling

Motivated by the need to solve ecological problems (climate change, habitat fragmentation and biological invasions), there has been increasing interest in species distribution models (SDMs). Predictions from these models inform conservation policy, invasive species management and disease-control measures. However, predictions are subject to uncertainty, the degree and source of which is often unrecognized. Here, we review the SDM literature in the context of uncertainty, focusing on three main classes of SDM: niche-based models, demographic models and process-based models. We identify sources of uncertainty for each class and discuss how uncertainty can be minimized or included in the modelling process to give realistic measures of confidence around predictions. Because this has typically not been performed, we conclude that uncertainty in SDMs has often been underestimated and a false precision assigned to predictions of geographical distribution. We identify areas where development of new statistical tools will improve predictions from distribution models, notably the development of hierarchical models that link different types of distribution model and their attendant uncertainties across spatial scales. Finally, we discuss the need to develop more defensible methods for assessing predictive performance, quantifying model goodness-of-fit and for assessing the significance of model covariates.