roll_quantile: Rolling Quantiles

In roll: Rolling and Expanding Statistics

Description

A function for computing the rolling and expanding quantiles of time-series data.

Usage

roll_quantile(x, width, weights = rep(1, width), p = 0.5,
  min_obs = width, complete_obs = FALSE, na_restore = FALSE,
  online = FALSE)

Arguments

x	vector or matrix. Rows are observations and columns are variables.
width	integer. Window size.
weights	vector. Weights for each observation within a window.
p	numeric. Probability between zero and one.
min_obs	integer. Minimum number of observations required to have a value within a window, otherwise result is NA.
complete_obs	logical. If TRUE then rows containing any missing values are removed, if FALSE then each value is used.
na_restore	logical. Should missing values be restored?
online	logical. Process observations using an online algorithm.

Details

The methodology for computing the quantiles is based on the inverse of the empirical distribution function with averaging at discontinuities (see "Definition 2" in Hyndman and Fan, 1996).

Value

An object of the same class and dimension as x with the rolling and expanding quantiles.

References

Hyndman, R.J. and Fan, Y. (1996). "Sample quantiles in statistical packages." American Statistician, 50(4), 361-365.

Examples

n <- 15
x <- rnorm(n)
weights <- 0.9 ^ (n:1)

# rolling quantiles with complete windows
roll_quantile(x, width = 5)

# rolling quantiles with partial windows
roll_quantile(x, width = 5, min_obs = 1)

# expanding quantiles with partial windows
roll_quantile(x, width = n, min_obs = 1)

# expanding quantiles with partial windows and weights
roll_quantile(x, width = n, min_obs = 1, weights = weights)

Example output

 [1]          NA          NA          NA          NA  0.03572019  0.35291800
 [7]  0.86705507  0.86705507  0.86705507  0.86705507  0.64867291 -0.65462721
[13] -0.65462721 -0.65462721  0.33727924
 [1] -0.18835217 -0.07631599  0.03572019  0.19431909  0.03572019  0.35291800
 [7]  0.86705507  0.86705507  0.86705507  0.86705507  0.64867291 -0.65462721
[13] -0.65462721 -0.65462721  0.33727924
 [1] -0.18835217 -0.07631599  0.03572019  0.19431909  0.03572019  0.19431909
 [7]  0.35291800  0.60998654  0.64867291  0.50079545  0.35291800  0.19431909
[13]  0.35291800  0.34509862  0.35291800
 [1] -0.18835217  0.03572019  0.03572019  0.35291800  0.03572019  0.35291800
 [7]  0.35291800  0.86705507  0.64867291  0.64867291  0.35291800  0.03572019
[13]  0.35291800  0.33727924  0.35291800

roll documentation built on July 13, 2020, 5:09 p.m.