mannKendall: Mann-Kendall test for trend

Description Usage Arguments Details See Also Examples

Description

Return the results of the Mann-Kendall test for trends. In case of autocorrelation, the p-value of the test is obtained by block bootstrap. This function is a wrapper for the function MannKendall.

Usage

1
2
3
4
5
6
7
mannKendall(obj, ...)

## S3 method for class 'atSite'
mannKendall(obj, lag.k = NULL, lag.max = 3, nboot = 1000,
  use.rank = TRUE, ...)

mannKendallPeaks(obj, ...)

Arguments

obj

Output of function atSite or fitPot.

...

Other parameter pass to tsboot.

lag.k, lag.max

Lag size for the bootstrap blocks. If lag.k is NULL it is automatically selected and force to be inferior to lag.max.

nboot

Number of bootstrap sample.

use.rank

Does ranks should be used to transform data into uniform [0,1] variable or the model parameter should be used. Necessary for the selection of lag.k.

Details

The automatic selection of the block bootstrap is chosen according to the best order of an autoregressive model using the AIC criteria and Yule-Walker estimates. See ar for more details. Previously, the ranks of the original data are transformed into normal variable using either the estimated distribution or the ranks.

The function mannKendallPeaks is performing the Mann-Kendall on the magnitude of the peaks for a POT analysis.

See Also

atSite, MannKendall, ar.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
library(lmomco)

# TS with trend
xset <- rlmomco(100, vec2par(c(100,20,0), 'gev')) + (1:100)/5

ax <- atSite(xset)
mannKendall(ax)

mannKendall(ax, lag.k = 0)
mannKendall(ax, lag.k = 3)

martindurocher/floodRFA documentation built on June 5, 2019, 8:44 p.m.