# Ceteris Lab downloadable Python script
# Course: Fundamentals of Python for Financial Econometrics

import numpy as np, statsmodels.api as sm
rng=np.random.default_rng(25); x=rng.normal(size=800); p=1/(1+np.exp(-(-.5+1.2*x))); y=rng.binomial(1,p)
fit=sm.Logit(y,sm.add_constant(x)).fit(disp=False)
print(fit.params); print(fit.get_margeff().summary())