Double Inverse‐Weighted Estimation of Cumulative Treatment Effects Under Nonproportional Hazards and Dependent Censoring

Summary In medical studies of time‐to‐event data, nonproportional hazards and dependent censoring are very common issues when estimating the treatment effect. A traditional method for dealing with time‐dependent treatment effects is to model the time‐dependence parametrically. Limitations of this approach include the difficulty to verify the correctness of the specified functional form and the fact that, in the presence of a treatment effect that varies over time, investigators are usually interested in the cumulative as opposed to instantaneous treatment effect. In many applications, censoring time is not independent of event time. Therefore, we propose methods for estimating the cumulative treatment effect in the presence of nonproportional hazards and dependent censoring. Three measures are proposed, including the ratio of cumulative hazards, relative risk, and difference in restricted mean lifetime. For each measure, we propose a double inverse‐weighted estimator, constructed by first using inverse probability of treatment weighting (IPTW) to balance the treatment‐specific covariate distributions, then using inverse probability of censoring weighting (IPCW) to overcome the dependent censoring. The proposed estimators are shown to be consistent and asymptotically normal. We study their finite‐sample properties through simulation. The proposed methods are used to compare kidney wait‐list mortality by race.     

Regression analysis for recurrent events data under dependent censoring

Summary Recurrent events data are commonly seen in longitudinal follow‐up studies. Dependent censoring often occurs due to death or exclusion from the study related to the disease process. In this article, we assume flexible marginal regression models on the recurrence process and the dependent censoring time without specifying their dependence structure. The proposed model generalizes the approach by Ghosh and Lin (2003, Biometrics 59, 877–885). The technique of artificial censoring provides a way to maintain the homogeneity of the hypothetical error variables under dependent censoring. Here we propose to apply this technique to two Gehan‐type statistics. One considers only order information for pairs whereas the other utilizes additional information of observed censoring times available for recurrence data. A model‐checking procedure is also proposed to assess the adequacy of the fitted model. The proposed estimators have good asymptotic properties. Their finite‐sample performances are examined via simulations. Finally, the proposed methods are applied to analyze the AIDS linked to the intravenous experiences cohort data.     

A quantile‐slicing approach for sufficient dimension reduction with censored responses

Sufficient dimension reduction (SDR) that effectively reduces the predictor dimension in regression has been popular in high‐dimensional data analysis. Under the presence of censoring, however, most existing SDR methods suffer. In this article, we propose a new algorithm to perform SDR with censored responses based on the quantile‐slicing scheme recently proposed by Kim et al. First, we estimate the conditional quantile function of the true survival time via the censored kernel quantile regression (Shin et al.) and then slice the data based on the estimated censored regression quantiles instead of the responses. Both simulated and real data analysis demonstrate promising performance of the proposed method.     

Analysis of multivariate recurrent event data with time‐dependent covariates and informative censoring

Multivariate recurrent event data are usually encountered in many clinical and longitudinal studies in which each study subject may experience multiple recurrent events. For the analysis of such data, most existing approaches have been proposed under the assumption that the censoring times are noninformative, which may not be true especially when the observation of recurrent events is terminated by a failure event. In this article, we consider regression analysis of multivariate recurrent event data with both time‐dependent and time‐independent covariates where the censoring times and the recurrent event process are allowed to be correlated via a frailty. The proposed joint model is flexible where both the distributions of censoring and frailty variables are left unspecified. We propose a pairwise pseudolikelihood approach and an estimating equation‐based approach for estimating coefficients of time‐dependent and time‐independent covariates, respectively. The large sample properties of the proposed estimates are established, while the finite‐sample properties are demonstrated by simulation studies. The proposed methods are applied to the analysis of a set of bivariate recurrent event data from a study of platelet transfusion reactions.     

Time‐dependent classification accuracy curve under marker‐dependent sampling

Evaluating the classification accuracy of a candidate biomarker signaling the onset of disease or disease status is essential for medical decision making. A good biomarker would accurately identify the patients who are likely to progress or die at a particular time in the future or who are in urgent need for active treatments. To assess the performance of a candidate biomarker, the receiver operating characteristic (ROC) curve and the area under the ROC curve (AUC) are commonly used. In many cases, the standard simple random sampling (SRS) design used for biomarker validation studies is costly and inefficient. In order to improve the efficiency and reduce the cost of biomarker validation, marker‐dependent sampling (MDS) may be used. In a MDS design, the selection of patients to assess true survival time is dependent on the result of a biomarker assay. In this article, we introduce a nonparametric estimator for time‐dependent AUC under a MDS design. The consistency and the asymptotic normality of the proposed estimator is established. Simulation shows the unbiasedness of the proposed estimator and a significant efficiency gain of the MDS design over the SRS design.     

Estimating differences in restricted mean lifetime using observational data subject to dependent censoring

In epidemiologic studies of time to an event, mean lifetime is often of direct interest. We propose methods to estimate group- (e.g., treatment-) specific differences in restricted mean lifetime for studies where treatment is not randomized and lifetimes are subject to both dependent and independent censoring. The proposed methods may be viewed as a hybrid of two general approaches to accounting for confounders. Specifically, treatment-specific proportional hazards models are employed to account for baseline covariates, while inverse probability of censoring weighting is used to accommodate time-dependent predictors of censoring. The average causal effect is then obtained by averaging over differences in fitted values based on the proportional hazards models. Large-sample properties of the proposed estimators are derived and simulation studies are conducted to assess their finite-sample applicability. We apply the proposed methods to liver wait list mortality data from the Scientific Registry of Transplant Recipients.     

Mind the gap: scale‐dependent and dispersal‐mediated response to forest fragmentation

Species richness of habitat fragments is affected by spatial isolation. However, the scale of this phenomenon, and its interactions with the species’ seed dispersal potential has remained underexplored. By integrating seed trap and species distribution data, Koh et al., in this issue of the Journal of Vegetation Science, make a compelling case for scale‐dependent species’ responses to forest fragmentation.     

A comparison of abundance estimates from extended batch‐marking and Jolly–Seber‐type experiments

Little attention has been paid to the use of multi‐sample batch‐marking studies, as it is generally assumed that an individual's capture history is necessary for fully efficient estimates. However, recently, Huggins et al. ( 2010 ) present a pseudo‐likelihood for a multi‐sample batch‐marking study where they used estimating equations to solve for survival and capture probabilities and then derived abundance estimates using a Horvitz–Thompson‐type estimator. We have developed and maximized the likelihood for batch‐marking studies. We use data simulated from a Jolly–Seber‐type study and convert this to what would have been obtained from an extended batch‐marking study. We compare our abundance estimates obtained from the Crosbie–Manly–Arnason–Schwarz (CMAS) model with those of the extended batch‐marking model to determine the efficiency of collecting and analyzing batch‐marking data. We found that estimates of abundance were similar for all three estimators: CMAS, Huggins, and our likelihood. Gains are made when using unique identifiers and employing the CMAS model in terms of precision; however, the likelihood typically had lower mean square error than the pseudo‐likelihood method of Huggins et al. ( 2010 ). When faced with designing a batch‐marking study, researchers can be confident in obtaining unbiased abundance estimators. Furthermore, they can design studies in order to reduce mean square error by manipulating capture probabilities and sample size.     

Estimating species richness: still a long way off!

Xu et al., in this issue of the Journal of Vegetation Science, compare several species richness estimators. All the non‐parametric estimators, such as Chao and jackknife estimators, underestimated the true number, whereas all the area‐based models, based on species–area curves, overestimated it. No reliable method yet exists to predict the number of species in an area that is appreciably larger than the one(s) sampled.     

Importance of plant‐ and microbe‐driven metabolic pathways for plant defence

Expression of plant phenotypes can depend on both plant genomes and interactions between plants and the microbes living in, on and near their roots. We understand a growing number of the mechanistic links between plant genotypes and phenotypes, such as defence against herbivory (see brief review in Hubbard et al., 2019), yet the links between root microbiomes and the comprehensive swathe of plant phenotypes they affect (Friesen et al., 2011) remain less clear. In this issue of Molecular Ecology, Hubbard et al. (2019) follow microbe‐ and plant‐driven changes in plant defence against hervibory from molecular underpinnings to ecological consequences, contrasting both the metabolites affected and the magnitude of defensive impact. Naively, we might expect plant genomes to drive more variation in phenotype than the root microbiome, but Hubbard et al. (2019) find the opposite, implying profound consequences for plant trait evolution and ecological interactions.     

Time to abandon the loss of dispersal ability hypothesis in island plants: A comment on García‐Verdugo,Mairal, Monroy,Sajeva and Caujapé‐Castells (2017)

García‐Verdugo et al. (2017) recently tested the loss of dispersal ability hypothesis in a wind‐dispersed shrub from Southern Europe. Although the hypothesis has guided research for over 150 years, García‐Verdugo et al. (2017) results failed to substantiate its central prediction—the loss of seed dispersal potential in island populations. Here, I highlight several additional limitations of the hypothesis. First, García‐Verdugo et al. (2017) results are not unusual. Empirical support for the hypothesis is equivocal. Second, when reduced dispersal potential is documented, it may often evolve as a passive by‐product of selection for large seeds, for reasons that are wholly unrelated to their dispersal. Third, the hypothesis does not readily apply to all plant dispersal modes, particularly plants that produce fleshy fruits. These issues advocate a fresh approach to the study of how selection shapes the evolution of dispersal potential on islands.     

Ecologically flexible endemics dominate Indo‐Pacific bird communities

Reeve et al. (2016, Ecography, 39 , 990‐997) found that ecologically flexible endemics dominate Indo‐Pacific bird communities. This negative relationship between local abundance and global range size contrasts strongly with the positive range size‐abundance relationship “rule,” which would predict community dominance by globally widespread species. Theuerkauf et al. (2017, Journal of Biogeography, 44 , 2161–2163) provide new data from New Caledonia which they claim invalidate our study. They find positive relationships between local abundance and local range size, which they attribute to endemic species having narrower habitat niches than globally widespread species. We reanalysed their data using global range sizes, corroborating the pattern we originally reported: negative relationships between local abundance and global range size, driven by a subset of adaptable endemic species. We stress the importance of being explicit about the scale of ecological mechanisms, and ensuring that the scale of analysis matches the scale of interpretation.     

Inference for low‐ and high‐dimensional multigroup repeated measures designs with unequal covariance matrices

We propose tests for main and simple treatment effects, time effects, as well as treatment by time interactions in possibly high‐dimensional multigroup repeated measures designs. The proposed inference procedures extend the work by Brunner et al. (2012) from two to several treatment groups and remain valid for unbalanced data and under unequal covariance matrices. In addition to showing consistency when sample size and dimension tend to infinity at the same rate, we provide finite sample approximations and evaluate their performance in a simulation study, demonstrating better maintenance of the nominal α‐level than the popular Box‐Greenhouse–Geisser and Huynh–Feldt methods, and a gain in power for informatively increasing dimension. Application is illustrated using electroencephalography (EEG) data from a neurological study involving patients with Alzheimer's disease and other cognitive impairments.     

On the validity of time‐dependent AUC estimation in the presence of cure fraction

During the last decades, several approaches have been proposed to estimate the time‐dependent area under the receiver operating characteristic curve (AUC) of risk tools derived from survival data. The validity of these estimators relies on some regularity assumptions among which a survival function being proper. In practice, this assumption is not always satisfied because a fraction of the population may not be susceptible to experience the event of interest even for long follow‐up. Studying the sensitivity of the proposed estimators to the violation of this assumption is of substantial interest. In this paper, we investigate the performance of a nonparametric simple estimator, developed for classical survival data, in the case when the population exhibits a cure fraction. Motivated from the current practice of deriving risk tools in oncology and cardiovascular disease prevention, we also assess the loss, in terms of predictive performance, when deriving risk tools from survival models that do not acknowledge the presence of cure. The simulation results show that the investigated method is valid even under the presence of cure. They also show that risk tools derived from survival models that ignore the presence of cure have smaller AUC compared to those derived from survival models that acknowledge the presence of cure. This was also attested with a real data analysis from a breast cancer study.     

Semiparametric analysis of recurrent events data in the presence of dependent censoring

Dependent censoring occurs in longitudinal studies of recurrent events when the censoring time depends on the potentially unobserved recurrent event times. To perform regression analysis in this setting, we propose a semiparametric joint model that formulates the marginal distributions of the recurrent event process and dependent censoring time through scale-change models, while leaving the distributional form and dependence structure unspecified. We derive consistent and asymptotically normal estimators for the regression parameters. We also develop graphical and numerical methods for assessing the adequacy of the proposed model. The finite-sample behavior of the new inference procedures is evaluated through simulation studies. An application to recurrent hospitalization data taken from a study of intravenous drug users is provided.     

Semiparametric regression estimation in the presence of dependent censoring

We propose a semiparametric estimation procedure for estimatingthe regression of an outcome Y, measured at the end of a fixedfollow-up period, on baseline explanatory variables X, measuredprior to start of follow-up, in the presence of dependent censoringgiven X. The proposed estimators are consistent when the dataare ‘missing at random’ but not ‘missing completelyat random’ (Rubin, 1976), and do not require full specificationof the complete data likelihood. Specifically, we assume thatthe probability of censoring at time t is independent of theoutcome Y conditional on the recorded history up to t of a vectorof time-dependent covariates that are correlated with Y. Ourestimators can be used to adjust for dependent censoring andnonrandom noncompliance in randomised trials studying the effectof a treatment on the mean of a response variable of interest.Even with independent censoring, our methods allow the investigatorto increase efficiency by exploiting the correlation of theoutcome with a vector of time-dependent covariates.     

Minimum‐Norm Estimation for Binormal Receiver Operating Characteristic (ROC) Curves

The receiver operating characteristic (ROC) curve is often used to assess the usefulness of a diagnostic test. We present a new method to estimate the parameters of a popular semi‐parametric ROC model, called the binormal model. Our method is based on minimization of the functional distance between two estimators of an unknown transformation postulated by the model, and has a simple, closed‐form solution. We study the asymptotics of our estimators, show via simulation that they compare favorably with existing estimators, and illustrate how covariates may be incorporated into the norm minimization framework.     

Robust approach to combining multiple markers to improve surrogacy

Identifying effective and valid surrogate markers to make inference about a treatment effect on long-term outcomes is an important step in improving the efficiency of clinical trials. Replacing a long-term outcome with short-term and/or cheaper surrogate markers can potentially shorten study duration and reduce trial costs. There is sizable statistical literature on methods to quantify the effectiveness of a single surrogate marker. Both parametric and nonparametric approaches have been well developed for different outcome types. However, when there are multiple markers available, methods for combining markers to construct a composite marker with improved surrogacy remain limited. In this paper, building on top of the optimal transformation framework of Wang et al. (2020), we propose a novel calibrated model fusion approach to optimally combine multiple markers to improve surrogacy. Specifically, we obtain two initial estimates of optimal composite scores of the markers based on two sets of models with one set approximating the underlying data distribution and the other directly approximating the optimal transformation function. We then estimate an optimal calibrated combination of the two estimated scores which ensures both validity of the final combined score and optimality with respect to the proportion of treatment effect explained by the final combined score. This approach is unique in that it identifies an optimal combination of the multiple surrogates without strictly relying on parametric assumptions while borrowing modeling strategies to avoid fully nonparametric estimation which is subject to the curse of dimensionality. Our identified optimal transformation can also be used to directly quantify the surrogacy of this identified combined score. Theoretical properties of the proposed estimators are derived, and the finite sample performance of the proposed method is evaluated through simulation studies. We further illustrate the proposed method using data from the Diabetes Prevention Program study.     

Giant taxon‐character matrices: the future of morphological systematics

Simões et al. ( 2016 ) argued that large matrices are linked to the construction of “problematic” characters, and that those characters negatively affect tree topology. In their re‐evaluation of two squamate datasets, however, Simões et al. ( 2016 ) simply eliminated what they termed “problematic” characters, rather than recode them. This practice ignores potential sources of phylogenetic information and, if it were to be more widely followed, would inhibit the advancement of the field of systematics. Here, we defend the necessity and inevitability of large morphological (phenomic) datasets and discuss best practices for morphological data collection in contemporary phylogenetics.     

A Parametric ROC Model‐Based Approach for Evaluating the Predictiveness of Continuous Markers in Case–Control Studies

Summary The predictiveness curve shows the population distribution of risk endowed by a marker or risk prediction model. It provides a means for assessing the model's capacity for stratifying the population according to risk. Methods for making inference about the predictiveness curve have been developed using cross‐sectional or cohort data. Here we consider inference based on case–control studies, which are far more common in practice. We investigate the relationship between the ROC curve and the predictiveness curve. Insights about their relationship provide alternative ROC interpretations for the predictiveness curve and for a previously proposed summary index of it. Next the relationship motivates ROC based methods for estimating the predictiveness curve. An important advantage of these methods over previously proposed methods is that they are rank invariant. In addition they provide a way of combining information across populations that have similar ROC curves but varying prevalence of the outcome. We apply the methods to prostate‐specific antigen (PSA), a marker for predicting risk of prostate cancer.