partial.mk.test {trend} | R Documentation |
Partial Mann-Kendall Trend Test
Description
Performs a partial Mann-Kendall Trend Test
Usage
partial.mk.test(x, y, alternative = c("two.sided", "greater", "less"))
Arguments
x |
a "vector" or "ts" object that contains the variable, which is tested for trend (i.e. correlated with time) |
y |
a "vector" or "ts" object that contains the variable, which effect on "x" is partialled out |
alternative |
character, the alternative method; defaults to "two.sided" |
Details
According to Libiseller and Grimvall (2002), the test statistic
for x
with its covariate y
is
where the correlation is calculated as:
The conditional covariance between and
is
with
and
Value
A list with class "htest"
method |
a character string indicating the chosen test |
data.name |
a character string giving the name(s) of the data |
statistic |
the value of the test statistic |
estimate |
the Mann-Kendall score S, the variance varS and the correlation between x and y |
alternative |
a character string describing the alternative hypothesis |
p.value |
the p-value of the test |
null.value |
the null hypothesis |
Note
Current Version is for complete observations only. The test statistic is not corrected for ties.
References
Libiseller, C. and Grimvall, A., (2002). Performance of partial Mann-Kendall tests for trend detection in the presence of covariates. Environmetrics 13, 71–84, doi:10.1002/env.507.
See Also
Examples
data(maxau)
s <- maxau[,"s"]; Q <- maxau[,"Q"]
partial.mk.test(s,Q)