## ---- echo=FALSE---------------------------------------------------------
knitr::opts_chunk$set(fig.width=7, fig.height=5)

## ----Northwest_Megafires, message=FALSE----------------------------------
suppressPackageStartupMessages(library(PWFSLSmoke))

PacNW <- Northwest_Megafires
PacNW_24 <- monitor_rollingMean(PacNW, width=24)

## ----map-----------------------------------------------------------------
monitor_map(PacNW_24, slice=max)
addAQILegend(title="Max AQI", cex=0.7)

## ----leaflet-------------------------------------------------------------
# Commented out for the vignette
#monitor_leaflet(PacNW_24, slice=max)

## ----dygraph-------------------------------------------------------------
NezPerceIDs <- c("160571012_01","160690012_01","160690013_01",
                 "160690014_01","160490003_01","160491012_01")
NezPerce <- monitor_subset(PacNW, monitorIDs=NezPerceIDs)
monitor_timeseriesPlot(NezPerce, style='gnats')
addAQILines()
addAQILegend(cex=0.7)

## ----august--------------------------------------------------------------
PacNW <- monitor_subset(PacNW, tlim=c(20150801,20150831),
                        timezone="America/Los_Angeles")
PacNW_24 <- monitor_subset(PacNW_24, tlim=c(20150801,20150831),
                           timezone="America/Los_Angeles")
NezPerce <- monitor_subset(NezPerce, tlim=c(20150801,20150831),
                           timezone="America/Los_Angeles")

## ----dailyBarplot--------------------------------------------------------
layout(matrix(seq(6)))
par(mar=c(1,1,1,1))
for (monitorID in NezPerceIDs) {
  siteName <- NezPerce$meta[monitorID,'siteName']
  monitor_dailyBarplot(NezPerce, monitorID=monitorID, main=siteName, axes=FALSE) 
}
par(mar=c(5,4,4,2)+.1)
layout(1)

## ----acute---------------------------------------------------------------
data <- PacNW$data[,-1] # omit 'datetime' column
maxPM25 <- apply(data, 2, max, na.rm=TRUE) # maximum value at each site
worstAcute <- names(sort(maxPM25, decreasing=TRUE))[1:6] # monitorIDs for the six worst sites in PacNW
dplyr::intersect(worstAcute, NezPerceIDs)
PacNW$meta[worstAcute[1],c('siteName','countyName','stateCode')]

## ------------------------------------------------------------------------
PacNW_dailyAvg <- monitor_dailyStatistic(PacNW, FUN=mean, minHours=20)
data <- PacNW_dailyAvg$data[,-1]
unhealthyDays <- apply(data, 2, function(x) { sum(x >= AQI$breaks_24[4], na.rm=TRUE) })
worstChronic <- names(sort(unhealthyDays, decreasing=TRUE))[1:6]
dplyr::intersect(worstChronic, NezPerceIDs)
PacNW$meta[worstChronic[1],c('siteName','countyName','stateCode')]

## ---- echo=FALSE---------------------------------------------------------
knitr::opts_chunk$set(fig.width=5, fig.height=5)

## ----googleMap-----------------------------------------------------------
fireLons <- c(-118.461,-117.679,-120.039,-119.002,-119.662)
fireLats <- c(48.756,46.11,47.814,48.338,48.519)
monitor_esriMap(PacNW_dailyAvg, zoom=7, centerLon=-118, centerLat=47, slice=max)
addIcon('redFlame', fireLons, fireLats, expansion=.002)
addAQILegend(cex=0.8)
title("August, 2015", line=-1.5, cex.main=2)

## ---- echo=FALSE---------------------------------------------------------
knitr::opts_chunk$set(fig.width=7, fig.height=5)

## ----tribalMonitors_daily------------------------------------------------
Omak <- monitor_subset(PacNW, monitorIDs="530470013_01")
Kamiah <- monitor_subset(PacNW, monitorIDs="160490003_01")
layout(matrix(seq(2)))
monitor_dailyBarplot(Omak, main="August 2015 Daiy AQI -- Omak, WA",
                         labels_x_nudge=0.8, labels_y_nudge=250)
monitor_dailyBarplot(Kamiah, main="August Daily AQI -- Kamiah, ID",
                         labels_x_nudge=0.8, labels_y_nudge=250)
layout(1)

