Optimal Two-Phase Sampling Designs for Generalized Raking Estimators With Multiple Parameters of Interest.

Large observational datasets, including those derived from electronic health records, are a valuable resource for medical research but are often affected by missingness, measurement error, and misclassification. Two-phase sampling with generalized raking (GR) estimation is an efficient and robust approach to statistical inference in such settings. In this approach, variables that are unavailable or measured with error in a large Phase 1 cohort are obtained with higher-quality measurements in a Phase 2 subsample. Previous research has studied optimal Phase 2 sampling designs for inverse-probability-weighted (IPW) estimators in nonadaptive, multiparameter settings, and for GR estimators in single-parameter settings. In this work, we extend these results by deriving optimal adaptive, multiwave sampling designs for IPW and GR estimators when multiple parameters are of interest. We propose several practical allocation strategies, examine the design choices that arise when implementing them, and evaluate their performance through extensive simulations and an example using data from people living with HIV attending the Vanderbilt Comprehensive Care Clinic (VCCC). Our results show that independently optimizing allocation for each parameter improves efficiency over traditional case-control sampling. We also derive an integer-valued, A-optimal allocation method that typically outperforms independent optimization. Notably, we find that optimal designs for GR can differ substantially from those for IPW, and that this distinction can meaningfully affect estimator efficiency in the multiple-parameter setting. These findings offer practical guidance for future two-phase studies involving incomplete or error-prone data.