qm.empqm.doy: Short The actual QM function considering the DOY without...

Description Usage Arguments Details Value Author(s) Examples

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

Description

qm.empqm.doy is the actual QM function considering the day-of-year (DOY) and not considering frequency adaptation for precipitation.

Usage

1
2
qm.empqm.doy(series, doy.series, cdf.mod, cor.fun, method = c("linear",
  "bin"), var, output.diag = FALSE, minq, maxq, incq)

Arguments

series

Time series of modelled data, e.g. transient (control+) scenario series.

doy.series

Corresponding series of DOYs (output of function qm.doystring).

cdf.mod

CDF of the modelled time series in the calibration period, estimated with function qm.cdf.doy(...) (matrix[365,nquantiles]).

method

QM method: "binary" -> correction for closest quantile is used, "linear" -> linear interpolation of correction function between lower and upper quantile.

var

Variable to correct.

output.diag

Diagnostic QM console output (TRUE) or not (FALSE); default: FALSE.

minq

Minimum quantile for correction function [0.01 .. 0.99]; note: should not be 0 (for correct handling of extremes)!

maxq

Maximum quantile for correction function [0.01 .. 0.99]; note: should not be 1 (for correct handling of extremes)!

incq

Quantile increment for correction function (bin size).

corr.fun

Additive correction function (cdf A - cdf B), obtained by function qm.corfun.doy(...).

Details

The function estimates the index (quantile) in which an observation falls with respect to the calibration period. Based on this quantile (and the DOY) a correction is applied, the final quantile map. For the 29th February the same correction function as for 1st March (DOY=60) is applied. Old and new extremes beyond the min and max percentile considered are corrected according to the correction of the min and max percentile, respectively (first and last quantile considered therefore have to be 0.01 and 0.99 -> check in function qm.doqm).

Value

List of 3: $qm.input.series: Input series to QM (input argument series). $qm.corrected.series: Corrected series. $quantile.index: A series of quantile indices wrt. simulated series in calibration period.

Author(s)

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

Examples

1
2
3
4
5
6
## Not run: 
# Carry out DOY-dependent QM without frequency adaptation. Standard
# quantiles. Variable: precipitation. Linear QM, diagnostic terminal output.
qm.empqm.doy(series, doy.series, cdf.modelled, correction.function, 'linear', 'pr', TRUE, 0.01, 0.00, 0.01)

## End(Not run)

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