Evaluating large-scale propensity score performance through real-world and synthetic data experiments

Abstract

<div class="section"> <a class="named-anchor" id="s1"> <!-- named anchor --> </a> <h5 class="section-title" id="d266634e152">Background</h5> <p id="d266634e154">Propensity score adjustment is a popular approach for confounding control in observational studies. Reliable frameworks are needed to determine relative propensity score performance in large-scale studies, and to establish optimal propensity score model selection methods. </p> </div><div class="section"> <a class="named-anchor" id="s2"> <!-- named anchor --> </a> <h5 class="section-title" id="d266634e157">Methods</h5> <p id="d266634e159">We detail a propensity score evaluation framework that includes synthetic and real-world data experiments. Our synthetic experimental design extends the ‘plasmode’ framework and simulates survival data under known effect sizes, and our real-world experiments use a set of negative control outcomes with presumed null effect sizes. In reproductions of two published cohort studies, we compare two propensity score estimation methods that contrast in their model selection approach: <span class="inline-formula"> <math id="IM1" overflow="scroll"> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </math> </span>-regularized regression that conducts a penalized likelihood regression, and the ‘high-dimensional propensity score’ (hdPS) that employs a univariate covariate screen. We evaluate methods on a range of outcome-dependent and outcome-independent metrics. </p> </div><div class="section"> <a class="named-anchor" id="s3"> <!-- named anchor --> </a> <h5 class="section-title" id="d266634e174">Results</h5> <p id="d266634e176"> <span class="inline-formula"> <math id="IM2" overflow="scroll"> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </math> </span>-regularization propensity score methods achieve superior model fit, covariate balance and negative control bias reduction compared with the hdPS. Simulation results are mixed and fluctuate with simulation parameters, revealing a limitation of simulation under the proportional hazards framework. Including regularization with the hdPS reduces commonly reported non-convergence issues but has little effect on propensity score performance. </p> </div><div class="section"> <a class="named-anchor" id="s4"> <!-- named anchor --> </a> <h5 class="section-title" id="d266634e191">Conclusions</h5> <p id="d266634e193"> <span class="inline-formula"> <math id="IM3" overflow="scroll"> <msub> <mrow> <mi>L</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </math> </span>-regularization incorporates all covariates simultaneously into the propensity score model and offers propensity score performance superior to the hdPS marginal screen. </p> </div>