st.chisqprob = lambda chisq, df: st.chi2.sf(chisq, df)cols = df_constant.columns[:-1]model = sm.Logit(df.TenYearCHD, df_constant[cols])result = model.fit()result.summary()