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.

Limitations of Remotely-sensed Aerosol as a Spatial Proxy for Fine Particulate Matter

Recent research highlights the promise of remotely-sensed aerosol optical depth (AOD) as a proxy for ground-level PM2.5. Particular interest lies in the information on spatial heterogeneity potentially provided by AOD, with important application to estimating and monitoring pollution exposure for public health purposes. Given the temporal and spatio-temporal correlations reported between AOD and PM2.5 , it is tempting to interpret the spatial patterns in AOD as reflecting patterns in PM2.5 . Here we find only limited spatial associations of AOD from three satellite retrievals with PM2.5 over the eastern U.S. at the daily and yearly levels in 2004. We then use statistical modeling to show that the patterns in monthly average AOD poorly reflect patterns in PM2.5 because of systematic, spatially-correlated error in AOD as a proxy for PM2.5 . Furthermore, when we include AOD as a predictor of monthly PM2.5 in a statistical prediction model, AOD provides little additional information to improve predictions of PM2.5 when included in a model that already accounts for land use, emission sources, meteorology and regional variability. These results suggest caution in using spatial variation in AOD to stand in for spatial variation in ground-level PM2.5 in epidemiological analyses and indicate that when PM2.5 monitoring is available, careful statistical modeling outperforms the use of AOD.

A BAYESIAN HIERARCHICAL FRAMEWORK FOR SPATIAL MODELING OF fMRI DATA

Functional neuroimaging techniques enable investigations into the neural basis of human cognition, emotions, and behaviors. In practice, applications of functional magnetic resonance imaging (fMRI) have provided novel insights into the neuropathophysiology of major psychiatric,neurological, and substance abuse disorders, as well as into the neural responses to their treatments. Modern activation studies often compare localized task-induced changes in brain activity between experimental groups. One may also extend voxel-level analyses by simultaneously considering the ensemble of voxels constituting an anatomically defined region of interest (ROI) or by considering means or quantiles of the ROI. In this work we present a Bayesian extension of voxel-level analyses that offers several notable benefits. First, it combines whole-brain voxel-by-voxel modeling and ROI analyses within a unified framework. Secondly, an unstructured variance/covariance for regional mean parameters allows for the study of inter-regional functional connectivity, provided enough subjects are available to allow for accurate estimation. Finally, an exchangeable correlation structure within regions allows for the consideration of intra-regional functional connectivity. We perform estimation for our model using Markov Chain Monte Carlo (MCMC) techniques implemented via Gibbs sampling which, despite the high throughput nature of the data, can be executed quickly (less than 30 minutes). We apply our Bayesian hierarchical model to two novel fMRI data sets: one considering inhibitory control in cocaine-dependent men and the second considering verbal memory in subjects at high risk for Alzheimer’s disease. The unifying hierarchical model presented in this manuscript is shown to enhance the interpretation content of these data sets.

BIVARIATE BINOMIAL SPATIAL MODELLING LOA loa PREVALENCE IN TROPICAL AFRICA

We present a state-of-the-art application of smoothing for dependent bivariate binomial spatial data to Loa loa prevalence mapping in West Africa. This application is special because it starts with the non-spatial calibration of survey instruments, continues with the spatial model building and assessment and ends with robust, tested software that will be used by the field scientists of the World Health Organization for online prevalence map updating. From a statistical perspective several important methodological issues were addressed: (a) building spatial models that are complex enough to capture the structure of the data but remain computationally usable; (b)reducing the computational burden in the handling of very large covariate data sets; (c) devising methods for comparing spatial prediction methods for a given exceedance policy threshold.

Endothelial survival factors and spatial completion, but not pericyte coverage of retinal capillaries determine vessel plasticity

Pericyte loss and capillary regression are characteristic for incipient diabetic retinopathy. Pericyte recruitment is involved in vessel maturation, and ligand-receptor systems contributing to pericyte recruitment are survival factors for endothelial cells in pericyte-free in vitro systems. We studied pericyte recruitment in relation to the susceptibility toward hyperoxia-induced vascular remodeling using the pericyte reporter X-LacZ mouse and the mouse model of retinopathy of prematurity (ROP). Pericytes were found in close proximity to vessels, both during formation of the superficial and the deep capillary layers. When exposure of mice to the ROP was delayed by 24 h, i.e., after the deep retinal layer had formed [at postnatal (p) day 8], preretinal neovascularizations were substantially diminished at p18. Mice with a delayed ROP exposure had 50% reduced avascular zones. Formation of the deep capillary layers at p8 was associated with a combined up-regulation of angiopoietin-1 and PDGF-B, while VEGF was almost unchanged during the transition from a susceptible to a resistant capillary network. Inhibition of Tie-2 function either by soluble Tie-2 or by a sulindac analog, an inhibitor of Tie-2 phosphorylation, resensitized retinal vessels to neovascularizations due to a reduction of the deep capillary network. Inhibition of Tie-2 function had no effect on pericyte recruitment. Our data indicate that the final maturation of the retinal vasculature and its resistance to regressive signals such as hyperoxia depend on the completion of the multilayer structure, in particular the deep capillary layers, and are independent of the coverage by pericytes.

Spatial Assessment of Erosion and its Impact on Soil Fertility in the Tajik Foothills

Efficient planning of soil conservation measures requires, first, to understand the impact of soil erosion on soil fertility with regard to local land cover classes; and second, to identify hot spots of soil erosion and bright spots of soil conservation in a spatially explicit manner. Soil organic carbon (SOC) is an important indicator of soil fertility. The aim of this study was to conduct a spatial assessment of erosion and its impact on SOC for specific land cover classes. Input data consisted of extensive ground truth, a digital elevation model and Landsat 7 imagery from two different seasons. Soil spectral reflectance readings were taken from soil samples in the laboratory and calibrated with results of SOC chemical analysis using regression tree modelling. The resulting model statistics for soil degradation assessments are promising (R2=0.71, RMSEV=0.32). Since the area includes rugged terrain and small agricultural plots, the decision tree models allowed mapping of land cover classes, soil erosion incidence and SOC content classes at an acceptable level of accuracy for preliminary studies. The various datasets were linked in the hot-bright spot matrix, which was developed to combine soil erosion incidence information and SOC content levels (for uniform land cover classes) in a scatter plot. The quarters of the plot show different stages of degradation, from well conserved land to hot spots of soil degradation. The approach helps to gain a better understanding of the impact of soil erosion on soil fertility and to identify hot and bright spots in a spatially explicit manner. The results show distinctly lower SOC content levels on large parts of the test areas, where annual crop cultivation was dominant in the 1990s and where cultivation has now been abandoned. On the other hand, there are strong indications that afforestations and fruit orchards established in the 1980s have been successful in conserving soil resources.

Golden-winged Warbler habitat model validation for northern Wisconsin and central Minnesota

Between 1966 and 2003, the Golden-winged Warbler (Vermivora chrysoptera) experienced declines of 3.4% per year in large parts of the breeding range and has been identified by Partners in Flight as one of 28 land birds requiring expedient action to prevent its continued decline. It is currently being considered for listing under the Endangered Species Act. A major step in advancing our understanding of the status and habitat preferences of Golden-winged Warbler populations in the Upper Midwest was initiated by the publication of new predictive spatially explicit Golden-winged Warbler habitat models for the northern Midwest. Here, I use original data on observed Golden-winged Warbler abundances in Wisconsin and Minnesota to compare two population models: the hierarchical spatial count (HSC) model with the Habitat Suitability Index (HSI) model. I assessed how well the field data compared to the model predictions and found that within Wisconsin, the HSC model performed slightly better than the HSI model whereas both models performed relatively equally in Minnesota. For the HSC model, I found a 10% error of commission in Wisconsin and a 24.2% error of commission for Minnesota. Similarly, the HSI model has a 23% error of commission in Minnesota; in Wisconsin due to limited areas where the HSI model predicted absences, there was incomplete data and I was unable to determine the error of commission for the HSI model. These are sites where the model predicted presences and the Golden-winged Warbler did not occur. To compare predicted abundance from the two models, a 3x3 contingency table was used. I found that when overlapped, the models do not complement one another in identifying Golden-winged Warbler presences. To calculate discrepancy between the models, the error of commission shows that the HSI model has only a 6.8% chance of correctly classifying absences in the HSC model. The HSC model has only 3.3% chance of correctly classifying absences in the HSI model. These findings highlight the importance of grasses for nesting, shrubs used for cover and foraging, and trees for song perches and foraging as key habitat characteristics for breeding territory occupancy by singing males.

Utilizing Gaussian-Schell model beams to mitigate atmospheric turbulence in free space optical communications / by Kyle R. Drexler.

Turbulence affects traditional free space optical communication by causing speckle to appear in the received beam profile. This occurs due to changes in the refractive index of the atmosphere that are caused by fluctuations in temperature and pressure, resulting in an inhomogeneous medium. The Gaussian-Schell model of partial coherence has been suggested as a means of mitigating these atmospheric inhomogeneities on the transmission side. This dissertation analyzed the Gaussian-Schell model of partial coherence by verifying the Gaussian-Schell model in the far-field, investigated the number of independent phase control screens necessary to approach the ideal Gaussian-Schell model, and showed experimentally that the Gaussian-Schell model of partial coherence is achievable in the far-field using a liquid crystal spatial light modulator. A method for optimizing the statistical properties of the Gaussian-Schell model was developed to maximize the coherence of the field while ensuring that it does not exhibit the same statistics as a fully coherent source. Finally a technique to estimate the minimum spatial resolution necessary in a spatial light modulator was developed to effectively propagate the Gaussian-Schell model through a range of atmospheric turbulence strengths. This work showed that regardless of turbulence strength or receiver aperture, transmitting the Gaussian-Schell model of partial coherence instead of a fully coherent source will yield a reduction in the intensity fluctuations of the received field. By measuring the variance of the intensity fluctuations and the received mean, it is shown through the scintillation index that using the Gaussian-Schell model of partial coherence is a simple and straight forward method to mitigate atmospheric turbulence instead of traditional adaptive optics in free space optical communications.

The cranial bone window model: studying angiogenesis of primary and secondary bone tumors by intravital microscopy

The successful treatment of primary and secondary bone tumors in a huge number of cases remains one of the major unsolved challenges in modern medicine. Malignant primary bone tumor growth predominantly occurs in younger people, whereas older people predominantly suffer from secondary bone tumors since up to 85% of the most frequently occurring malignant solid tumors, such as lung, mammary, and prostate carcinomas, metastasize into the bone. It is well known that a tumor's course may be altered by its surrounding tissue. For this reason, reported here is the protocol for the surgical preparation of a cranial bone window in mice as well as the method to implant tumors in this bone window for further investigations of angiogenesis and other microcirculatory parameters in orthotopically growing primary or secondary bone tumors using intravital microscopy. Intravital microscopy represents an internationally accepted and sophisticated experimental method to study angiogenesis, microcirculation, and many other parameters in a wide variety of neoplastic and nonneoplastic tissues. Since most physiologic and pathophysiologic processes are active and dynamic events, one of the major strengths of chronic animal models using intravital microscopy is the possibility of monitoring the regions of interest in vivo continuously up to several weeks with high spatial and temporal resolution. In addition, after the termination of experiments, tissue samples can be excised easily and further examined by various in vitro methods such as histology, immunohistochemistry, and molecular biology.

Visualization of alveolar recruitment in a porcine model of unilateral lung lavage using 3He-MRI

BACKGROUND: In the acute respiratory distress syndrome potentially recruitable lung volume is currently discussed. (3)He-magnetic resonance imaging ((3)He-MRI) offers the possibility to visualize alveolar recruitment directly. METHODS: With the approval of the state animal care committee, unilateral lung damage was induced in seven anesthetized pigs by saline lavage of the right lungs. The left lung served as an intraindividual control (healthy lung). Unilateral lung damage was confirmed by conventional proton MRI and spiral-CT scanning. The total aerated lung volume was determined both at a positive end-expiratory pressure (PEEP) of 0 and 10 mbar from three-dimensionally reconstructed (3)He images, both for healthy and damaged lungs. The fractional increase of aerated volume in damaged and healthy lungs, followed by a PEEP increase from 0 to 10 mbar, was compared. RESULTS: Aerated gas space was visualized with a high spatial resolution in the three-dimensionally reconstructed (3)He-MR images, and aeration defects in the lavaged lung matched the regional distribution of atelectasis in proton MRI. After recruitment and PEEP increase, the aerated volume increased significantly both in healthy lungs from 415 ml [270-445] (median [min-max]) to 481 ml [347-523] and in lavaged lungs from 264 ml [71-424] to 424 ml [129-520]. The fractional increase in lavaged lungs was significantly larger than that in healthy lungs (healthy: 17% [11-38] vs. lavage: 42% [14-90] (P=0.031). CONCLUSION: The (3)He-MRI signal might offer an experimental approach to discriminate atelectatic vs. poor aerated lung areas in a lung damage animal model. Our results confirm the presence of potential recruitable lung volume by either alveolar collapse or alveolar flooding, in accordance with previous reports by computed tomography.

Assessing the Relevance of Anomie Theory for Explaining Spatial Variation in Lethal Criminal Violence: An Aggregate-Level Analysis of Homicide within the United States

One of the most influential statements in the anomie theory tradition has been Merton’s argument that the volume of instrumental property crime should be higher where there is a greater imbalance between the degree of commitment to monetary success goals and the degree of commitment to legitimate means of pursing such goals. Contemporary anomie theories stimulated by Merton’s perspective, most notably Messner and Rosenfeld’s institutional anomie theory, have expanded the scope conditions by emphasizing lethal criminal violence as an outcome to which anomie theory is highly relevant, and virtually all contemporary empirical studies have focused on applying the perspective to explaining spatial variation in homicide rates. In the present paper, we argue that current explications of Merton’s theory and IAT have not adequately conveyed the relevance of the core features of the anomie perspective to lethal violence. We propose an expanded anomie model in which an unbalanced pecuniary value system – the core causal variable in Merton’s theory and IAT – translates into higher levels of homicide primarily in indirect ways by increasing levels of firearm prevalence, drug market activity, and property crime, and by enhancing the degree to which these factors stimulate lethal outcomes. Using aggregate-level data collected during the mid-to-late 1970s for a sample of relatively large social aggregates within the U.S., we find a significant effect on homicide rates of an interaction term reflecting high levels of commitment to monetary success goals and low levels of commitment to legitimate means. Virtually all of this effect is accounted for by higher levels of property crime and drug market activity that occur in areas with an unbalanced pecuniary value system. Our analysis also reveals that property crime is more apt to lead to homicide under conditions of high levels of structural disadvantage. These and other findings underscore the potential value of elaborating the anomie perspective to explicitly account for lethal violence.

Reanalysis-driven climate simulation over CORDEX North America domain using the Canadian Regional Climate Model, version 5: model performance evaluation

The performance of reanalysis-driven Canadian Regional Climate Model, version 5 (CRCM5) in reproducing the present climate over the North American COordinated Regional climate Downscaling EXperiment domain for the 1989–2008 period has been assessed in comparison with several observation-based datasets. The model reproduces satisfactorily the near-surface temperature and precipitation characteristics over most part of North America. Coastal and mountainous zones remain problematic: a cold bias (2–6 °C) prevails over Rocky Mountains in summertime and all year-round over Mexico; winter precipitation in mountainous coastal regions is overestimated. The precipitation patterns related to the North American Monsoon are well reproduced, except on its northern limit. The spatial and temporal structure of the Great Plains Low-Level Jet is well reproduced by the model; however, the night-time precipitation maximum in the jet area is underestimated. The performance of CRCM5 was assessed against earlier CRCM versions and other RCMs. CRCM5 is shown to have been substantially improved compared to CRCM3 and CRCM4 in terms of seasonal mean statistics, and to be comparable to other modern RCMs.

The Role of Online Platforms for Media Markets – Two‐Dimensional Spatial Competition in a Two‐Sided Market

We analyze the market for online and offline media in a model of two-dimensional spatial competition where media outlets sell content and advertising space. Consumer preferences are distributed along the style and type of news coverage where the distance costs may vary across dimensions. For integrated provision of online and offline platforms we show that entering the online market reduces average profits and may even constitute a prisoner's dilemma. Specialized provision may yield polarization in the style and type dimensions. This is in contrast to the maximum–minimum differentiation result previously established in the literature on multidimensional horizontal competition. We show that maximal differentiation in both dimensions occurs due to the discrete nature of the type dimension and asymmetric advertising markets.

Evaluation of control and surveillance strategies for classical swine fever using a simulation model

Classical swine fever (CSF) outbreaks can cause enormous losses in naïve pig populations. How to best minimize the economic damage and number of culled animals caused by CSF is therefore an important research area. The baseline CSF control strategy in the European Union and Switzerland consists of culling all animals in infected herds, movement restrictions for animals, material and people within a given distance to the infected herd and epidemiological tracing of transmission contacts. Additional disease control measures such as pre-emptive culling or vaccination have been recommended based on the results from several simulation models; however, these models were parameterized for areas with high animal densities. The objective of this study was to explore whether pre-emptive culling and emergency vaccination should also be recommended in low- to moderate-density areas such as Switzerland. Additionally, we studied the influence of initial outbreak conditions on outbreak severity to improve the efficiency of disease prevention and surveillance. A spatial, stochastic, individual-animal-based simulation model using all registered Swiss pig premises in 2009 (n=9770) was implemented to quantify these relationships. The model simulates within-herd and between-herd transmission (direct and indirect contacts and local area spread). By varying the four parameters (a) control measures, (b) index herd type (breeding, fattening, weaning or mixed herd), (c) detection delay for secondary cases during an outbreak and (d) contact tracing probability, 112 distinct scenarios were simulated. To assess the impact of scenarios on outbreak severity, daily transmission rates were compared between scenarios. Compared with the baseline strategy (stamping out and movement restrictions) vaccination and pre-emptive culling neither reduced outbreak size nor duration. Outbreaks starting in a herd with weaning piglets or fattening pigs caused higher losses regarding to the number of culled premises and were longer lasting than those starting in the two other index herd types. Similarly, larger transmission rates were estimated for these index herd type outbreaks. A longer detection delay resulted in more culled premises and longer duration and better transmission tracing increased the number of short outbreaks. Based on the simulation results, baseline control strategies seem sufficient to control CSF in low-medium animal-dense areas. Early detection of outbreaks is crucial and risk-based surveillance should be focused on weaning piglet and fattening pig premises.

Present state of global wetland extent and wetland methane modelling: methodology of a model inter-comparison project (WETCHIMP)

The Wetland and Wetland CH4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH4) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO2) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2013). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extent and CH4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH4 fluxes from wetland systems. Even though modelling of wetland extent and CH4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models.