nseVarTd {hydromad}R Documentation

Time-delay corrected performance measure

Description

adjVarTd coalesces modelled flow peaks to observed flow peaks for each event separately. nseVarTd calculates Nash-Sutcliffe efficiency on the result using nseStat. Depending on the quality of the coalescing, this better indicate performance ignoring timing error.

Usage

nseVarTd(obs, mod, event,...)
adjVarTd(obs,mod,event,...)

Arguments

obs

observed data vector

mod

model-predicted data vector corresponding to obs.

event

zoo object of events, as returned by eventseq

...

Additional arguments to nseStat and estimateDelay.

Details

The success of this method in minimising the effect of timing error depends on how well modelled and observed peaks can be coalesced. This depends on:

Also note that large numbers of events will run slower.

Value

For nseVarTd, a single numeric value. For adjVarTd, a zoo object with the original modelled and observed data, the adjusted model output and the lag estimated for each event.

Author(s)

Joseph Guillaume

See Also

hydromad.stats,nseStat,objFunVal

Examples


data(Murrindindi)
x <- Murrindindi[1:100]
x <- merge(x,X=lag(x$Q,2))

event <- eventseq(x$P, thresh = 5, inthresh = 3.5, indur = 7, continue = TRUE)

nseStat(x$Q,x$X)
nseVarTd(x$Q,x$X,event,lag.max=3)

## Avoiding passing lag.max
hydromad.getOption("max.delay") ## Current setting - default is 10
hydromad.options(max.delay=3)

nseVarTd(x$Q,x$X,event)
hmadstat("r.sq.vartd")(x$Q,x$X,event=event)
[Package hydromad version 0.9-18 Index]