R如何從趨勢包計算森斜率函數的截距？

Question

我想估計長期的降水趨勢（百分比）。 第一次，我從“趨勢包”中計算sen的斜率，但是它錯過了使用該線性方程式來查找每年的預測值的截距：

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

您能幫我如何計算f（S）方程的插值（B），以計算每年的S參數的預測值？

預先感謝您的幫助 ！

Answer 1

我不了解R'trend'軟件包不會生成Sen斜率的截距

您可以嘗試使用“ zyp”包，該包根據值的中位數生成截距 。

這是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