powered by
model selection function taken from S3IB,
saut(Lv, pen, Kseq, n, seuil = sqrt(n)/log(n), biggest = TRUE)
likelihood
penalty
number of changes
number of datapoints
threshold
heuristic (biggest jump or slope)
a selected number of chagnes