Eu tenho a seguinte regressão linear:

import statsmodels.formula.api as sm

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

Eu quero adicionar um termo quadrático para b neste modelo.

Existe uma maneira simples de fazer isso com o statsmodels.ols? Existe um pacote melhor que eu deveria estar usando para conseguir isso?