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

Question

我想估计长期的降水趋势（百分比）。第一次，我从 "trend package" 计算 sen 的斜率，但它错过了应用此线性方程来查找每年的预测值的截距：

f(S) = QS + B (Q : sen's slope, S : parameter, B : intercept)

我用这个例子来说明我的问题：

# 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 包的文档，这是一个 dataframe 的例子

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

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

R How to calculate the intercept of sen's slope function from trend package?

statistics

r

function

intercept

trend