Try Intermediate Regression with statsmodels in Python
Download PDF