drop_cov_from_formula adds drops some covariates from a
formula; covs should be a list of variable names
Examples
ff <- y ~ x + w + z
drop_cov_from_formula(list("w", "z"), ff)
#> y ~ x
#> <environment: 0x6519f90d0200>
drop_cov_from_formula("z", ff)
#> y ~ x + w
#> <environment: 0x6519f9113b80>
