R如何计算趋势包中森的斜率函数的截距？

Question

我想估计长期的降水趋势（百分比）。在第一次，我从“趋势包”计算森林的斜率，但是它错过了截距以应用这个线性方程来找到每年的预测值：

  

f（S）= QS + B（Q：sen的斜率，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）方程的intercpt（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