Multithreshold change plane model: Estimation theory and applications in subgroup identification

Subgroup analysis for longitudinal data based on a partial linear varying coefficient model with a change plane

Subgroup analysis has become an important tool to characterize the treatment effect heterogeneity, and finally towards precision medicine. On the other hand, longitudinal study is widespread in many fields, but subgroup analysis for this data type is still limited. In this article, we study a partial linear varying coefficient model with a change plane, in which the subgroups are defined based on linear combination of grouping variables, and the time-varying effects in different subgroups are estimated to capture the dynamic association between predictors and response. The varying coefficients are approximated by basis functions and the group indicator function is smoothed by kernel function, which are included in the generalized estimating equation for estimation. Asymptotic properties of the estimators for the varying coefficients, the constant coefficients and the change plane coefficients are established. Simulations are conducted to demonstrate the flexibility, efficiency and robustness of the proposed method. Based on the Standard and New Antiepileptic Drugs study, we successfully identify a subgroup in which patients are sensitive to the newer drug in a specific period of time.

Change-plane analysis for subgroup detection with a continuous treatment

Detecting and characterizing subgroups with differential effects of a binary treatment has been widely studied and led to improvements in patient outcomes and population risk management. Under the setting of a continuous treatment, however, such investigations remain scarce. We propose a semiparametric change-plane model and consequently a doubly robust test statistic for assessing the existence of two subgroups with differential treatment effects under a continuous treatment. The proposed testing procedure is valid when either the baseline function for the covariate effects or the generalized propensity score function for the continuous treatment is correctly specified. The asymptotic distributions of the test statistic under the null and local alternative hypotheses are established. When the null hypothesis of no subgroup is rejected, the change-plane parameters that define the subgroups can be estimated. This paper provides a unified framework of the change-plane method to handle various types of outcomes, including the exponential family of distributions and time-to-event outcomes. Additional extensions with nonparametric estimation approaches are also provided. We evaluate the performance of our proposed methods through extensive simulation studies under various scenarios. An application to the Health Effects of Arsenic Longitudinal Study with a continuous environmental exposure of arsenic is presented.

Change plane model averaging for subgroup identification

Central to personalized medicine and tailored therapies is discovering the subpopulations that account for treatment effect heterogeneity and are likely to benefit more from given interventions. In this article, we introduce a change plane model averaging method to identify subgroups characterized by linear combinations of predictive variables and multiple cut-offs. We first fit a sequence of statistical models, each incorporating the thresholding effect of one particular covariate. The estimation of submodels is accomplished through an iterative integration of a change point detection method and numerical optimization algorithms. A frequentist model averaging approach is then employed to linearly combine the submodels with optimal weights. Our approach can deal with high-dimensional settings involving enormous potential grouping variables by adopting the sparsity-inducing penalties. Simulation studies are conducted to investigate the prediction and subgrouping performance of the proposed method, with a comparison to various competing subgroup detection methods. Our method is applied to a dataset from a warfarin pharmacogenetics study, producing some new findings.

A reference-free R-learner for treatment recommendation

Assigning optimal treatments to individual patients based on their characteristics is the ultimate goal of precision medicine. Deriving evidence-based recommendations from observational data while considering the causal treatment effects and patient heterogeneity is a challenging task, especially in situations of multiple treatment options. Herein, we propose a reference-free R-learner based on a simplex algorithm for treatment recommendation. We showed through extensive simulation that the proposed method produced accurate recommendations that corresponded to optimal treatment outcomes, regardless of the reference group. We used the method to analyze data from the Systolic Blood Pressure Intervention Trial (SPRINT) and achieved recommendations consistent with the current clinical guidelines.

Multiply robust subgroup analysis based on a single-index threshold linear marginal model for longitudinal data with dropouts

Identifying subpopulations that may be sensitive to the specific treatment is an important step toward precision medicine. On the other hand, longitudinal data with dropouts is common in medical research, and subgroup analysis for this data type is still limited. In this paper, we consider a single-index threshold linear marginal model, which can be used simultaneously to identify subgroups with differential treatment effects based on linear combination of the selected biomarkers, estimate the treatment effects in different subgroups based on regression coefficients, and test the significance of the difference in treatment effects based on treatment-subgroup interaction. The regression parameters are estimated by solving a penalized smoothed generalized estimating equation and the selection bias caused by missingness is corrected by a multiply robust weighting matrix, which allows multiple missingness models to be taken account into estimation. The proposed estimator remains consistent when any model for missingness is correctly specified. Under regularity conditions, the asymptotic normality of the estimator is established. Simulation studies confirm the desirable finite-sample performance of the proposed method. As an application, we analyze the data from a clinical trial on pancreatic cancer.

Personalized treatment selection using observational data

Estimating the optimal treatment regime based on individual patient characteristics has been a topic of discussion in many forums. Advanced computational power has added momentum to this discussion over the last two decades and practitioners have been advocating the use of new methods in determining the best treatment. Treatments that are geared toward the 'best' outcome for a patient based on his/her genetic markers and characteristics are of high importance. In this article, we develop an approach to predict the optimal personalized treatment based on observational data. We have used inverse probability of treatment weighted machine learning methods to obtain score functions to predict the optimal treatment. Extensive simulation studies showed that our proposed method has desirable performance in selecting the optimal treatment. We provided a case study to examine the Statin use on cognitive function to illustrate the use of our proposed method.