for fid, grp in grps:
    result = ols(y=grp.loc[:, ('MEAN')], x=grp.loc[:, (['Accum_Prcp',
        'Accum_HDD'])])
text_file = open('Output {}.txt'.format(fid), 'w')
Print OLS regression summary to text file