# NOT RUN {
# Spokane summer of 2015
airnow <- airnow_load(20150701,20150930)
airnow <- monitor_rollingMean(airnow, width=3)
MonroeSt <- monitor_subset(airnow, monitorIDs="530630047_01")
EBroadway <- monitor_subset(airnow, monitorIDs="530639997_01")
rocList <- skill_ROC(EBroadway, MonroeSt, t1Range=c(0,100), t2=55)
roc <- rocList$roc
auc <- rocList$auc
plot(roc$TPR ~ roc$FPR, type='S')
title(paste0('Area Under Curve = ', format(auc,digits=3)))
# }
Run the code above in your browser using DataLab