Coherent modeling of longitudinal causal effects on binary outcomes

Analyses of biomedical studies often necessitate modeling longitudinal causal effects. The current focus on personalized medicine and effect heterogeneity makes this task even more challenging. Toward this end, structural nested mean models (SNMMs) are fundamental tools for studying heterogeneous treatment effects in longitudinal studies. However, when outcomes are binary, current methods for estimating multiplicative and additive SNMM parameters suffer from variation dependence between the causal parameters and the noncausal nuisance parameters. This leads to a series of difficulties in interpretation, estimation, and computation. These difficulties have hindered the uptake of SNMMs in biomedical practice, where binary outcomes are very common. We solve the variation dependence problem for the binary multiplicative SNMM via a reparameterization of the noncausal nuisance parameters. Our novel nuisance parameters are variation independent of the causal parameters, and hence allow for coherent modeling of heterogeneous effects from longitudinal studies with binary outcomes. Our parameterization also provides a key building block for flexible doubly robust estimation of the causal parameters. Along the way, we prove that an additive SNMM with binary outcomes does not admit a variation independent parameterization, thereby justifying the restriction to multiplicative SNMMs.     

Combining parametric and nonparametric models to estimate treatment effects in observational studies

Performing causal inference in observational studies requires we assume confounding variables are correctly adjusted for. In settings with few discrete-valued confounders, standard models can be employed. However, as the number of confounders increases these models become less feasible as there are fewer observations available for each unique combination of confounding variables. In this paper, we propose a new model for estimating treatment effects in observational studies that incorporates both parametric and nonparametric outcome models. By conceptually splitting the data, we can combine these models while maintaining a conjugate framework, allowing us to avoid the use of Markov chain Monte Carlo (MCMC) methods. Approximations using the central limit theorem and random sampling allow our method to be scaled to high-dimensional confounders. Through simulation studies we show our method can be competitive with benchmark models while maintaining efficient computation, and illustrate the method on a large epidemiological health survey.     

On-line state recognition in a yeast fed-batch culture using error vectors

The physiological states with respect to cell growth and ethanol production in a yeast fed-batch culture expressed in linguistic form could be recognized on-line by fuzzy inferencing based on error vectors. The error vector was newly defined here in a macroscopic elemental balance equation. The physiological states for cell growth and ethanol production were characterized by error vectors using many experimental data from fed-batch cultures. Fuzzy membership functions were constructed from the frequency distributions of the error vectors and state recognition was performed by fuzzy inferencing. In particular, an unusual physiological state for a yeast cultivation, in which aerobic ethanol production was accompanied by very low cell growth, could be recognized accurately. According to the results of the state recognition, an energy parameter, the P/O ratio in the metabolic reaction model was adaptively estimated, and the cell growth was successfully evaluated with the estimated P/O. (c) 1995 John Wiley & Sons, Inc.     

Assessing the response of forest productivity to climate extremes in Switzerland using model–data fusion

The response of forest productivity to climate extremes strongly depends on ambient environmental and site conditions. To better understand these relationships at a regional scale, we used nearly 800 observation years from 271 permanent long‐term forest monitoring plots across Switzerland, obtained between 1980 and 2017. We assimilated these data into the 3‐PG forest ecosystem model using Bayesian inference, reducing the bias of model predictions from 14% to 5% for forest stem carbon stocks and from 45% to 9% for stem carbon stock changes. We then estimated the productivity of forests dominated by Picea abies and Fagus sylvatica for the period of 1960–2018, and tested for productivity shifts in response to climate along elevational gradient and in extreme years. Simulated net primary productivity (NPP) decreased with elevation (2.86 ± 0.006 Mg C ha^?1 year^?1 km^?1 for P. abies and 0.93 ± 0.010 Mg C ha^?1 year^?1 km^?1 for F. sylvatica). During warm–dry extremes, simulated NPP for both species increased at higher and decreased at lower elevations, with reductions in NPP of more than 25% for up to 21% of the potential species distribution range in Switzerland. Reduced plant water availability had a stronger effect on NPP than temperature during warm‐dry extremes. Importantly, cold–dry extremes had negative impacts on regional forest NPP comparable to warm–dry extremes. Overall, our calibrated model suggests that the response of forest productivity to climate extremes is more complex than simple shift toward higher elevation. Such robust estimates of NPP are key for increasing our understanding of forests ecosystems carbon dynamics under climate extremes.