我在R中使用bfast来适应时间序列。我需要在它找到的每个段的斜率上提取斜率和错误。我想出了前者而不是后者。这是我如何得到斜坡：
fit<-bfast(my_ts,season="none",max.iter=2)
plot(fit,type="trend")
niter<-length(fit$output)
slopes<-coef(fit$output[[niter]]$bp.Vt)[,2]
slopes
但我似乎无法找到在斜坡上得到错误的方法。我发现有人建议使用
plot(fit,type="trend",ANOVA=T)
而不仅仅是绘图（拟合），但是a）显然，它只打印出图上的斜率+误差，我需要将它们输出到文件中，并且b）它对我来说无论如何都不适用，我明白了：
Error in out$bp.Wt$datatsp : $ operator is invalid for atomic vectors
In addition: Warning message:
In anova.lm(model) :
ANOVA F-tests on an essentially perfect fit are unreliable
您的帮助将不胜感激！