modelos de estatísticas python - termo quadrático em regressão
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?