cumulate: augments a time series with net and gross sums of one data...

In delftfews: delftfews R extensions

Description

given the univariate 'input' timeseries, performs a discrete integration of its non-zero stretches of data.

Usage

cumulate(input, gap = 1, integration.method = 3, with.partials = FALSE, ...)
## S3 method for class 'zoo'
cumulate(input, gap = 1, integration.method = 3, with.partials =
FALSE, units="secs", skip.first = FALSE, ...)

Arguments

input	zoo; an univariate zoo series, indexed on POSIXct timestamps.
gap	numeric; the size of the minimal gap that interrupts gross stretches.
integration.method	numeric; refers to the overview presented in http://portal.acm.org/citation.cfm?id=578374, figure 7.2. 1: rectangular (top left), 2: rectangular (midpoint), 3: trapezoid, 4: simpson's..
units	character; this is passed to as.double for converting differences in timestamps. all units acceptable there can be used here.
with.partials	logic; whether the output should contain the partials used in the calculations. in this case, this will be the fifth column in the output.
skip.first	logic: whether an event at the very first timestamp should be skipped. you might consider an event on the first timestamp as suspect, as you probably missed the true beginning, occurring before the start of the available data.
...	placeholder for future extensions

Details

stretches are integrated "net" and "gross", where the gaps between "gross" stretches are at most 'gap' wide.

numeric values thus computed can be seen as events and they are located at the timestamp of the earliest non-zero entry of the stretch. all other entries are 'NA'.

it's best to specify 'data.description$missVal = NULL' when passing the result of 'cumulate.zoo to 'write.PI'. 'NA' events (no event) will be omitted.

Examples

input <- c(0, 0, 0, 5, 6, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0)

inputAsDelftfews <- timeseries(20580000, by=5, length.out=length(input), input=input)
cumulate(inputAsDelftfews, gap=3)

inputAsZoo <- zoo(input, order.by = structure(seq(20580000, 20580075, 5), class = c("POSIXct","POSIXt")))
cumulate(inputAsZoo, gap=3)

delftfews documentation built on May 2, 2019, 4:48 p.m.