The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

Readme file of VCCP package

VCCP

The goal of VCCP is to implement Vine Copula Change Points (VCCP) model to estimate the number and location of multiple change points in the functional connectivity structure of multivariate time series. Possible binary segmentation methods of the VCCP model include NBS, OBS, MOSUM, and WBS. You can also check the significance of detected candidates by Stationary Bootstrap or Vuong test method.

Installation

You can install the released version of VCCP from CRAN with:

install.packages("VCCP")

And the development version from GitHub with:

# install.packages("devtools")
devtools::install_github("xinxiong0238/VCCP")

Example

These are some basic examples which show you how to detect changes in the functional connectivity of multivariate time series using different binary segmentation methods and inference tests.

library(vccp)
## basic example code
data <- cbind(1:180, mvn.sim.2.cps(180, 8, seed = 101))
T <- 180

result.NV <- vccp.fun(data, method = "NBS", delta = 30, test = "V")
#> Binary search, round 1 ...
#> Binary search, round 2 ...
#> Binary search, round 3 ...
#> Perform Vuong test on candidates...
getTestPlot(result.NV)
title("VCCP: NBS + Vuong")

result.NB <- vccp.fun(data, method = "NBS", delta = 30, test = "B")
#> Binary search, round 1 ...
#> Binary search, round 2 ...
#> Binary search, round 3 ...
#> Perform stationary bootstrap test on candidates...
getTestPlot(result.NB)
title("VCCP: NBS + Stationary Bootstrap")

result.MV <- vccp.fun(data, method = "MOSUM", delta = 30, test = "V")
#> MOSUM search ...
#> Perform Vuong test on candidates...

getTestPlot(result.MV)
title("VCCP: MOSUM + Vuong")

result.MB <- vccp.fun(data, method = "MOSUM", delta = 30, test = "B")
#> MOSUM search ...
#> Perform stationary bootstrap test on candidates...

getTestPlot(result.MB)
title("VCCP: MOSUM + Stationary Bootstrap")

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.