skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Propensity Score Modeling in Electronic Health Records with Time-to-Event Endpoints: Application to Kidney Transplantation
For large observational studies lacking a control group (unlike randomized controlled trials, RCT), propensity scores (PS) are often the method of choice to account for pre-treatment confounding in baseline characteristics, and thereby avoid substantial bias in treatment estimation. A vast majority of PS techniques focus on average treatment effect estimation, without any clear consensus on how to account for confounders, especially in a multiple treatment setting. Furthermore, for time-to event outcomes, the analytical framework is further complicated in presence of high censoring rates (sometimes, due to non-susceptibility of study units to a disease), imbalance between treatment groups, and clustered nature of the data (where, survival outcomes appear in groups). Motivated by a right-censored kidney transplantation dataset derived from the United Network of Organ Sharing (UNOS), we investigate and compare two recent promising PS procedures, (a) the generalized boosted model (GBM), and (b) the covariate-balancing propensity score (CBPS), in an attempt to decouple the causal effects of treatments (here, study subgroups, such as hepatitis C virus (HCV) positive/negative donors, and positive/negative recipients) on time to death of kidney recipients due to kidney failure, post transplantation. For estimation, we employ a 2-step procedure which addresses various complexities observed in the UNOS database within a unified paradigm. First, to adjust for the large number of confounders on the multiple sub-groups, we fit multinomial PS models via procedures (a) and (b). In the next stage, the estimated PS is incorporated into the likelihood of a semi-parametric cure rate Cox proportional hazard frailty model via inverse probability of treatment weighting, adjusted for multi-center clustering and excess censoring, Our data analysis reveals a more informative and superior performance of the full model in terms of treatment effect estimation, over sub-models that relaxes the various features of the event time dataset.  more » « less
Award ID(s):
1811245
PAR ID:
10337419
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Data Science
ISSN:
1680-743X
Page Range / eLocation ID:
188 to 208
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Propensity score weighting is a tool for causal inference to adjust for measured confounders in observational studies. In practice, data often present complex structures, such as clustering, which make propensity score modeling and estimation challenging. In addition, for clustered data, there may be unmeasured cluster-level covariates that are related to both the treatment assignment and outcome. When such unmeasured cluster-specific confounders exist and are omitted in the propensity score model, the subsequent propensity score adjustment may be biased. In this article, we propose a calibration technique for propensity score estimation under the latent ignorable treatment assignment mechanism, i. e., the treatment-outcome relationship is unconfounded given the observed covariates and the latent cluster-specific confounders. We impose novel balance constraints which imply exact balance of the observed confounders and the unobserved cluster-level confounders between the treatment groups. We show that the proposed calibrated propensity score weighting estimator is doubly robust in that it is consistent for the average treatment effect if either the propensity score model is correctly specified or the outcome follows a linear mixed effects model. Moreover, the proposed weighting method can be combined with sampling weights for an integrated solution to handle confounding and sampling designs for causal inference with clustered survey data. In simulation studies, we show that the proposed estimator is superior to other competitors. We estimate the effect of School Body Mass Index Screening on prevalence of overweight and obesity for elementary schools in Pennsylvania. 
    more » « less
  2. Abstract BackgroundWe hypothesized that alemtuzumab use is safe in pediatric kidney transplant recipients (KTRs) with equivalent long‐term outcomes compared to other induction agents. MethodsUsing pediatric kidney transplant recipient data in the UNOS database between January 1, 2000, and June 30, 2022, multivariate logistic regression, multivariable Cox regression, and survival analyses were utilized to estimate the likelihoods of 1st‐year and all‐time hospitalizations, acute rejection, CMV infection, delayed graft function (DGF), graft loss, and patient mortality among recipients of three common induction regimens (ATG, alemtuzumab, and basiliximab). ResultsThere were no differences in acute rejection or graft failure among induction or maintenance regimens. Basiliximab was associated with lower odds of DGF in deceased donor recipients (OR 0.77 [0.60–0.99],p = .04). Mortality was increased in patients treated with steroid‐containing maintenance (HR 1.3 [1.005–1.7]p = .045). Alemtuzumab induction correlated with less risk of CMV infection than ATG (OR 0.76 [0.59–0.99],p = .039). Steroid‐containing maintenance conferred lower rate of PTLD compared to steroid‐free maintenance (HR 0.59 [0.4–0.8]p = .001). Alemtuzumab was associated with less risk of hospitalization within 1 year (OR 0.79 [0.67–0.95]p = .012) and 5 years (HR 0.54 [0.46–0.65]p < .001) of transplantation. Steroid maintenance also decreased 5 years hospitalization risk (HR 0.78 [0.69–0.89]p < .001). ConclusionsPediatric KTRs may be safely treated with alemtuzumab induction without increased risk of acute rejection, DGF, graft loss, or patient mortality. The decreased risk of CMV infections and lower hospitalization rates compared to other agents make alemtuzumab an attractive choice for induction in pediatric KTRs, especially in those who cannot tolerate ATG. 
    more » « less
  3. BACKGROUND: Lung transplantation is the gold standard for a carefully selected patient population with end-stage lung disease. We sought to create a unique risk stratification model using only preoperative recipient data to predict one-year postoperative mortality during our pre-transplant assessment. METHODS: Data of lung transplant recipients at Houston Methodist Hospital (HMH) from 1/2009 to 12/2014 were extracted from the United Network for Organ Sharing (UNOS) database. Patients were randomly divided into development and validation cohorts. Cox proportional-hazards models were conducted. Variables associated with 1-year mortality post-transplant were assigned weights based on the beta coefficients, and risk scores were derived. Patients were stratified into low-, medium- and high-risk categories. Our model was validated using the validation dataset and data from other US transplant centers in the UNOS database RESULTS: We randomized 633 lung recipients from HMH into the development (n=317 patients) and validation cohort (n=316). One-year survival after transplant was significantly different among risk groups: 95% (low-risk), 84% (medium-risk), and 72% (high-risk) (p<0.001) with a C-statistic of 0.74. Patient survival in the validation cohort was also significantly different among risk groups (85%, 77% and 65%, respectively, p<0.001). Validation of the model with the UNOS dataset included 9,920 patients and found 1-year survival to be 91%, 86% and 82%, respectively (p < 0.001). CONCLUSIONS: Using only recipient data collected at the time of pre-listing evaluation, our simple scoring system has good discrimination power and can be a practical tool in the assessment and selection of potential lung transplant recipients. 
    more » « less
  4. Extensive literature has been proposed for the analysis of correlated survival data. Subjects within a cluster share some common characteristics, e.g., genetic and environmental factors, so their time-to-event outcomes are correlated. The frailty model under proportional hazards assumption has been widely applied for the analysis of clustered survival outcomes. However, the prediction performance of this method can be less satisfactory when the risk factors have complicated effects, e.g., nonlinear and interactive. To deal with these issues, we propose a neural network frailty Cox model that replaces the linear risk function with the output of a feed-forward neural network. The estimation is based on quasi-likelihood using Laplace approximation. A simulation study suggests that the proposed method has the best performance compared with existing methods. The method is applied to the clustered time-to-failure prediction within the kidney transplantation facility using the national kidney transplant registry data from the U.S. Organ Procurement and Transplantation Network. All computer programs are available at https://github.com/rivenzhou/deep_learning_clustered. 
    more » « less
  5. Key PointsA2 to B incompatible transplantation is not fully practiced in the country, and further policies should encourage centers to perform more blood incompatible transplants.Centers that currently practice A2 to B incompatible transplants should give priority to blood type B patients who are willing to accept an A organ. This will benefit Asian and Black patients. BackgroundThe rate of A2 to B incompatible (ABO-i) kidney transplant continues to be low despite measures in the new kidney allocation system (KAS) to facilitate such transplants. This study shows how the number of ABO-i transplants could increase if KAS policies were used to their fullest extent through a boost in ABO-i priority points. MethodTransplant outcomes were predicted using the Kidney Pancreas Simulated Allocation Model, preloaded with national data of 2010. We used this simulation to compare KAS with a new intervention in which priority equal to cPRA=100 has been given to blood type B candidates who are willing to accept an A blood type organ. ResultsThe number of Black recipients increased by 375 (from 35% of the total recipient population to 38.7%), the number of blood type B Blacks increased by 65 (from 8% of the total recipient population to 9%), and the number of blood type B Black patients receiving blood type A kidneys increased by 49 (from 2% of the total recipient population to 2.5%). The same change occurred for Asians, particularly blood type B Asians (from 0.54% of the total recipient population to 0.7%). The average wait time notably decreased by 27 days for blood type B Black patients. In the proposed scenario, 263 blood type B Black patients received a blood type A organ (2.5% of the total recipient population) while only 181 (1.1%) of such transplants were performed in 2021. These results signify a considerable opportunity loss of ABO-i transplants for Black patients. ConclusionsIf this policy was universally adopted, we would expect to see an overall increase in A2 to B transplantation, but in reality, not all centers perform ABO-i transplantation. Thus, adopting this policy would incentivize other centers to perform more subtyping of A-type kidneys, and it would increase access to organs for blood type B Asian and Black patients in centers where ABO-i transplantation already takes place. 
    more » « less