changes to analysis.py, using statsmodels