R中使用greenbrown/bfast包计算时序趋势出现无效区间数报错

问题原因 #

• 报错来自bfast 1.6.1版本与greenbrown 2.4.3的兼容性问题：bfast 1.6.1近期更新了季节分箱的底层计算逻辑，当输入frequency = 1的无季节周期年度时间序列时，分箱数计算会返回无效值，触发cut函数的invalid number of intervals报错。
• greenbrown的Trend()函数默认会自动调用bfast的季节检测逻辑，你传入的是年度数据、不需要做季节分解，没有手动关闭季节检测就会触发这个问题。

修复方案 #

最优解法（无需修改包版本） #

直接在调用Trend()时指定关闭季节检测即可，你的数据是年度序列，本身不存在季节分量，该参数设置完全符合数据特征：

Trend(mdata_ts, Seasonality = FALSE)

运行后即可正常输出年度趋势结果。

备选解法（回退bfast版本） #

如果后续处理带季节周期的序列也遇到同类报错，可以将bfast回退到更新前的稳定版本：

# 未安装devtools先执行安装
# install.packages("devtools")
devtools::install_version("bfast", version = "1.5.9", repos = "https://cran.r-project.org")

回退完成后重启R会话，重新加载包即可恢复之前的运行逻辑。

内容的提问来源于stack exchange，提问作者geo_dd