R Как рассчитать точку пересечения функции наклона Сена из пакета трендов?

Я хочу оценить тенденции осадков (в процентах) в долгосрочной перспективе. В первый раз я рассчитываю наклон сена из «пакета трендов», но он пропускает точку пересечения, чтобы применить это линейное уравнение для нахождения прогнозируемого значения каждого года:

f(S) = QS + B (Q: наклон Сена, S: параметр, B: точка пересечения)

Я использую этот пример, чтобы объяснить мою проблему:

# import data
require(trend)
data("maxau")
maxau

# Estimate Annual trend with sen's slope

s <- maxau[,"s"]
sens.slope(s) # equal to  -0.2876139 / year

Можете ли вы помочь мне, как рассчитать точку пересечения (B) уравнения f (S), чтобы рассчитать прогнозируемое значение параметра S для каждого года?

Заранее спасибо за вашу помощь !


person tazrart    schedule 22.01.2018    source источник
comment
Не могли бы вы помочь мне, как рассчитать точку пересечения наклона сена из пакета трендов!   -  person tazrart    schedule 24.01.2018


Ответы (1)


я не понял, что пакет R 'trend' не генерирует пересечение наклона Сена

вы можете попробовать использовать пакет 'zyp', который генерирует перехват на основе медианы значений formula.

вот документация пакета zyp и это пример с фреймворком данных

library(zyp)
df=data.frame(x=c(1,2,3,4,5),y=c(3,6,8,1,9))
zyp.sen(y~x,df)

документация: https://cran.r-project.org/web/packages/zyp/zyp.pdf

person Italo Moletto Lobos    schedule 03.11.2018