The goal of ocp is to implement Bayesian Online Changepoint Detection, as described: https://arxiv.org/abs/0710.3742
This is a basic example of how to use the function “onlineCPD” on simulated univariate Gaussian data as input.
library(ocp)
# the true changepoint locations including the first and last point
<- c(1, 51, 71, 121)
truecps#simulate the data
set.seed(1)
<- c(rnorm(n=diff(truecps)[1], mean=0, sd=2),
uvgrnorm(n=diff(truecps)[2], mean=20, sd=4),
rnorm(n=diff(truecps)[3], mean=10, sd=3))
<- onlineCPD(uvg) ocpd_output