r - 寻找断点前后的斜率

Question

我目前正在 R 中使用 bfast 包来分解时间序列。我很好奇是否有可能在断点之前和之后(直到结束或另一个断点)提取趋势段的斜率？

下面的示例取自引用手册。



生成它的代码在这里。

require(bfast)
require(strucchange)
require(sandwich)
require(forecast)
require(raster)
require(sp)
fit <- bfast(harvest, season="harmonic", max.iter=2)
plot(fit, type="trend")

问题是结果输出对象不是自动的，换句话说，我可以找到断点之间的值并尝试根据这些趋势值制定斜率，但是这个过程非常耗时，因为我必须手动找到起点/end 断点值并提取其间的每个值。所以理想情况下，我想找到一种更简单的方法来识别多个时间序列的所有可用断点之前和之后的趋势斜率(蓝线)。

任何帮助表示赞赏，谢谢。

Answer 1

plot(fit, ANOVA = TRUE)

将为您提供“每个已识别趋势段的斜率和显着性值”。

有趣的是，当您指定时它确实有效:type="trend"