Guest User

Untitled

a guest
Dec 12th, 2018
73
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.77 KB | None | 0 0
  1. desc_list = df_grouped['Medicare_Desc'].unique()
  2. brand_list = df_grouped['RptBrand'].unique()
  3.  
  4. dfs={}
  5. for b,d in [(b,d) for b in brand_list for d in desc_list]:
  6. series = df_grouped[(df_grouped['RptBrand'] == b) &
  7. (df_grouped['Medicare_Desc']== d)]
  8. series_trend = series.filter(items=['UtilPer1000'])
  9.  
  10. try:
  11. hp_cycle, hp_trend = sm.tsa.filters.hpfilter(series_trend, lamb=100)
  12. hp_trend.rename(columns = {'UtilPer1000': 'Trend'}, inplace=True)
  13.  
  14. rolling_avg = hp_trend.rolling(12).mean()
  15. rolling_avg.rename(columns = {'Trend': 'Rolling_avg'}, inplace=True)
  16.  
  17. except:
  18. pass
  19.  
  20. series = series.drop(['UtilPer1000', 'RptBrand'], axis=1)
  21. dfs[b] = pd.concat([series, hp_trend, rolling_avg], axis=1)
  22.  
  23. data=pd.concat(dfs, axis=0)
Add Comment
Please, Sign In to add comment