Affiliations: [a] Division of Biostatistics, Department of Population Health, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, USA | [b] Department of Applied Mathematics and Cybernetics, Tomsk State University, Lenin St., 36, Tomsk, 634050, Russia. E-mail: dmit@fpmk.tsu.ru | [c] Department of Statistics and Department of Biostatistics, University of Kentucky, Patterson Office Tower 839, Lexington, KY 40506-0027, USA. E-mail: kryscio@email.uky.edu
Abstract: This paper derives two classes of minimum variance estimators in the presence of missing data. One estimator represents an alternative expression to Weighted Least Square (WLS). Another estimator, hierarchical WLS or HWLS is identical to WLS for normal and categorical missing data. However, WLS and HWLS are different in general and this manuscript focuses on this new class of estimators. Application of HWLS to random dropout survival data yields the well-known Kaplan-Meier estimator. Simulation studies compare WLS, HWLS, and Expectation Maximization (EM) for the difference in means in a bivariate normal model. Another set of simulations compares WLS, HWLS, and Multiple Imputation (MI) under normality assumption for the difference in medians in a bivariate exponential model. Both Missing Completely At Random (MCAR) and Missing At Random (MAR) data are considered in these simulation studies. Finally, HWLS is applied to estimating parameters of a proportional hazards model in a well known mouse leukemia dataset with missing covariates.
Keywords: Parameter estimation, weighted least squares, missing and censored data, hierarchical structures
