Description Usage Arguments Details See Also Examples
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
.
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, ...)
|
obj |
Output of function |
... |
Other parameter pass to |
lag.k, lag.max |
Lag size for the bootstrap blocks.
If |
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 |
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.
atSite
, MannKendall
, ar
.
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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.