Я использую функцию "zyp.trend.dataframe" пакета "zyp", чтобы вычислить наклон и получить значение p.value, равное "sig" (P-значение Кендалла, вычисленное для окончательных временных рядов с отклонениями), но я замечаю этопоследний отличается (sig = 0,2206714) от вычисленного в этих пакетах (p-значение = 0,06029): Манн-Кендалл, ModifiedMK или Trend (см. мой пример ниже)?
Знаете ли вы, как рассчитатьp.value (sig) пакета "zyp"?
Спасибо за ваше разъяснение
library(zyp)
zyp.trend.vector(c(0, 1, 3, 4, 2,5),method="yuepilon",conf.intervals=TRUE)
lbound trend trend ubound tau sig nruns autocor valid_frac линейный перехват
-0,5000000 1,0000000 6,0000000 1,5000000 0,6000000 0,2206714 1,0000000 -0,1666667 1,0000000 0,8285714 -1,0000000
library(Kendall)
MannKendall(c(0, 1, 3, 4, 2, 5))
тау = 0,733, 2-стороннее значение = 0,060289
library(modifiedmk)
mkttest(c(0, 1, 3, 4, 2, 5))
Z-значение угла наклонаS) P-значение Тау 1,87867287 1,00000000 11,00000000 28,33333333 0,06028917 0,73333333
library(trend)
sens.slope(c(0, 1, 3, 4, 2, 5))
Данные об уклоне Сена: c (0, 1, 3, 4, 2, 5) z = 1,8787, n = 6, p-значение= 0,06029 альтернативный чгипотеза: истина z не равна 0 95-процентный доверительный интервал: -0,5 2,0 выборочные оценки: наклон Сена 1