A varying‐coefficient generalized odds rate model with time‐varying exposure: An application to fitness and cardiovascular disease mortality

Varying‐coefficient models have become a common tool to determine whether and how the association between an exposure and an outcome changes over a continuous measure. These models are complicated when the exposure itself is time‐varying and subjected to measurement error. For example, it is well known that longitudinal physical fitness has an impact on cardiovascular disease (CVD) mortality. It is not known, however, how the effect of longitudinal physical fitness on CVD mortality varies with age. In this paper, we propose a varying‐coefficient generalized odds rate model that allows flexible estimation of age‐modified effects of longitudinal physical fitness on CVD mortality. In our model, the longitudinal physical fitness is measured with error and modeled using a mixed‐effects model, and its associated age‐varying coefficient function is represented by cubic B‐splines. An expectation‐maximization algorithm is developed to estimate the parameters in the joint models of longitudinal physical fitness and CVD mortality. A modified pseudoadaptive Gaussian‐Hermite quadrature method is adopted to compute the integrals with respect to random effects involved in the E‐step. The performance of the proposed method is evaluated through extensive simulation studies and is further illustrated with an application to cohort data from the Aerobic Center Longitudinal Study.     

A cure‐rate model for Q‐learning: Estimating an adaptive immunosuppressant treatment strategy for allogeneic hematopoietic cell transplant patients

Cancers treated by transplantation are often curative, but immunosuppressive drugs are required to prevent and (if needed) to treat graft‐versus‐host disease. Estimation of an optimal adaptive treatment strategy when treatment at either one of two stages of treatment may lead to a cure has not yet been considered. Using a sample of 9563 patients treated for blood and bone cancers by allogeneic hematopoietic cell transplantation drawn from the Center for Blood and Marrow Transplant Research database, we provide a case study of a novel approach to Q‐learning for survival data in the presence of a potentially curative treatment, and demonstrate the results differ substantially from an implementation of Q‐learning that fails to account for the cure‐rate.     

A characterization of missingness at random in a generalized shared‐parameter joint modeling framework for longitudinal and time‐to‐event data,and sensitivity analysis

We consider a conceptual correspondence between the missing data setting, and joint modeling of longitudinal and time‐to‐event outcomes. Based on this, we formulate an extended shared random effects joint model. Based on this, we provide a characterization of missing at random, which is in line with that in the missing data setting. The ideas are illustrated using data from a study on liver cirrhosis, contrasting the new framework with conventional joint models.     

A new mixed‐effects regression model for the analysis of zero‐modified hierarchical count data

Count data sets are traditionally analyzed using the ordinary Poisson distribution. However, such a model has its applicability limited as it can be somewhat restrictive to handle specific data structures. In this case, it arises the need for obtaining alternative models that accommodate, for example, (a) zero‐modification (inflation or deflation at the frequency of zeros), (b) overdispersion, and (c) individual heterogeneity arising from clustering or repeated (correlated) measurements made on the same subject. Cases (a)–(b) and (b)–(c) are often treated together in the statistical literature with several practical applications, but models supporting all at once are less common. Hence, this paper's primary goal was to jointly address these issues by deriving a mixed‐effects regression model based on the hurdle version of the Poisson–Lindley distribution. In this framework, the zero‐modification is incorporated by assuming that a binary probability model determines which outcomes are zero‐valued, and a zero‐truncated process is responsible for generating positive observations. Approximate posterior inferences for the model parameters were obtained from a fully Bayesian approach based on the Adaptive Metropolis algorithm. Intensive Monte Carlo simulation studies were performed to assess the empirical properties of the Bayesian estimators. The proposed model was considered for the analysis of a real data set, and its competitiveness regarding some well‐established mixed‐effects models for count data was evaluated. A sensitivity analysis to detect observations that may impact parameter estimates was performed based on standard divergence measures. The Bayesian ‐value and the randomized quantile residuals were considered for model diagnostics.     

Semiparametric transformation models for interval‐censored data in the presence of a cure fraction

Mixed case interval‐censored data arise when the event of interest is known only to occur within an interval induced by a sequence of random examination times. Such data are commonly encountered in disease research with longitudinal follow‐up. Furthermore, the medical treatment has progressed over the last decade with an increasing proportion of patients being cured for many types of diseases. Thus, interest has grown in cure models for survival data which hypothesize a certain proportion of subjects in the population are not expected to experience the events of interest. In this article, we consider a two‐component mixture cure model for regression analysis of mixed case interval‐censored data. The first component is a logistic regression model that describes the cure rate, and the second component is a semiparametric transformation model that describes the distribution of event time for the uncured subjects. We propose semiparametric maximum likelihood estimation for the considered model. We develop an EM type algorithm for obtaining the semiparametric maximum likelihood estimators (SPMLE) of regression parameters and establish their consistency, efficiency, and asymptotic normality. Extensive simulation studies indicate that the SPMLE performs satisfactorily in a wide variety of settings. The proposed method is illustrated by the analysis of the hypobaric decompression sickness data from National Aeronautics and Space Administration.     

Zero‐modified Poisson model: Bayesian approach,influence diagnostics,and an application to a Brazilian leptospirosis notification data

In this paper, a Bayesian method for inference is developed for the zero‐modified Poisson (ZMP) regression model. This model is very flexible for analyzing count data without requiring any information about inflation or deflation of zeros in the sample. A general class of prior densities based on an information matrix is considered for the model parameters. A sensitivity study to detect influential cases that can change the results is performed based on the Kullback–Leibler divergence. Simulation studies are presented in order to illustrate the performance of the developed methodology. Two real datasets on leptospirosis notification in Bahia State (Brazil) are analyzed using the proposed methodology for the ZMP model.     

Using a second‐order differential model to fit data without baselines in protein isothermal chemical denaturation

In vitro protein stability studies are commonly conducted via thermal or chemical denaturation/renaturation of protein. Conventional data analyses on the protein unfolding/(re)folding require well‐defined pre‐ and post‐transition baselines to evaluate Gibbs free‐energy change associated with the protein unfolding/(re)folding. This evaluation becomes problematic when there is insufficient data for determining the pre‐ or post‐transition baselines. In this study, fitting on such partial data obtained in protein chemical denaturation is established by introducing second‐order differential (SOD) analysis to overcome the limitations that the conventional fitting method has. By reducing numbers of the baseline‐related fitting parameters, the SOD analysis can successfully fit incomplete chemical denaturation data sets with high agreement to the conventional evaluation on the equivalent completed data, where the conventional fitting fails in analyzing them. This SOD fitting for the abbreviated isothermal chemical denaturation further fulfills data analysis methods on the insufficient data sets conducted in the two prevalent protein stability studies.     

Joint structural and physiological control on the interannual variation in productivity in a temperate grassland: A data‐model comparison

Given the important contributions of semiarid region to global land carbon cycle, accurate modeling of the interannual variability (IAV) of terrestrial gross primary productivity (GPP) is important but remains challenging. By decomposing GPP into leaf area index (LAI) and photosynthesis per leaf area (i.e., GPP_leaf), we investigated the IAV of GPP and the mechanisms responsible in a temperate grassland of northwestern China. We further assessed six ecosystem models for their capabilities in reproducing the observed IAV of GPP in a temperate grassland from 2004 to 2011 in China. We observed that the responses to LAI and GPP_leaf to soil water significantly contributed to IAV of GPP at the grassland ecosystem. Two of six models with prescribed LAI simulated of the observed IAV of GPP quite well, but still underestimated the variance of GPP_leaf, therefore the variance of GPP. In comparison, simulated pattern by the other four models with prognostic LAI differed significantly from the observed IAV of GPP. Only some models with prognostic LAI can capture the observed sharp decline of GPP in drought years. Further analysis indicated that accurately representing the responses of GPP_leaf and leaf stomatal conductance to soil moisture are critical for the models to reproduce the observed IAV of GPP_leaf. Our framework also identified that the contributions of LAI and GPP_leaf to the observed IAV of GPP were relatively independent. We conclude that our framework of decomposing GPP into LAI and GPP_leaf has a significant potential for facilitating future model intercomparison, benchmarking and optimization should be adopted for future data‐model comparisons.     

Cognitive flexibility deficits in a mouse model for the absence of full‐length dystrophin

Duchenne muscular dystrophy (DMD) is a progressive muscle‐wasting disorder, caused by mutations in the DMD gene and the resulting lack of dystrophin. The DMD gene has seven promoters, giving rise to multiple full‐length and shorter isoforms. Besides the expression of dystrophin in muscles, the majority of dystrophin isoforms is expressed in brain and dystrophinopathy can lead to cognitive deficits, including intellectual impairments and deficits in executive function. In contrast to the muscle pathology, the impact of the lack of dystrophin on the brain is not very well studied. Here, we study the behavioral consequences of a lack of full‐length dystrophin isoforms in mdx mice, particularly with regard to domains of executive functions and anxiety. We observed a deficit in cognitive flexibility in mdx mice in the absence of motor dysfunction or general learning impairments using two independent behavioral tests. In addition, increased anxiety was observed, but its expression depended on the context. Overall, these results suggest that the absence of full‐length dystrophin in mice has specific behavioral effects that compare well to deficits observed in DMD patients.     

Scale and uncertainty in modeled soil organic carbon stock changes for US croplands using a process‐based model

Process‐based model analyses are often used to estimate changes in soil organic carbon (SOC), particularly at regional to continental scales. However, uncertainties are rarely evaluated, and so it is difficult to determine how much confidence can be placed in the results. Our objective was to quantify uncertainties across multiple scales in a process‐based model analysis, and provide 95% confidence intervals for the estimates. Specifically, we used the Century ecosystem model to estimate changes in SOC stocks for US croplands during the 1990s, addressing uncertainties in model inputs, structure and scaling of results from point locations to regions and the entire country. Overall, SOC stocks increased in US croplands by 14.6 Tg C yr^?1 from 1990 to 1995 and 17.5 Tg C yr^?1 during 1995 to 2000, and uncertainties were ±22% and ±16% for the two time periods, respectively. Uncertainties were inversely related to spatial scale, with median uncertainties at the regional scale estimated at ±118% and ±114% during the early and latter part of 1990s, and even higher at the site scale with estimates at ±739% and ±674% for the time periods, respectively. This relationship appeared to be driven by the amount of the SOC stock change; changes in stocks that exceeded 200 Gg C yr^?1 represented a threshold where uncertainties were always lower than ±100%. Consequently, the amount of uncertainty in estimates derived from process‐based models will partly depend on the level of SOC accumulation or loss. In general, the majority of uncertainty was associated with model structure in this application, and so attaining higher levels of precision in the estimates will largely depend on improving the model algorithms and parameterization, as well as increasing the number of measurement sites used to evaluate the structural uncertainty.