My blog post details the process of loading the esoph dataset into Python using statsmodels.
The ipython notebook file in this repository contains the full analysis. You can see, for example, how to make the following figures:
agalea91/esoph_linear_model
Importing esoph into python from R datasets and building multivariate linear models
Jupyter Notebook