estimateDelayFrac {hydromad}R Documentation

Estimate the dead time between input and output, with a fractional component (redistribution of the input)

Description

Optimises the delay TDopt using lagFrac to maximise the correlation between the delayed input time series and (rises in) the corresponding time series

Usage

estimateDelayFrac(DATA, rises = TRUE, lag.max = hydromad.getOption("max.delay"))

Arguments

DATA

a ts-like object with named components:

U

input (forcing) time series.

Q

output (response) time series.

rises

use only rises in the output to estimate delay.

lag.max

largest delay (in time steps) to consider.

Details

estimateDelay may be preferred if there's no good reason to have a fractional lag/redistribution of the input across days.

Value

The estimated delay as an integer number of time steps.

Author(s)

Joseph Guillaume

See Also

estimateDelay, lagFrac

Examples

L<-0.6 ## Lag of 0.6
P <- c(2,0,5,1,6,10,0,0,0)
V1<-lagFrac(P,L)

estimateDelay(cbind(P,V1), rises = FALSE)
estimateDelayFrac(cbind(U=P,Q=V1),lag.max=5,rises=FALSE)
[Package hydromad version 0.9-18 Index]