R bfastmonitor,如何提取断点的日历日期?
R bfastmonitor, how to extract the calendar date of the breakpoint?
我在光栅块上使用 bfastmonitor 来检测时间序列中的断点。
我正在为这样的区域使用 bfmSpOver:
bfm_area <- bfmSpOver(s, area_shape, fun=mean, start=c(2006, 1), formula = response~harmon+trend, order = 1, h=0.25 , 历史 = c("ROC"), 水平 = 0.05)
我的输出是:
bfm_area$breakpoint
[1] 2006.603
现在2006.603对应的是什么日期?到目前为止,我认为 bfm 输出以年份、儒略日显示,但这与此输出不匹配?
难道是 2006 年 8 月 8 日? 365 * 0.603 ~ 2006年第220天.
lubridate
包有一个函数 date_decimal
可以用来将十进制日期转换为日期对象。
我在光栅块上使用 bfastmonitor 来检测时间序列中的断点。
我正在为这样的区域使用 bfmSpOver:
bfm_area <- bfmSpOver(s, area_shape, fun=mean, start=c(2006, 1), formula = response~harmon+trend, order = 1, h=0.25 , 历史 = c("ROC"), 水平 = 0.05)
我的输出是:
bfm_area$breakpoint
[1] 2006.603
现在2006.603对应的是什么日期?到目前为止,我认为 bfm 输出以年份、儒略日显示,但这与此输出不匹配?
难道是 2006 年 8 月 8 日? 365 * 0.603 ~ 2006年第220天.
lubridate
包有一个函数 date_decimal
可以用来将十进制日期转换为日期对象。