檢索繪制的變更點的坐標

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 。