## Create a time series of length 200 with three mean shifts at 50, 100, 150.
data = simgen(2, 1);
X = data$X[1, ]; ## time series
meta = data$meta; ## locations of metadata times
## Parameter estimation in the configuration where changepoints are time 50 and 99.
results = bcpmeta.parameters(X, meta = meta, eta = c(50, 99), trend = FALSE);
## Plot
cp.plot(results, meta.loc = -0.42, cex = 1.5);
Run the code above in your browser using DataLab