Distributed Inference for Spatial Extremes Modeling in High Dimensions.

Extreme environmental events frequently exhibit spatial and temporal dependence. These data are often modeled using max stable processes (MSPs) that are computationally prohibitive to fit for as few as a dozen observations. Supposed computationally-efficient approaches like the composite likelihood remain computationally burdensome with a few hundred observations. In this paper, we propose a spatial partitioning approach based on local modeling of subsets of the spatial domain that delivers computationally and statistically efficient inference. Marginal and dependence parameters of the MSP are estimated locally on subsets of observations using censored pairwise composite likelihood, and combined using a modified generalized method of moments procedure. The proposed distributed approach is extended to estimate inverted MSP models, and to estimate spatially varying coefficient models to deliver computationally efficient modeling of spatial variation in marginal parameters. We demonstrate consistency and asymptotic normality of estimators, and show empirically that our approach leads to statistically efficient estimation of model parameters. We illustrate the flexibility and practicability of our approach through simulations and the analysis of streamflow data from the U.S. Geological Survey.

Associations between Flood Risk and United States Census Tract-Level Health Outcomes.

Human-induced climate change has led to more frequent and severe flooding throughout the globe. We examined the association between flood risk and the prevalence of coronary heart disease, high blood pressure, asthma, and poor mental health in the UnitedStates, while taking into account different levels of social vulnerability. We aggregated flood risk variables from First Street Foundation by census tract and used principal component analysis to derive a set of five interpretable flood risk factors. The dependent variables were census-tract level disease prevalences generated by the Centers for Disease Control and Prevention. Bayesian spatial conditional autoregressive models were fit on this data to quantify the relationship between flood risk and health outcomes under different stratifications of social vulnerability. We showed that three flood risk principal components had small but significant associations with each of the health outcomes, across the different stratifications of social vulnerability. Our analysis gives the first United States-wide estimates of the associated effects of flood risk on specific health outcomes. We also show that social vulnerability is an important moderator of the relationship between flood risk and health outcomes. Our approach can be extended to other ecological studies that examine the health impacts of climate hazards.

A diverse portfolio of marine protected areas can better advance global conservation and equity.

Marine protected areas (MPAs) are widely used for ocean conservation, yet the relative impacts of various types of MPAs are poorly understood. We estimated impacts on fish biomass from no-take and multiple-use (fished) MPAs, employing a rigorous matched counterfactual design with a global dataset of >14,000 surveys in and around 216 MPAs. Both no-take and multiple-use MPAs generated positive conservation outcomes relative to no protection (58.2% and 12.6% fish biomass increases, respectively), with smaller estimated differences between the two MPA types when controlling for additional confounding factors (8.3% increase). Relative performance depended on context and management: no-take MPAs performed better in areas of high human pressure but similar to multiple-use in remote locations. Multiple-use MPA performance was low in high-pressure areas but improved significantly with better management, producing similar outcomes to no-take MPAs when adequately staffed and appropriate use regulations were applied. For priority conservation areas where no-take restrictions are not possible or ethical, our findings show that a portfolio of well-designed and well-managed multiple-use MPAs represents a viable and potentially equitable pathway to advance local and global conservation.

Residential Structural Racism and Prevalence of Chronic Health Conditions.

Importance: Studies elucidating determinants of residential neighborhood-level health inequities are needed. Objective: To quantify associations of structural racism indicators with neighborhood prevalence of chronic kidney disease (CKD), diabetes, and hypertension. Design, Setting, and Participants: This cross-sectional study used public data (2012-2018) and deidentified electronic health records (2017-2018) to describe the burden of structural racism and the prevalence of CKD, diabetes, and hypertension in 150 residential neighborhoods in Durham County, North Carolina, from US census block groups and quantified their associations using bayesian models accounting for spatial correlations and residents' age. Data were analyzed from January 2021 to May 2023. Exposures: Global (neighborhood percentage of White residents, economic-racial segregation, and area deprivation) and discrete (neighborhood child care centers, bus stops, tree cover, reported violent crime, impervious areas, evictions, election participation, income, poverty, education, unemployment, health insurance coverage, and police shootings) indicators of structural racism. Main Outcomes and Measures: Outcomes of interest were neighborhood prevalence of CKD, diabetes, and hypertension. Results: A total of 150 neighborhoods with a median (IQR) of 1708 (1109-2489) residents; median (IQR) of 2% (0%-6%) Asian residents, 30% (16%-56%) Black residents, 10% (4%-20%) Hispanic or Latino residents, 0% (0%-1%) Indigenous residents, and 44% (18%-70%) White residents; and median (IQR) residential income of $54 531 ($37 729.25-$78 895.25) were included in analyses. In models evaluating global indicators, greater burden of structural racism was associated with greater prevalence of CKD, diabetes, and hypertension (eg, per 1-SD decrease in neighborhood White population percentage: CKD prevalence ratio [PR], 1.27; 95% highest density interval [HDI], 1.18-1.35; diabetes PR, 1.43; 95% HDI, 1.37-1.52; hypertension PR, 1.19; 95% HDI, 1.14-1.25). Similarly in models evaluating discrete indicators, greater burden of structural racism was associated with greater neighborhood prevalence of CKD, diabetes, and hypertension (eg, per 1-SD increase in reported violent crime: CKD PR, 1.15; 95% HDI, 1.07-1.23; diabetes PR, 1.20; 95% HDI, 1.13-1.28; hypertension PR, 1.08; 95% HDI, 1.02-1.14). Conclusions and Relevance: This cross-sectional study found several global and discrete structural racism indicators associated with increased prevalence of health conditions in residential neighborhoods. Although inferences from this cross-sectional and ecological study warrant caution, they may help guide the development of future community health interventions.

A cross-sectional analysis of medical conditions and environmental factors associated with fractional exhaled nitric oxide (FeNO) in women and children from the ISA birth cohort, Costa Rica.

BACKGROUND: Fractional exhaled nitric oxide (FeNO) is a marker of airway inflammation. Elevated FeNO has been associated with environmental exposures, however, studies from tropical countries are limited. Using data from the Infants' Environmental Health Study (ISA) birth cohort, we evaluated medical conditions and environmental exposures' association with elevated FeNO. METHODS: We performed a cross-sectional analysis of 277 women and 293 8-year old children who participated in the 8-year post-partum visit in 2019. We measured FeNO and collected information on medical conditions and environmental exposures including smoke from waste burning, work in banana plantations, and home pesticide use. We defined elevated FeNO as >25 ppb for women and >20 ppb for children. To evaluate factors associated with elevated FeNO, we used logistic regression models adjusted for obesity in women and unadjusted in children. RESULTS: Overall elevated FeNO was common (20% of women, 13% of children). Rhinitis diagnosis was significantly associated with elevated FeNO in both women (odds ratio (OR): 3.67 95% Confidence Interval (CI): 1.81,7.35) and children (OR: 8.18 95%CI: 3.15, 21.22); wheeze was associated with elevated FeNO in women (OR: 4.50 95% CI: 2.25, 8.99). Environmental exposures were associated with elevated FeNO, but not significantly. Waste burning was associated with elevated FeNO in both women (OR: 1.58 95%CI 0.68, 4.15) and children (OR: 2.49 95%CI:0.82, 10.79). Para-occupational pesticide exposures were associated with elevated FeNO in women and children. For women, having a partner working in agriculture was associated with elevated FeNO (OR: 1.61 95%CI:0.77, 3.58) and for children, maternal work in agriculture was associated with elevated FeNO. (OR 2.08 95%CI 0.86, 4.67) CONCLUSION: Rhinitis and wheeze were associated with elevated FeNO in this rural, agricultural population. Smoke from waste burning as well as para-occupational pesticide exposure may contribute to elevated FeNO in rural communities.

A nonparametric test of group distributional differences for hierarchically clustered functional data.

Biological sex and gender are critical variables in biomedical research, but are complicated by the presence of sex-specific natural hormone cycles, such as the estrous cycle in female rodents, typically divided into phases. A common feature of these cycles are fluctuating hormone levels that induce sex differences in many behaviors controlled by the electrophysiology of neurons, such as neuronal membrane potential in response to electrical stimulus, typically summarized using a priori defined metrics. In this paper, we propose a method to test for differences in the electrophysiological properties across estrous cycle phase without first defining a metric of interest. We do this by modeling membrane potential data in the frequency domain as realizations of a bivariate process, also depending on the electrical stimulus, by adopting existing methods for longitudinal functional data. We are then able to extract the main features of the bivariate signals through a set of basis function coefficients. We use these coefficients for testing, adapting methods for multivariate data to account for an induced hierarchical structure that is a product of the experimental design. We illustrate the performance of the proposed approach in simulations and then apply the method to experimental data.

Bayesian nonparametric quantile process regression and estimation of marginal quantile effects.

Flexible estimation of multiple conditional quantiles is of interest in numerous applications, such as studying the effect of pregnancy-related factors on low and high birth weight. We propose a Bayesian nonparametric method to simultaneously estimate noncrossing, nonlinear quantile curves. We expand the conditional distribution function of the response in I-spline basis functions where the covariate-dependent coefficients are modeled using neural networks. By leveraging the approximation power of splines and neural networks, our model can approximate any continuous quantile function. Compared to existing models, our model estimates all rather than a finite subset of quantiles, scales well to high dimensions, and accounts for estimation uncertainty. While the model is arbitrarily flexible, interpretable marginal quantile effects are estimated using accumulative local effect plots and variable importance measures. A simulation study shows that our model can better recover quantiles of the response distribution when the data are sparse, and an analysis of birth weight data is presented.

Respiratory and allergic outcomes among 5-year-old children exposed to pesticides.

BACKGROUND: Little is known about the effects of pesticides on children's respiratory and allergic outcomes. We evaluated associations of prenatal and current pesticide exposures with respiratory and allergic outcomes in children from the Infants' Environmental Health Study in Costa Rica. METHODS: Among 5-year-old children (n=303), we measured prenatal and current specific gravity-corrected urinary metabolite concentrations of insecticides (chlorpyrifos, pyrethroids), fungicides (mancozeb, pyrimethanil, thiabendazole) and 2,4-D. We collected information from caregivers on respiratory (ever doctor-diagnosed asthma and lower respiratory tract infections (LRTI), wheeze and cough during last 12 months) and allergic (nasal allergies, itchy rash, ever eczema) outcomes. We fitted separate multivariable logistic regression models for high (≥75th percentile (P75)) vs low (

Spectral adjustment for spatial confounding.

Adjusting for an unmeasured confounder is generally an intractable problem, but in the spatial setting it may be possible under certain conditions. We derive necessary conditions on the coherence between the exposure and the unmeasured confounder that ensure the effect of exposure is estimable. We specify our model and assumptions in the spectral domain to allow for different degrees of confounding at different spatial resolutions. One assumption that ensures identifiability is that confounding present at global scales dissipates at local scales. We show that this assumption in the spectral domain is equivalent to adjusting for global-scale confounding in the spatial domain by adding a spatially smoothed version of the exposure to the mean of the response variable. Within this general framework, we propose a sequence of confounder adjustment methods that range from parametric adjustments based on the Matérn coherence function to more robust semiparametric methods that use smoothing splines. These ideas are applied to areal and geostatistical data for both simulated and real datasets.

Estimating intervention effects on infectious disease control: The effect of community mobility reduction on Coronavirus spread.

Understanding the effects of interventions, such as restrictions on community and large group gatherings, is critical to controlling the spread of COVID-19. Susceptible-Infectious-Recovered (SIR) models are traditionally used to forecast the infection rates but do not provide insights into the causal effects of interventions. We propose a spatiotemporal model that estimates the causal effect of changes in community mobility (intervention) on infection rates. Using an approximation to the SIR model and incorporating spatiotemporal dependence, the proposed model estimates a direct and indirect (spillover) effect of intervention. Under an interference and treatment ignorability assumption, this model is able to estimate causal intervention effects, and additionally allows for spatial interference between locations. Reductions in community mobility were measured by cell phone movement data. The results suggest that the reductions in mobility decrease Coronavirus cases 4 to 7 weeks after the intervention.

A discrete-time survival model for porcine epidemic diarrhoea virus.

Since the arrival of porcine epidemic diarrhea virus (PEDV) in the United States in 2013, elimination and control programmes have had partial success. The dynamics of its spread are hard to quantify, though previous work has shown that local transmission and the transfer of pigs within production systems are most associated with the spread of PEDV. Our work relies on the history of PEDV infections in a region of the southeastern United States. This infection data is complemented by farm-level features and extensive industry data on the movement of both pigs and vehicles. We implement a discrete-time survival model and evaluate different approaches to modelling the local-transmission and network effects. We find strong evidence in that the local-transmission and pig-movement effects are associated with the spread of PEDV, even while controlling for seasonality, farm-level features and the possible spread of disease by vehicles. Our fully Bayesian model permits full uncertainty quantification of these effects. Our farm-level out-of-sample predictions have a receiver-operating characteristic area under the curve (AUC) of 0.779 and a precision-recall AUC of 0.097. The quantification of these effects in a comprehensive model allows stakeholders to make more informed decisions about disease prevention efforts.

Pyrimethanil and chlorpyrifos air concentrations and pregnant women's urinary metabolites in the Infants' Environmental Health Study (ISA), Costa Rica.

BACKGROUND: Only few studies have compared environmental pesticide air concentrations with specific urinary metabolites to evaluate pathways of exposure. Therefore, we compared pyrimethanil and chlorpyrifos concentrations in air with urinary 4-hydroxypyrimethanil (OHP, metabolite of pyrimethanil) and 3,5,6-trichloro-2-pyridinol (TCPy, metabolite of chlorpyrifos) among pregnant women from the Infant's Environmental Health Study (ISA) in Matina County, Costa Rica. METHODS: During pregnancy, we obtained repeat urinary samples from 448 women enrolled in the ISA study. We extrapolated pyrimethanil and chlorpyrifos concentrations measured with passive air samplers (PAS) (n = 48, from 12 schools), across space and time using a Bayesian spatiotemporal model. We subsequently compared these concentrationswith urinary OHP and TCPy in 915 samples from 448 women, usingseparatemixed models andconsidering several covariables. RESULTS: A 10% increase in air pyrimethanil (ng/m3) was associated with a 5.7% (95% confidence interval (CI 4.6, 6.8) increase in OHP (µg/L). Women living further from banana plantations had lower OHP: -0.7% (95% CI -1.2, -0.3) for each 10% increase in distance (meters) as well as women who ate rice and beans ≥15 times a week -23% (95% CI -38, -4). In addition, each 1 ng/m3 increase in chlorpyrifos in air was associated with a 1.5% (95% CI 0.2, 2.8) increase in TCPy (µg/L), and women working in agriculture tended to have increased TCPy (21%, 95% CI -2, 49). CONCLUSION: The Bayesian spatiotemporal models were useful to estimate pyrimethanil and chlorpyrifos air concentrations across space and time. Our results suggest inhalation of pyrimethanil and chlorpyrifos is a pathway of environmental exposure. PAS seems a useful technique to monitor environmental current-use pesticide exposures. For future studies, we recommend increasing the number of locations of environmental air measurements, obtaining all air and urine measurements during the same month, and, ideally, including dermal exposure estimates as well.

Soil Properties and Moisture Synergistically Influence Nontuberculous Mycobacterial Prevalence in Natural Environments of Hawai'i.

Nontuberculous mycobacteria (NTM) are opportunistic pathogens that cause chronic pulmonary disease (PD). NTM infections are thought to be acquired from the environment; however, the basal environmental factors that drive and sustain NTM prevalence are not well understood. The highest prevalence of NTM PD cases in the United States is reported from Hawai'i, which is unique in its climate and soil composition, providing an opportunity to investigate the environmental drivers of NTM prevalence. We used microbiological sampling and spatial logistic regression complemented with fine-scale soil mineralogy to model the probability of NTM presence across the natural landscape of Hawai'i. Over 7 years, we collected and microbiologically cultured 771 samples from 422 geographic sites in natural areas across the Hawaiian Islands for the presence of NTM. NTM were detected in 210 of these samples (27%), with Mycobacterium abscessus being the most frequently isolated species. The probability of NTM presence was highest in expansive soils (those that swell with water) with a high water balance (>1-m difference between rainfall and evapotranspiration) and rich in Fe-oxides/hydroxides. We observed a positive association between NTM presence and iron in wet soils, supporting past studies, but no such association in dry soils. High soil-water balance may facilitate underground movement of NTM into the aquifer system, potentially compounded by expansive capabilities allowing crack formation under drought conditions, representing further possible avenues for aquifer infiltration. These results suggest both precipitation and soil properties are mechanisms by which surface NTM may reach the human water supply. IMPORTANCE Nontuberculous mycobacteria (NTM) are ubiquitous in the environment, being found commonly in soils and natural bodies of freshwater. However, little is known about the environmental niches of NTM and how they relate to NTM prevalence in homes and other human-dominated areas. To characterize NTM environmental associations, we collected and cultured 771 samples from 422 geographic sites in natural areas across Hawai'i, the U.S. state with the highest prevalence of NTM pulmonary disease. We show that the environmental niches of NTM are most associated with highly expansive, moist soils containing high levels of iron oxides/hydroxides. Understanding the factors associated with NTM presence in the natural environment will be crucial for identifying potential mechanisms and risk factors associated with NTM infiltration into water supplies, which are ultimately piped into homes where most exposure risk is thought to occur.

Environmental exposures contribute to respiratory and allergic symptoms among women living in the banana growing regions of Costa Rica.

OBJECTIVES: This research evaluates whether environmental exposures (pesticides and smoke) influence respiratory and allergic outcomes in women living in a tropical, agricultural environment. METHODS: We used data from 266 mothers from the Infants' Environmental Health cohort study in Costa Rica. We evaluated environmental exposures in women by measuring seven pesticide and two polycyclic aromatic hydrocarbons metabolites in urine samples. We defined 'high exposure' as having a metabolite value in the top 75th percentile. We collected survey data on respiratory and allergic outcomes in mothers as well as on pesticides and other environmental exposures. Using logistic regression models adjusted for obesity, we assessed the associations of pesticide exposure with multiple outcomes (wheeze, doctor-diagnosed asthma, high (≥2) asthma score based on symptoms, rhinitis, eczema and itchy rash). RESULTS: Current pesticide use in the home was positively associated with diagnosed asthma (OR=1.99 (95% CI=1.05 to 3.87)). High urinary levels of 5-hydroxythiabendazole (thiabendazole metabolite) and living in a neighbourhood with frequent smoke from waste burning were associated with a high asthma score (OR=1.84 (95%CI=1.05 to 3.25) and OR=2.31 (95%CI=1.11 to 5.16), respectively). Women who worked in agriculture had a significantly lower prevalence of rhinitis (0.19 (0.01 to 0.93)), but were more likely to report eczema (OR=2.54 (95%CI=1.33 to 4.89)) and an itchy rash (OR=3.17 (95%CI=1.24 to 7.73)). CONCLUSIONS: While limited by sample size, these findings suggest that environmental exposure to both pesticides and smoke may impact respiratory and skin-related allergic outcomes in women.

A SPATIAL CAUSAL ANALYSIS OF WILDLAND FIRE-CONTRIBUTED PM2.5 USING NUMERICAL MODEL OUTPUT.

Wildland fire smoke contains hazardous levels of fine particulate matter (PM2.5), a pollutant shown to adversely effect health. Estimating fire attributable PM2.5 concentrations is key to quantifying the impact on air quality and subsequent health burden. This is a challenging problem since only total PM2.5 is measured at monitoring stations and both fire-attributable PM2.5 and PM2.5 from all other sources are correlated in space and time. We propose a framework for estimating fire-contributed PM2.5 and PM2.5 from all other sources using a novel causal inference framework and bias-adjusted chemical model representations of PM2.5 under counterfactual scenarios. The chemical model representation of PM2.5 for this analysis is simulated using Community Multiscale Air Quality Modeling System (CMAQ), run with and without fire emissions across the contiguous U.S. for the 2008-2012 wildfire seasons. The CMAQ output is calibrated with observations from monitoring sites for the same spatial domain and time period. We use a Bayesian model that accounts for spatial variation to estimate the effect of wildland fires on PM2.5 and state assumptions under which the estimate has a valid causal interpretation. Our results include estimates of the contributions of wildfire smoke to PM2.5 for the contiguous U.S. Additionally, we compute the health burden associated with the PM2.5 attributable to wildfire smoke.

Geostatistical modeling of positive-definite matrices: An application to diffusion tensor imaging.

Geostatistical modeling for continuous point-referenced data has extensively been applied to neuroimaging because it produces efficient and valid statistical inference. However, diffusion tensor imaging (DTI), a neuroimaging technique characterizing the brain's anatomical structure, produces a positive-definite (p.d.) matrix for each voxel. Currently, only a few geostatistical models for p.d. matrices have been proposed because introducing spatial dependence among p.d. matrices properly is challenging. In this paper, we use the spatial Wishart process, a spatial stochastic process (random field), where each p.d. matrix-variate random variable marginally follows a Wishart distribution, and spatial dependence between random matrices is induced by latent Gaussian processes. This process is valid on an uncountable collection of spatial locations and is almost-surely continuous, leading to a reasonable way of modeling spatial dependence. Motivated by a DTI data set of cocaine users, we propose a spatial matrix-variate regression model based on the spatial Wishart process. A problematic issue is that the spatial Wishart process has no closed-form density function. Hence, we propose an approximation method to obtain a feasible Cholesky decomposition model, which we show to be asymptotically equivalent to the spatial Wishart process model. A local likelihood approximation method is also applied to achieve fast computation. The simulation studies and real data application demonstrate that the Cholesky decomposition process model produces reliable inference and improved performance, compared to other methods.

A spatiotemporal recommendation engine for malaria control.

Malaria is an infectious disease affecting a large population across the world, and interventions need to be efficiently applied to reduce the burden of malaria. We develop a framework to help policy-makers decide how to allocate limited resources in realtime for malaria control. We formalize a policy for the resource allocation as a sequence of decisions, one per intervention decision, that map up-to-date disease related information to a resource allocation. An optimal policy must control the spread of the disease while being interpretable and viewed as equitable to stakeholders. We construct an interpretable class of resource allocation policies that can accommodate allocation of resources residing in a continuous domain and combine a hierarchical Bayesian spatiotemporal model for disease transmission with a policy-search algorithm to estimate an optimal policy for resource allocation within the pre-specified class. The estimated optimal policy under the proposed framework improves the cumulative long-term outcome compared with naive approaches in both simulation experiments and application to malaria interventions in the Democratic Republic of the Congo.