use formula
handle missing data in covariates with liptak