import numpy as np
import wooldridge as woo
import statsmodels.formula.api as smf
ceosal1 = woo.dataWoo('ceosal1')
# estimate log-log model:
reg = smf.ols(formula='np.log(salary) ~ np.log(sales)', data=ceosal1)
results = reg.fit()
b = results.params
print(f'b: \n{b}\n')