Запуск мышей с формулой в качестве переменной: мгновенная оценка вместо последующей оценки?

Пакет R mice содержит следующий пример:

library("mice")
imp <- mice(nhanes)
fit <- with(data=imp,exp=lm(bmi~hyp+chl))

Мне нужен гибкий вызов with(), например:

model_formula <- bmi~hyp+chl
fit <- with(data=imp,exp=lm(model_formula))

Но это выкидывает Error in eval(predvars, data, env) : object 'bmi' not found. Искал похожие проблемы. Я обнаружил проблему со шкафом Помогите понять ошибку в функции, которую я определил в R. У меня сложилось впечатление, что при записи exp=lm(model_formula) выражение lm(model_formula) вычисляется мгновенно, но при записи exp = lm(bmi~hyp+chl) оно вычисляется не сразу - вместо этого определение будет происходить в функции with.mice()? И если да, то как я могу предотвратить мгновенную оценку?