add a where argument

[markvanderloo]

impute_lm(y~x|z, where x > 0)

would first select records where x>0, fit the model, impute, and return the whole dataset.