Я пытаюсь изменить формулы модели, используя update.formula
из базового пакета stats
, но добавление нового термина продукта не соответствует ожиданиям.Например,
f <- y~x
update.formula(f, .~.*z)
Я ожидал, что это вернется:
y~x*z
Но он вернул:
y ~ x + z + x:z
Может кто-нибудь объяснить это?Есть ли способ получить ожидаемый результат?