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