Robust and Replicable Statistical Methods for Incomplete Data with Complex Observational Patterns