检索绘制changepoints

Question

我已经绘制在R改为点作为

values <- c(1, 2, 3, 4, 5, 6, 7, 8, 3, 1, 2, 3) 
values.ts = ts(values, frequency = 12, start = c(2017, 1)) 
chpoints = cpt.mean(values.ts, method="PELT") 
cpts(chpoints) 
plot(chpoints) 

和绘制图形的坐标

的cpts(chpoints)返回[1] 4 8。

问题1：什么是[1] 4 8？

问题2：如何获得绘制的变化点坐标，红色线的x，y坐标？

Answer 1

问题1：cpts(chpoints)返回段末端位置的向量。以下示例：

​​

将添加两条垂直线，指示具有“稳定”平均值的段。没有。的细分 - 1个变化点（[1] 4 8），因为最后的“变化点”是您时间系列的最后一个值。

问题2：日期坐标的第i个破发可以发现：

index(values.ts)[cpts(chpoints)[i]] 

'稳定' 的平均值被存储在这里：

param.est(chpoints) # 2.50 6.50 2.25 

我还是不明白为什么会出现是变化点之后的空白（红线不包括整个时间序列）。但是，如果我们手动计算平均值：

mean(values.ts[1:cpts(chpoints)[1]]) # 
mean(values.ts[(cpts(chpoints)[1]+1):cpts(chpoints)[2]]) 
mean(values.ts[(cpts(chpoints)[2]+1):length(values.ts)]) 

我们收到2.50 6.50 2.25如上。