monitor_dailyStatistic: Calculate daily statistics

Description Usage Arguments Details Value Note Examples

View source: R/monitor_dailyStatistic.R

Description

Calculates daily statistics for each monitor in ws_monitor.

Usage

1
2
3
4
5
6
7
monitor_dailyStatistic(
  ws_monitor,
  FUN = get("mean"),
  dayStart = "midnight",
  na.rm = TRUE,
  minHours = 18
)

Arguments

ws_monitor

ws_monitor object

FUN

function used to collapse a day's worth of data into a single number for each monitor in the ws_monitor object

dayStart

one of sunset|midnight|sunrise

na.rm

logical value indicating whether NA values should be ignored

minHours

minimum number of valid data hours required to calculate each daily statistic

Details

Sunrise and sunset times are calculated based on the first monitor encountered. This should be accurate enough for all use cases involving co-located monitors. Monitors from different regions should have daily statistics calculated separately.

Value

A ws_monitor object with daily statistics for the local timezone.

Note

Note that the incoming ws_monitor object should have UTC (GMT) times and that this function calculates daily statistics based on local (clock) time. If you choose a date range based on UTC times this may result in an insufficient number of hours in the first and last daily records of the returned ws_monitor object.

The returned ws_monitor object has a daily time axis where each datetime is set to the beginning of each day, 00:00:00, local time.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
N_M <- monitor_subset(Northwest_Megafires, tlim=c(20150801,20150831))
WinthropID <- '530470010_01'
TwispID <- '530470009_01'
MethowValley <- monitor_subset(N_M,
                               tlim=c(20150801,20150831),
                               monitorIDs=c(WinthropID,TwispID))
MethowValley_dailyMean <- monitor_dailyStatistic(MethowValley,
                                                 FUN=get('mean'),
                                                 dayStart='midnight')
# Get the full Y scale
monitor_timeseriesPlot(MethowValley, style='gnats', col='transparent')
monitor_timeseriesPlot(MethowValley, monitorID=TwispID,
                       style='gnats', col='forestgreen', add=TRUE)
monitor_timeseriesPlot(MethowValley, monitorID=WinthropID,
                       style='gnats', col='purple', add=TRUE)
monitor_timeseriesPlot(MethowValley_dailyMean, monitorID=TwispID,
                       type='s', lwd=2, col='forestgreen', add=TRUE)
monitor_timeseriesPlot(MethowValley_dailyMean, monitorID=WinthropID,
                       type='s', lwd=2, col='purple', add=TRUE)
addAQILines()
addAQILegend("topleft", lwd=1, pch=NULL)
title("Winthrop & Twisp, Washington Daily Mean PM2.5, 2015")

PWFSLSmoke documentation built on Nov. 23, 2021, 5:06 p.m.