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

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