¿Cómo agregar una restricción de suma a cero a GLM en Python?
Tengo un modelo configurado en Python usando el modelo statsglm
función pero ahora quiero agregar una restricción de suma a cero al modelo.
El modelo se define de la siguiente manera:
import statsmodels.formula.api as smf
model = smf.glm(formula="A ~ B + C + D", data=data, family=sm.families.Poisson()).fit()
En R, para agregar la restricción, simplemente haría algo como esto:
model <- glm(A ~ B + C + D –1, family=poisson(), data=data, contrasts=list(C="contr.sum", D="contr.sum"))
Eso agrega la restricción de suma a cero a C y D, pero no estoy seguro de cómo lograr lo mismo en Python.
He visto que hay unfit_constraint()
método disponible, pero no estoy muy seguro de cómo usarlo o si es incluso lo correcto para lograr lo que necesito.
¿Alguien puede ofrecer algún consejo para aplicar esta restricción?