for i in range(0,12):
stat='SITE=="{stat_name}"'
sta=stat.format(stat_name=station[i])
print(sta)
#print(data.index.levels[0])
#print(data.query(sta)[value])
plotData(axes[i],data.index.levels[0],data.query(sta)[value],ptype)#注意下站点对不对应
#axes[i].set_title(station[i],fontsize=4,loc='left')
axes[i].annotate(station[i], xy=(0.8, 1.0), xycoords="axes fraction",fontsize=6)
for i in range(0,12):
sta=station[i]
data.loc[sta,'rmse_100'], data.loc[sta,'mape_100'],data.loc[sta,'R_100']= statis(np.array(df.loc[sta,'pm25_100']), np.array(df.loc[sta,'pm25']))
data.loc[sta,'rmse_50'],data.loc[sta,'mape_50'],data.loc[sta,'R_50'] = statis(np.array(df.loc[sta,'pm25_50']), np.array(df.loc[sta,'pm25']))
data.loc[sta,'rmse_10'],data.loc[sta,'mape_10'],data.loc[sta,'R_10'] = statis(np.array(df.loc[sta,'pm25_10']), np.array(df.loc[sta,'pm25']))
data.loc[sta,'rmse_1'],data.loc[sta,'mape_1'],data.loc[sta,'R_1'] = statis(np.array(df.loc[sta,'pm25_1']), np.array(df.loc[sta,'pm25']))
data.loc[sta,'rmse_100use'],data.loc[sta,'mape_100use'],data.loc[sta,'R_100use'] = statis(np.array(df.loc[sta,'pm25_100use']), np.array(df.loc[sta,'pm25']))
data.loc[sta,'rmse_50use'],data.loc[sta,'mape_50use'],data.loc[sta,'R_50use'] = statis(np.array(df.loc[sta,'pm25_50use']), np.array(df.loc[sta,'pm25']))
