У меня есть следующая линейная регрессия:

import statsmodels.formula.api as sm

model = sm.ols(formula = 'a ~ b + c', data = data).fit()

Я хочу добавить квадратный термин для б в этой модели.

Есть ли простой способ сделать это с помощью statsmodels.ols? Есть ли лучший пакет, который я должен использовать, чтобы достичь этого?