import statsmodels.api as sm import numpy as np import statsmodels.formula.api as smf dat = sm.datasets.get_rdataset("Guerry", "HistData").data results = smf.ols('Lottery ~ Literacy + np.log(Pop1831)', data=dat).fit() print(results.summary())
