# NOT RUN {
## Simulate MVN data with 2 change points
data = cbind(1:180, mvn.sim.2.cps(180,8,seed=101))
## Change point detection using VCCP (it may take several minutes to complete...)
result = vccp.fun(data, "NBS", test = "V")
## Plot the result
getTestPlot(result)
result.2 = vccp.fun(data, "NBS", test = "B")
## Plot the result
getTestPlot(result.2)
# }
Run the code above in your browser using DataLab