RESUMO
Predicting the response at an unobserved location is a fundamental problem in spatial statistics. Given the difficulty in modeling spatial dependence, especially in non-stationary cases, model-based prediction intervals are at risk of misspecification bias that can negatively affect their validity. Here we present a new approach for model-free nonparametric spatial prediction based on the conformal prediction machinery. Our key observation is that spatial data can be treated as exactly or approximately exchangeable in a wide range of settings. In particular, under an infill asymptotic regime, we prove that the response values are, in a certain sense, locally approximately exchangeable for a broad class of spatial processes, and we develop a local spatial conformal prediction algorithm that yields valid prediction intervals without strong model assumptions like stationarity. Numerical examples with both real and simulated data confirm that the proposed conformal prediction intervals are valid and generally more efficient than existing model-based procedures for large datasets across a range of non-stationary and non-Gaussian settings.
RESUMO
The dynamics that govern disease spread are hard to model because infections are functions of both the underlying pathogen as well as human or animal behavior. This challenge is increased when modeling how diseases spread between different spatial locations. Many proposed spatial epidemiological models require trade-offs to fit, either by abstracting away theoretical spread dynamics, fitting a deterministic model, or by requiring large computational resources for many simulations. We propose an approach that approximates the complex spatial spread dynamics with a Gaussian process. We first propose a flexible spatial extension to the well-known SIR stochastic process, and then we derive a moment-closure approximation to this stochastic process. This moment-closure approximation yields ordinary differential equations for the evolution of the means and covariances of the susceptibles and infectious through time. Because these ODEs are a bottleneck to fitting our model by MCMC, we approximate them using a low-rank emulator. This approximation serves as the basis for our hierarchical model for noisy, underreported counts of new infections by spatial location and time. We demonstrate using our model to conduct inference on simulated infections from the underlying, true spatial SIR jump process. We then apply our method to model counts of new Zika infections in Brazil from late 2015 through early 2016.
Assuntos
Simulação por Computador , Processos Estocásticos , Infecção por Zika virus , Humanos , Distribuição Normal , Infecção por Zika virus/epidemiologia , Infecção por Zika virus/transmissão , Modelos Epidemiológicos , Modelos Estatísticos , Cadeias de MarkovRESUMO
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.
RESUMO
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.
RESUMO
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.
Assuntos
Conservação dos Recursos Naturais , Pesqueiros , Animais , Humanos , Biomassa , Peixes , EcossistemaRESUMO
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.
Assuntos
Diabetes Mellitus , Hipertensão , Insuficiência Renal Crônica , Humanos , Estudos Transversais , Teorema de Bayes , Prevalência , Racismo Sistêmico , Doença Crônica , Hipertensão/epidemiologiaRESUMO
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.
Assuntos
Asma , Praguicidas , Rinite , Lactente , Humanos , Criança , Feminino , Estudos Transversais , Teste da Fração de Óxido Nítrico Exalado , Asma/epidemiologia , Coorte de Nascimento , Costa Rica , Óxido Nítrico/análise , Testes Respiratórios , Fumaça/efeitos adversos , ExpiraçãoRESUMO
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.
Assuntos
Hormônios , Caracteres Sexuais , Animais , Feminino , Masculino , Hormônios/fisiologia , Roedores/fisiologiaRESUMO
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 (Assuntos
Asma
, Exantema
, Praguicidas
, Piretrinas
, Infecções Respiratórias
, Rinite Alérgica
, Lactente
, Gravidez
, Feminino
, Humanos
, Pré-Escolar
, Praguicidas/toxicidade
, Tosse
RESUMO
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.
Assuntos
Redes Neurais de Computação , Feminino , Gravidez , Humanos , Análise de Regressão , Peso ao Nascer , Teorema de Bayes , Simulação por ComputadorRESUMO
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.
RESUMO
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.
RESUMO
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.
Assuntos
Infecções por Coronavirus , Vírus da Diarreia Epidêmica Suína , Doenças dos Suínos , Animais , Estados Unidos/epidemiologia , Suínos , Teorema de Bayes , Infecções por Coronavirus/epidemiologia , Infecções por Coronavirus/prevenção & controle , Infecções por Coronavirus/veterinária , Doenças dos Suínos/epidemiologia , Doenças dos Suínos/prevenção & controle , MovimentoRESUMO
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.
Assuntos
Clorpirifos , Inseticidas , Praguicidas , Humanos , Feminino , Lactente , Gravidez , Clorpirifos/urina , Gestantes , Costa Rica , Teorema de Bayes , Praguicidas/urina , Saúde Ambiental , Inseticidas/urinaRESUMO
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.
Assuntos
Pneumopatias , Infecções por Mycobacterium não Tuberculosas , Havaí/epidemiologia , Humanos , Ferro , Infecções por Mycobacterium não Tuberculosas/epidemiologia , Infecções por Mycobacterium não Tuberculosas/microbiologia , Micobactérias não Tuberculosas , Óxidos , Prevalência , Solo , Estados UnidosRESUMO
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.
RESUMO
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.