Project description:In spatial regression models, collinearity between covariates and spatial effects can lead to significant bias in effect estimates. This problem, known as spatial confounding, is encountered modeling forestry data to assess the effect of temperature on tree health. Reliable inference is difficult as results depend on whether or not spatial effects are included in the model. We propose a novel approach, spatial+, for dealing with spatial confounding when the covariate of interest is spatially dependent but not fully determined by spatial location. Using a thin plate spline model formulation we see that, in this case, the bias in covariate effect estimates is a direct result of spatial smoothing. Spatial+ reduces the sensitivity of the estimates to smoothing by replacing the covariates by their residuals after spatial dependence has been regressed away. Through asymptotic analysis we show that spatial+ avoids the bias problems of the spatial model. This is also demonstrated in a simulation study. Spatial+ is straightforward to implement using existing software and, as the response variable is the same as that of the spatial model, standard model selection criteria can be used for comparisons. A major advantage of the method is also that it extends to models with non-Gaussian response distributions. Finally, while our results are derived in a thin plate spline setting, the spatial+ methodology transfers easily to other spatial model formulations.
Project description: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.
Project description:The confounding between fixed effects and (spatial) random effects in a regression setup is termed spatial confounding. This topic continues to gain attention and has been studied extensively in recent years, given that failure to account for this may lead to a suboptimal inference. To mitigate this, a variety of projection-based approaches under the class of restricted spatial models are available in the context of generalized linear mixed models. However, these projection approaches cannot be directly extended to the spatial survival context via frailty models due to dimension incompatibility between the fixed and spatial random effects. In this work, we introduce a two-step approach to handle this, which involves (i) projecting the design matrix to the dimension of the spatial effect (via dimension reduction) and (ii) assuring that the random effect is orthogonal to this new design matrix (confounding alleviation). Under a fully Bayesian paradigm, we conduct fast estimation and inference using integrated nested Laplace approximation. Both simulation studies and application to a motivating data evaluating respiratory cancer survival in the US state of California reveal the advantages of our proposal in terms of model performance and confounding alleviation, compared to alternatives.
Project description:Unmeasured, spatially-structured factors can confound associations between spatial environmental exposures and health outcomes. Adding flexible splines to a regression model is a simple approach for spatial confounding adjustment, but the spline degrees of freedom do not provide an easily interpretable spatial scale. We describe a method for quantifying the extent of spatial confounding adjustment in terms of the Euclidean distance at which variation is removed. We develop this approach for confounding adjustment with splines and using Fourier and wavelet filtering. We demonstrate differences in the spatial scales these bases can represent and provide a comparison of methods for selecting the amount of confounding adjustment. We find the best performance for selecting the amount of adjustment using an information criterion evaluated on an outcome model without exposure. We apply this method to spatial adjustment in an analysis of fine particulate matter and blood pressure in a cohort of United States women.
Project description:We discuss the methods and results of the RESSTE team in the competition on spatial statistics for large datasets. In the first sub-competition, we implemented block approaches both for the estimation of the covariance parameters and for prediction using ordinary kriging. In the second sub-competition, a two-stage procedure was adopted. In the first stage, the marginal distribution is estimated neglecting spatial dependence, either according to the flexible Tuckey g and h distribution or nonparametrically. In the second stage, estimation of the covariance parameters and prediction are performed using Kriging. Vecchias's approximation implemented in the GpGp package proved to be very efficient. We then make some propositions for future competitions.Supplementary informationThe online version supplementary material available at 10.1007/s13253-021-00462-2.
Project description:Propensity score matching is a common tool for adjusting for observed confounding in observational studies, but is known to have limitations in the presence of unmeasured confounding. In many settings, researchers are confronted with spatially-indexed data where the relative locations of the observational units may serve as a useful proxy for unmeasured confounding that varies according to a spatial pattern. We develop a new method, termed distance adjusted propensity score matching (DAPSm) that incorporates information on units' spatial proximity into a propensity score matching procedure. We show that DAPSm can adjust for both observed and some forms of unobserved confounding and evaluate its performance relative to several other reasonable alternatives for incorporating spatial information into propensity score adjustment. The method is motivated by and applied to a comparative effectiveness investigation of power plant emission reduction technologies designed to reduce population exposure to ambient ozone pollution. Ultimately, DAPSm provides a framework for augmenting a "standard" propensity score analysis with information on spatial proximity and provides a transparent and principled way to assess the relative trade-offs of prioritizing observed confounding adjustment versus spatial proximity adjustment.
Project description:Spatial organization of signalling events of the phytohormone auxin is fundamental for maintaining a dynamic transition from plant stem cells to differentiated descendants. The cambium, the stem cell niche mediating wood formation, fundamentally depends on auxin signalling but its exact role and spatial organization is obscure. Here we show that, while auxin signalling levels increase in differentiating cambium descendants, a moderate level of signalling in cambial stem cells is essential for cambium activity. We identify the auxin-dependent transcription factor ARF5/MONOPTEROS to cell-autonomously restrict the number of stem cells by directly attenuating the activity of the stem cell-promoting WOX4 gene. In contrast, ARF3 and ARF4 function as cambium activators in a redundant fashion from outside of WOX4-expressing cells. Our results reveal an influence of auxin signalling on distinct cambium features by specific signalling components and allow the conceptual integration of plant stem cell systems with distinct anatomies.
Project description:BackgroundResidual confounding is challenging to detect. Recently, we described a method for detecting confounding and justified it primarily for time-series studies. The method depends on an indicator with 2 key characteristics: (1) it is conditionally independent (given measured exposures and covariates) of the outcome, in the absence of confounding, misspecification, and measurement errors; and (2) like the exposure, it is associated with confounders, possibly unmeasured. We proposed using future exposure levels as the indicator to detect residual confounding. This choice seems natural for time-series studies because future exposure cannot have caused the event, yet they could be spuriously related to it. A related question addressed here is whether an analogous indicator can be used to identify residual confounding in a study based on spatial, rather than temporal, contrasts.MethodsUsing directed acyclic graphs, we show that future air pollution levels may have the characteristics appropriate for an indicator of residual confounding in spatial studies of environmental exposures. We empirically evaluate performance for spatial studies using simulations.ResultsIn simulations based on a spatial study of ambient air pollution levels and birth weight in Atlanta, and using ambient air pollution 1 year after conception as the indicator, we were able to detect residual confounding. The discriminatory ability approached 100% for some factors intentionally omitted from the model, but was very weak for others.ConclusionThe simulations illustrate that an indicator based on future exposures can have excellent ability to detect residual confounding in spatial studies, although performance varied by situation.
Project description:Wood is a remarkable material with great cultural, economic, and biogeochemical importance. However, our understanding of its formation is poor. Key properties that have not been explained include the anatomy of growth rings (with consistent transitions from low-density earlywood to high density latewood), strong temperature-dependence of latewood density (used for historical temperature reconstructions), the regulation of cell size, and overall growth-temperature relationships in conifer and ring-porous tree species. We have developed a theoretical framework based on observations on Pinus sylvestris L. in northern Sweden. The observed anatomical properties emerge from our framework as a consequence of interactions in time and space between the production of new cells, the dynamics of developmental zone widths, and the distribution of carbohydrates across the developing wood. Here we find that the diffusion of carbohydrates is critical to determining final ring anatomy, potentially overturning current understanding of how wood formation responds to environmental variability and transforming our interpretation of tree rings as proxies of past climates.
Project description:We suggest an information-theoretic approach for measuring stylistic coordination in dialogues. The proposed measure has a simple predictive interpretation and can account for various confounding factors through proper conditioning. We revisit some of the previous studies that reported strong signatures of stylistic accommodation, and find that a significant part of the observed coordination can be attributed to a simple confounding effect--length coordination. Specifically, longer utterances tend to be followed by longer responses, which gives rise to spurious correlations in the other stylistic features. We propose a test to distinguish correlations in length due to contextual factors (topic of conversation, user verbosity, etc.) and turn-by-turn coordination. We also suggest a test to identify whether stylistic coordination persists even after accounting for length coordination and contextual factors.