Estimation and interpretation of incidence rate difference for recurrent events when the estimation model is misspecified

Recurrent events data are common in experimental and observational studies. It is often of interest to estimate the effect of an intervention on the incidence rate of the recurrent events. The incidence rate difference is a useful measure of intervention effect. A weighted least squares estimator of the incidence rate difference for recurrent events was recently proposed for an additive rate model in which both the baseline incidence rate and the covariate effects were constant over time. In this article, we relax this model assumption and examine the properties of the estimator under the additive and multiplicative rate models assumption in which the baseline incidence rate and covariate effects may vary over time. We show analytically and numerically that the estimator gives an appropriate summary measure of the time‐varying covariate effects. In particular, when the underlying covariate effects are additive and time‐varying, the estimator consistently estimates the weighted average of the covariate effects over time. When the underlying covariate effects are multiplicative and time‐varying, and if there is only one binary covariate indicating the intervention status, the estimator consistently estimates the weighted average of the underlying incidence rate difference between the intervention and control groups over time. We illustrate the method with data from a randomized vaccine trial.     

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.     

Smooth individual level covariates adjustment in disease mapping

Spatial models for disease mapping should ideally account for covariates measured both at individual and area levels. The newly available “indiCAR” model fits the popular conditional autoregresssive (CAR) model by accommodating both individual and group level covariates while adjusting for spatial correlation in the disease rates. This algorithm has been shown to be effective but assumes log‐linear associations between individual level covariates and outcome. In many studies, the relationship between individual level covariates and the outcome may be non‐log‐linear, and methods to track such nonlinearity between individual level covariate and outcome in spatial regression modeling are not well developed. In this paper, we propose a new algorithm, smooth‐indiCAR, to fit an extension to the popular conditional autoregresssive model that can accommodate both linear and nonlinear individual level covariate effects while adjusting for group level covariates and spatial correlation in the disease rates. In this formulation, the effect of a continuous individual level covariate is accommodated via penalized splines. We describe a two‐step estimation procedure to obtain reliable estimates of individual and group level covariate effects where both individual and group level covariate effects are estimated separately. This distributed computing framework enhances its application in the Big Data domain with a large number of individual/group level covariates. We evaluate the performance of smooth‐indiCAR through simulation. Our results indicate that the smooth‐indiCAR method provides reliable estimates of all regression and random effect parameters. We illustrate our proposed methodology with an analysis of data on neutropenia admissions in New South Wales (NSW), Australia.     

Family-specific approaches to the analysis of case-control family data

Case-control studies augmented by the values of responses and covariates from family members allow investigators to study the association between the response and genetics and environment by relating differences in the response directly to within-family differences in covariates. However, existing approaches for case-control family data parameterize covariate effects in terms of the marginal probability of response, the same effects that one estimates from standard case-control studies. This article focuses on the estimation of family-specific covariate effects and develops efficient methods to fit family-specific models such as binary mixed-effects models. We also extend the approach to cover any setting where one has a fully specified model for the vector of responses in a family. We illustrate our approach using data from a case-control family study of brain cancer and consider the use of weighted and conditional likelihood methods as alternatives.     

Estimating causal effects of air quality regulations using principal stratification for spatially correlated multivariate intermediate outcomes

Methods for causal inference regarding health effects of air quality regulations are met with unique challenges because (1) changes in air quality are intermediates on the causal pathway between regulation and health, (2) regulations typically affect multiple pollutants on the causal pathway towards health, and (3) regulating a given location can affect pollution at other locations, that is, there is interference between observations. We propose a principal stratification method designed to examine causal effects of a regulation on health that are and are not associated with causal effects of the regulation on air quality. A novel feature of our approach is the accommodation of a continuously scaled multivariate intermediate response vector representing multiple pollutants. Furthermore, we use a spatial hierarchical model for potential pollution concentrations and ultimately use estimates from this model to assess validity of assumptions regarding interference. We apply our method to estimate causal effects of the 1990 Clean Air Act Amendments among approximately 7 million Medicare enrollees living within 6 miles of a pollution monitor.     

Regression analysis of longitudinal binary data with time-dependent environmental covariates: bias and efficiency

Generalized estimating equations (Liang and Zeger, 1986) is a widely used, moment-based procedure to estimate marginal regression parameters. However, a subtle and often overlooked point is that valid inference requires the mean for the response at time t to be expressed properly as a function of the complete past, present, and future values of any time-varying covariate. For example, with environmental exposures it may be necessary to express the response as a function of multiple lagged values of the covariate series. Despite the fact that multiple lagged covariates may be predictive of outcomes, researchers often focus interest on parameters in a 'cross-sectional' model, where the response is expressed as a function of a single lag in the covariate series. Cross-sectional models yield parameters with simple interpretations and avoid issues of collinearity associated with multiple lagged values of a covariate. Pepe and Anderson (1994), showed that parameter estimates for time-varying covariates may be biased unless the mean, given all past, present, and future covariate values, is equal to the cross-sectional mean or unless independence estimating equations are used. Although working independence avoids potential bias, many authors have shown that a poor choice for the response correlation model can lead to highly inefficient parameter estimates. The purpose of this paper is to study the bias-efficiency trade-off associated with working correlation choices for application with binary response data. We investigate data characteristics or design features (e.g. cluster size, overall response association, functional form of the response association, covariate distribution, and others) that influence the small and large sample characteristics of parameter estimates obtained from several different weighting schemes or equivalently 'working' covariance models. We find that the impact of covariance model choice depends highly on the specific structure of the data features, and that key aspects should be examined before choosing a weighting scheme.     

胡萝卜吸附式低温干燥特性的研究

对胡萝卜吸附式低温干燥过程的干燥特性进行了试验研究。考察了干燥气体 (空气 )的湿度和风量以及物料粒度对胡萝卜干燥特性和复水比的影响 ,得到形状相似的干燥曲线 ,用自定义最小二乘法拟合均能获得较好的结果。结果表明 :吸附式低温干燥过程可使被干燥物料达到超干水平 (含水率 <5 % ) ;增加干燥气体流量对干燥过程进行有利 ,当气体流量从 2 0 0L/h增加到 4 0 0L/h时 ,胡萝卜含水率从 4 1 13%降低到 34 93% ;被干燥物料颗粒大小和形状对干燥过程有显著影响 ;经过吸附式干燥后的胡萝卜色泽鲜艳 ,无褐变 ,外表品质优于热风干燥。吸附式低温干燥后胡萝卜的复水比 (6 5 2 )显著高于热风干燥 (1 78)。     

Statistical Modeling of Agatston Score in Multi-Ethnic Study of Atherosclerosis (MESA)

The MESA (Multi-Ethnic Study of Atherosclerosis) is an ongoing study of the prevalence, risk factors, and progression of subclinical cardiovascular disease in a multi-ethnic cohort. It provides a valuable opportunity to examine the development and progression of CAC (coronary artery calcium), which is an important risk factor for the development of coronary heart disease. In MESA, about half of the CAC scores are zero and the rest are continuously distributed. Such data has been referred to as “zero-inflated data” and may be described using two-part models. Existing two-part model studies have limitations in that they usually consider parametric models only, make the assumption of known forms of the covariate effects, and focus only on the estimation property of the models. In this article, we investigate statistical modeling of CAC in MESA. Building on existing studies, we focus on two-part models. We investigate both parametric and semiparametric, and both proportional and nonproportional models. For various models, we study their estimation as well as prediction properties. We show that, to fully describe the relationship between covariates and CAC development, the semiparametric model with nonproportional covariate effects is needed. In contrast, for the purpose of prediction, the parametric model with proportional covariate effects is sufficient. This study provides a statistical basis for describing the behaviors of CAC and insights into its biological mechanisms.     

Analysis of exponential growth data for yoghurt cultures

Obtaining accurate estimates of maximum specific growth rate, growth yield, and product yield is important for many fermentation processes. A systematic procedure is presented to select the exponential growth region and estimate the maximum specific growth rate using the covariate adjustment method with all the available measured variables (i.e. biomass, substrate, and product). The procedure is applied to data collected during growth of pure and mixed cultures of Lactobacillus bulgaricus and Streptococcus thermophilus on 3% dry milk under anaerobic conditions. The estimation procedure gives good estimates with relatively narrow confidence intervals even though biomass concentration is measured by an indirect method. The estimated values of maximum specific growth rate range from 0.2805 h(-1) for S. thermophilus (ATCC-19258) to 0.4672 h(-1) for S. thermophilus (Microlife). Growth and product yields are estimated using regression analysis and the data for the exponential growth region. The growth yields are compared to their theoretical maximum values.     

A multistate model for bivariate interval-censored failure time data

SUMMARY: Interval-censored life-history data arise when the events of interest are only detectable at periodic assessments. When interest lies in the occurrence of two such events, bivariate-interval censored event time data are obtained. We describe how to fit a four-state Markov model useful for characterizing the association between two interval-censored event times when the assessment times for the two events may be generated by different inspection processes. The approach treats the two events symmetrically and enables one to fit multiplicative intensity models that give estimates of covariate effects as well as relative risks characterizing the association between the two events. An expectation-maximization (EM) algorithm is described for estimation in which the maximization step can be carried out with standard software. The method is illustrated by application to data from a trial of HIV patients where the events are the onset of viral shedding in the blood and urine among individuals infected with cytomegalovirus.     

A Bayesian approach to a logistic regression model with incomplete information

Summary .   We consider a set of independent Bernoulli trials with possibly different success probabilities that depend on covariate values. However, the available data consist only of aggregate numbers of successes among subsets of the trials along with all of the covariate values. We still wish to estimate the parameters of a modeled relationship between the covariates and the success probabilities, e.g., a logistic regression model. In this article, estimation of the parameters is made from a Bayesian perspective by using a Markov chain Monte Carlo algorithm based only on the available data. The proposed methodology is applied to both simulation studies and real data from a dose–response study of a toxic chemical, perchlorate.     

An overview of techniques for linking high-dimensional molecular data to time-to-event endpoints by risk prediction models

Analysis of molecular data promises identification of biomarkers for improving prognostic models, thus potentially enabling better patient management. For identifying such biomarkers, risk prediction models can be employed that link high-dimensional molecular covariate data to a clinical endpoint. In low-dimensional settings, a multitude of statistical techniques already exists for building such models, e.g. allowing for variable selection or for quantifying the added value of a new biomarker. We provide an overview of techniques for regularized estimation that transfer this toward high-dimensional settings, with a focus on models for time-to-event endpoints. Techniques for incorporating specific covariate structure are discussed, as well as techniques for dealing with more complex endpoints. Employing gene expression data from patients with diffuse large B-cell lymphoma, some typical modeling issues from low-dimensional settings are illustrated in a high-dimensional application. First, the performance of classical stepwise regression is compared to stage-wise regression, as implemented by a component-wise likelihood-based boosting approach. A second issues arises, when artificially transforming the response into a binary variable. The effects of the resulting loss of efficiency and potential bias in a high-dimensional setting are illustrated, and a link to competing risks models is provided. Finally, we discuss conditions for adequately quantifying the added value of high-dimensional gene expression measurements, both at the stage of model fitting and when performing evaluation.     

REML estimation of variance parameters in nonlinear mixed effects models using the SAEM algorithm

Nonlinear mixed effects models are now widely used in biometrical studies, especially in pharmacokinetic research or for the analysis of growth traits for agricultural and laboratory species. Most of these studies, however, are often based on ML estimation procedures, which are known to be biased downwards. A few REML extensions have been proposed, but only for approximated methods. The aim of this paper is to present a REML implementation for nonlinear mixed effects models within an exact estimation scheme, based on an integration of the fixed effects and a stochastic estimation procedure. This method was implemented via a stochastic EM, namely the SAEM algorithm. The simulation study showed that the proposed REML estimation procedure considerably reduced the bias observed with the ML estimation, as well as the residual mean squared error of the variance parameter estimations, especially in the unbalanced cases. ML and REML based estimators of fixed effects were also compared via simulation. Although the two kinds of estimates were very close in terms of bias and mean square error, predictions of individual profiles were clearly improved when using REML vs. ML. An application of this estimation procedure is presented for the modelling of growth in lines of chicken.     

Inverse probability weighting methods for Cox regression with right-truncated data

Right-truncated data arise when observations are ascertained retrospectively, and only subjects who experience the event of interest by the time of sampling are selected. Such a selection scheme, without adjustment, leads to biased estimation of covariate effects in the Cox proportional hazards model. The existing methods for fitting the Cox model to right-truncated data, which are based on the maximization of the likelihood or solving estimating equations with respect to both the baseline hazard function and the covariate effects, are numerically challenging. We consider two alternative simple methods based on inverse probability weighting (IPW) estimating equations, which allow consistent estimation of covariate effects under a positivity assumption and avoid estimation of baseline hazards. We discuss problems of identifiability and consistency that arise when positivity does not hold and show that although the partial tests for null effects based on these IPW methods can be used in some settings even in the absence of positivity, they are not valid in general. We propose adjusted estimating equations that incorporate the probability of observation when it is known from external sources, which results in consistent estimation. We compare the methods in simulations and apply them to the analyses of human immunodeficiency virus latency.     

Stochastic Dynamic Model for Estimation of Rate Constants and Their Variances from Noisy and Heterogeneous PET Measurements

Tissue heterogeneity, radioactive decay and measurement noise are the main error sources in compartmental modeling used to estimate the physiologic rate constants of various radiopharmaceuticals from a dynamic PET study. We introduce a new approach to this problem by modeling the tissue heterogeneity with random rate constants in compartment models. In addition, the Poisson nature of the radioactive decay is included as a Poisson random variable in the measurement equations. The estimation problem will be carried out using the maximum likelihood estimation. With this approach, we do not only get accurate mean estimates for the rate constants, but also estimates for tissue heterogeneity within the region of interest and other possibly unknown model parameters, e.g. instrument noise variance, as well. We also avoid the problem of the optimal weighting of the data related to the conventionally used weighted least-squares method. The new approach was tested with simulated time–activity curves from the conventional three compartment – three rate constants model with normally distributed rate constants and with a noise mixture of Poisson and normally distributed random variables. Our simulation results showed that this new model gave accurate estimates for the mean of the rate constants, the measurement noise parameter and also for the tissue heterogeneity, i.e. for the variance of the rate constants within the region of interest.     

Generalized additive distributed lag models: quantifying mortality displacement

There are a number of applied settings where a response is measured repeatedly over time, and the impact of a stimulus at one time is distributed over several subsequent response measures. In the motivating application the stimulus is an air pollutant such as airborne particulate matter and the response is mortality. However, several other variables (e.g. daily temperature) impact the response in a possibly non-linear fashion. To quantify the effect of the stimulus in the presence of covariate data we combine two established regression techniques: generalized additive models and distributed lag models. Generalized additive models extend multiple linear regression by allowing for continuous covariates to be modeled as smooth, but otherwise unspecified, functions. Distributed lag models aim to relate the outcome variable to lagged values of a time-dependent predictor in a parsimonious fashion. The resultant, which we call generalized additive distributed lag models, are seen to effectively quantify the so-called 'mortality displacement effect' in environmental epidemiology, as illustrated through air pollution/mortality data from Milan, Italy.     

Nonparametric tests for continuous covariate effects with multistate survival data

SUMMARY: In clinical trials and observational studies, it is often of scientific interest to evaluate the effects of covariates on complex multistate event probabilities. With discrete covariates, nonparametric tests may be constructed using estimates of the relevant quantities. With continuous covariates, a common approach is to arbitrarily discretize the covariates, which may lead to substantial information loss. Another strategy is to formulate the covariate effects in a regression model. Model-based tests may have either low power or be biased under misspecification. We propose nonparametric tests not requiring arbitrary discretization. The tests involve integrals of estimates continuously indexed by dichotomizations of the covariates. General asymptotic results are derived under null and alternative hypotheses, and verified using empirical process theory in several special cases. The tests are consistent under stochastic ordering, which arises naturally with multistate data. A novel nonparametric measure of covariate effect is studied as a natural byproduct of the testing procedure. Simulation studies and two real data analyses demonstrate the gains of the new testing procedure over those based either on categorization or on regression models.     

Maximum likelihood estimation for cytogenetic dose-response curves

In vitro dose-response curves are used to describe the relation between chromosome aberrations and radiation dose for human lymphocytes. The lymphocytes are exposed to low-LET radiation, and the resulting dicentric chromosome aberrations follow the Poisson distribution. The expected yield depends on both the magnitude and the temporal distribution of the dose. A general dose-response model that describes this relation has been presented by Kellerer and Rossi (1972, Current Topics on Radiation Research Quarterly 8, 85-158; 1978, Radiation Research 75, 471-488) using the theory of dual radiation action. Two special cases of practical interest are split-dose and continuous exposure experiments, and the resulting dose-time-response models are intrinsically nonlinear in the parameters. A general-purpose maximum likelihood estimation procedure is described, and estimation for the nonlinear models is illustrated with numerical examples from both experimental designs. Poisson regression analysis is used for estimation, hypothesis testing, and regression diagnostics. Results are discussed in the context of exposure assessment procedures for both acute and chronic human radiation exposure.     

Allowing for random errors in radiation dose estimates for the atomic bomb survivor data

The presence of random errors in the individual radiation dose estimates for the A-bomb survivors causes underestimation of radiation effects in dose-response analyses, and also distorts the shape of dose-response curves. Statistical methods are presented which will adjust for these biases, provided that a valid statistical model for the dose estimation errors is used. Emphasis is on clarifying some rather subtle statistical issues. For most of this development the distinction between radiation dose and exposure is not critical. The proposed methods involve downward adjustment of dose estimates, but this does not imply that the dosimetry system is faulty. Rather, this is a part of the dose-response analysis required to remove biases in the risk estimates. The primary focus of this report is on linear dose-response models, but methods for linear-quadratic models are also considered briefly. Some plausible models for the dose estimation errors are considered, which have typical errors in a range of 30-40% of the true values, and sensitivity analysis of the resulting bias corrections is provided. It is found that for these error models the resulting estimates of excess cancer risk based on linear models are about 6-17% greater than estimates that make no allowance for dose estimation errors. This increase in risk estimates is reduced to about 4-11% if, as has often been done recently, survivors with dose estimates above 4 Gy are eliminated from the analysis.     

Air pollution combustion emissions: characterization of causative agents and mechanisms associated with cancer, reproductive, and cardiovascular effects

Combustion emissions account for over half of the fine particle (PM(2.5)) air pollution and most of the primary particulate organic matter. Human exposure to combustion emissions including the associated airborne fine particles and mutagenic and carcinogenic constituents (e.g., polycyclic aromatic compounds (PAC), nitro-PAC) have been studied in populations in Europe, America, Asia, and increasingly in third-world counties. Bioassay-directed fractionation studies of particulate organic air pollution have identified mutagenic and carcinogenic polycyclic aromatic hydrocarbons (PAH), nitrated PAH, nitro-lactones, and lower molecular weight compounds from cooking. A number of these components are significant sources of human exposure to mutagenic and carcinogenic chemicals that may also cause oxidative and DNA damage that can lead to reproductive and cardiovascular effects. Chemical and physical tracers have been used to apportion outdoor and indoor and personal exposures to airborne particles between various combustion emissions and other sources. These sources include vehicles (e.g., diesel and gasoline vehicles), heating and power sources (e.g., including coal, oil, and biomass), indoor sources (e.g., cooking, heating, and tobacco smoke), as well as secondary organic aerosols and pollutants derived from long-range transport. Biomarkers of exposure, dose and susceptibility have been measured in populations exposed to air pollution combustion emissions. Biomarkers have included metabolic genotype, DNA adducts, PAH metabolites, and urinary mutagenic activity. A number of studies have shown a significant correlation of exposure to PM(2.5) with these biomarkers. In addition, stratification by genotype increased this correlation. New multivariate receptor models, recently used to determine the sources of ambient particles, are now being explored in the analysis of human exposure and biomarker data. Human studies of both short- and long-term exposures to combustion emissions and ambient fine particulate air pollution have been associated with measures of genetic damage. Long-term epidemiologic studies have reported an increased risk of all causes of mortality, cardiopulmonary mortality, and lung cancer mortality associated with increasing exposures to air pollution. Adverse reproductive effects (e.g., risk for low birth weight) have also recently been reported in Eastern Europe and North America. Although there is substantial evidence that PAH or substituted PAH may be causative agents in cancer and reproductive effects, an increasing number of studies investigating cardiopulmonary and cardiovascular effects are investigating these and other potential causative agents from air pollution combustion sources.