set.seed(1)
x=c(rnorm(50,0,1),rnorm(50,0,10),rnorm(50,0,5),rnorm(50,0,1))
out=cpt.var(x,penalty="Manual",value="2*log(n)",method="BinSeg",Q=5)
likelihood(out) # returns the raw likelihood (783.9144) and the likelihood + penalty (805.1076)
Run the code above in your browser using DataLab