R中的多类变更点检测

Question

我的分类预测模型的输出如下：

a <- c(1,1.1,1,1,2,0.9,1.1,1,1.1,1) ## Class A
b <- c(2,2.1,1.9,1.7,2,2,3,2,2,2)   ## Class B
c <- c(3,3,3.1,3.6,3.2,2,3.1,3,3,3) ## Class C
x <- data.frame(c(a,b,c))
x$color <- rep(c("red","green","blue"),times=c(10,10,10))

我试图找出班级类型改变的位置。在这种情况下，类类型在位置11和21处更改。

我尝试了以下套餐：CPM, EPC, QCC但它没有给我我想要的东西。

检测班级类型发生变化的最佳方法是什么？

Answer 1

怎么样？

library(strucchange)
( bp <- breakpoints( x[, 1] ~ 1) )
#    Optimal 3-segment partition: 
# 
# Call:
# breakpoints.formula(formula = x[, 1] ~ 1)
# 
# Breakpoints at observation number:
# 10 20 
# 
# Corresponding to breakdates:
# 0.3333333 0.6666667 

bp$breakpoints
# [1] 10 20

（via）