mkTrend: Modified Mann Kendall

View source: R/mkTrend_rcpp.R

mkTrend_rR Documentation

Modified Mann Kendall

Description

If valid observations <= 5, NA will be returned.

Usage

mkTrend_r(y, ci = 0.95, IsPlot = FALSE)

mkTrend(y, x = seq_along(y), ci = 0.95, IsPlot = FALSE)

Arguments

y

numeric vector

ci

critical value of autocorrelation

IsPlot

boolean

x

(optional) numeric vector

Details

mkTrend is 4-fold faster with .lm.fit.

Value

  • Z0 : The original (non corrected) Mann-Kendall test Z statistic.

  • pval0: The original (non corrected) Mann-Kendall test p-value

  • Z : The new Z statistic after applying the correction

  • pval : Corrected p-value after accounting for serial autocorrelation N/n*s Value of the correction factor, representing the quotient of the number of samples N divided by the effective sample size n*s

  • slp : Sen slope, The slope of the (linear) trend according to Sen test

Note

slp is significant, if pval < alpha.

Author(s)

Dongdong Kong

References

Hipel, K.W. and McLeod, A.I. (1994), Time Series Modelling of Water Resources and Environmental Systems. New York: Elsevier Science.

Libiseller, C. and Grimvall, A., (2002), Performance of partial Mann-Kendall tests for trend detection in the presence of covariates. Environmetrics 13, 71–84, \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1002/env.507")}.

See Also

fume::mktrend and trend::mk.test

Examples

x <- c(4.81, 4.17, 4.41, 3.59, 5.87, 3.83, 6.03, 4.89, 4.32, 4.69)
r <- mkTrend(x)
r_cpp <- mkTrend(x, IsPlot = TRUE)

rtrend documentation built on June 22, 2024, 11:39 a.m.