Skip to contents

Extracts the p-value from the Generalized Mann-Kendall trend test.

Usage

get_MKp(X, level = 0.1, time_dependency_option = "AR1")

Arguments

X

Numeric vector of values to test for trend

level

numeric in (0,1), level of the test.

time_dependency_option

string, option for handling temporal dependence. Available :

  1. 'INDE', assume independence (i.e. the standard MK test)

  2. 'AR1', assumes AR1 short-term dependence structure (i.e. Hamed and Rao's version of the MK test)

  3. 'LTP', assume long-term persistence (i.e. Hamed's version of the MK test)

Value

Numeric p-value between 0 and 1 indicating the probability of observing the trend by chance if no true trend exists.

Details

Interpretation guidelines:

  • p ≤ 0.01: Very strong evidence against null hypothesis

  • 0.01 < p ≤ 0.05: Strong evidence

  • 0.05 < p ≤ 0.1: Weak evidence

  • p > 0.1: Little or no evidence

See also

GeneralMannKendall for the full test implementation

Examples

data(temperature)
#> Warning: data set ‘temperature’ not found
# Get p-value for temperature trend
pval <- get_MKp(temperature$T)
#> Error: object 'temperature' not found
print(paste("Trend significance:", ifelse(pval < 0.05, "Significant", "Not significant")))
#> Error: object 'pval' not found