#Estimate Binomial GLM#First instantiate a GLM model objectmodel=GLM(y,X,family=Binomial())#Set family to Binomial family object for Binomial GLM#Then use the fit method to estimate coefficients and compute diagnosticsresults=model.fit()

#Estimated prameters, intercept is always the first column on the leftprint(results.params)