Редактировать: Это научит меня внимательно смотреть на вопрос, прежде чем я отвечу.Оператор BAR действительно для взаимодействия, а не для полиномиальных эффектов.
У Logistic пока нет средств для достижения этой цели, о которых я знаю, - но у glimmix есть экспериментальная методика, использующая утверждение эффекта.Например, это ..
effect MyPoly = polynomial(x1-x3/degree=2);
model y = MyPoly;
совпадает с
model y = x1 x2 x3 x1*x1 x1*x2 x1*x3 x2*x2 x2*x3 x3*x3;