为什么置信区间的绘制不起作用?

时间:2018-09-11 07:25:04

标签: r time-series intervals

我的数据包括时间和数量。我正在尝试使用 strucchange软件包来检测时间数据中的结构变化。为此,我遵循以下互联网页面上的说明:https://datascienceplus.com/structural-changes-in-global-warming/

我从级别结构更改部分开始。一切正常,除了当我想绘制观察和拟合的时间序列以及置信区间时,我得到以下结果:

enter image description here

期望的结果应该遍及整个x轴,而是集中在0左右。而且,我也不知道为什么我将“索引”作为x轴而不是时间数据。我不知道我在做什么错。

这是我的代码:

plot(df1, type ="l")
df1_win <- window(df1$amount)
lev_fit <- lm(df1_win ~ 1)
summary(lev_fit)
plot(df1_win)
lines(ts(fitted(lev_fit)))
df1_brk <- breakpoints(df1_win ~ 1, h = 0.1)
summary(df1_brk)
plot(df1_brk)
plot(df1_win)
lines(fitted(df1_brk, breaks = 5))
lines(confint(df1_brk, breaks = 5))
breakdates(df1_brk, breaks = 5)
coef(df1_brk, breaks = 5)

这是我的数据帧df1的摘录,因为实际上有74行。

df1

    date     amount

 2012-07-01 0.0000000
 2012-08-01 1.1111111
 2012-09-01 0.2985075
 2012-10-01 0.5141388
 2012-11-01 0.0000000
 2012-12-01 0.0000000
 2013-01-01 0.6849315
 2013-02-01 1.9762846
 2013-03-01 1.1799410
 2013-04-01 0.2881844
 2013-05-01 0.2617801
 2013-06-01 1.2285012
 2013-07-01 1.2285012
 2013-08-01 1.3539652
 2013-09-01 1.6694491
 2013-10-01 2.4000000
 2013-11-01 2.5065963
 2013-12-01 2.4869110
 2014-01-01 2.0497804
 2014-02-01 1.4044944
 2014-03-01 3.9443155
 2014-04-01 2.9748284
 2014-05-01 3.0623020
 2014-06-01 2.2044088
 2014-07-01 2.9686175
 2014-08-01 3.1304348
 2014-09-01 3.9028621
 2014-10-01 2.3942538
 2014-11-01 2.9021559
 2014-12-01 4.6280992
 2015-01-01 3.8616251
 2015-02-01 3.0252101
 2015-03-01 3.7565740
 2015-04-01 4.0977714 

我真的很感谢您提供帮助。

0 个答案:

没有答案