qm.cdf.doy: Compute CDF for each day-of-year with moving window

Description Usage Arguments Value Author(s) Examples

View source: R/qm.cdf.doy.R

Description

qm.cdf.doy computes the cumulative distribution function (CDF) for each day-of-year (DOY) with a moving window (+/- half window - 0.5). In case of precipitation also the wet day frequency and wet day CDF are computed. The 29th February is considered for all DOY ranges that include 1st March (DOY=60).

Usage

1
2
qm.cdf.doy(series, yearbegin, yearend, doywindow, minq, maxq, incq, var,
  pr.wet)

Arguments

series

Time series of modelled / observed data.

yearbegin

First year of series.

yearend

Last year of series.

doywindow

Total width of moving window [days].

minq

Minimum quantile [0..1].

maxq

Maximum quantile [0..1].

incq

Quantile increment (bin size).

var

Variable name ('pr' for precip).

pr.wet

Wet day threshold [mm/day].

Value

List of 3: $cdf.matrix: CDF for each DOY (matrix[365,nquantiles]). $f.wet: wet day frequency for each DOY (vector[365]); NA if not precipitation. $cdf.wet.vector: wet day CDF for each DOY (matrix[365,101]; all percentiles); NA if not precipitation.

Author(s)

Sven Kotlarski (MeteoSwiss), Jan Rajczak (ETH Zurich)

Examples

1
2
3
4
5
6
7
## Not run: 
# Compute DOY-dependent CDF for a daily series ranging from 1981 to 2010 and
# for min and max quantile of 0.01 and 0.99, respectively. Quantile
# increment: 0.01. Variable: Precipitation. Wet day threshold: 0.1 mm/day.
qm.cdf.doy(time.series,1981,2010,91,0.01,0.99,0.01,'pr',0.1)

## End(Not run)

SvenKotlarski/qmCH2018 documentation built on July 14, 2019, 7:39 p.m.