update with formula design; add covariates