import statsmodels.api as sm
import statsmodels.formula.api as smf
sm.show_versions()
data = pd.read_csv('sample.csv')
formula = 'OBJ ~ EXP1 + EXP2 + EXP3'
model = smf.glm(formula=formula, data=data, family=sm.families.Binomial()).fit()
print(model.summary())
print(model.conf_int(alpha=0.05))