检索绘制的变更点的坐标

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)返回段末端位置的向量。 按照您的示例：

plot(chpoints)
abline(v=index(values.ts)[cpts(chpoints)[1]], col="blue") # 4
abline(v=index(values.ts)[cpts(chpoints)[2]], col="green") # 8

将添加两条垂直线，以指示“稳定”均值段。 没有 段数-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 。