Multivariate Surrogate Synchrony (mvSUSY) estimates the synchrony within datasets that contain more than two time series. mvSUSY was developed from Surrogate Synchrony (SUSY) with respect to implementing surrogate controls, and extends synchrony estimation to multivariate data. ‘mvSUSY’ works as described in Meier & Tschacher (2021).
#install.packages("mvSUSY") ## not yet on CRAN
## development version
install.packages("mvSUSY", repos=c("https://wtschacher.github.io/mvSUSY/","https://cloud.r-project.org"))Note that the following example assumes that the source data are in a flat file and it has particular structure (column names in first row, whitespace as field separator). If you do not have such, then use the command in the comment below to mockup random data.
library(mvSUSY)
## read in data from a flat file
data = read.csv(file.choose(), header=TRUE, sep=" ", na.strings=".")
## mockup random data if needed
#data = as.data.frame(replicate(5, sample(10, 5000, TRUE)))
## compute mvSUSY using 'lambda_max' method
res = mvsusy(data, segment=10, Hz=10)
res
## plot
plot(res, type="eigenvalue")
plot(res, type="density")
plot(res, type="free scale")
plot(res, type="segment-wise")
plot(res, type="time series")
## compute mvSUSY using 'omega' method
res = mvsusy(data, segment=10, Hz=10, method="omega")
res
plot(res, type="density")
plot(res, type="free scale")
plot(res, type="segment-wise")
plot(res, type="time series")
## export to flat file via data.frame and write.csv
df = as.data.frame(res)
df
Need a high-speed mirror for your open-source project?
Contact our mirror admin team at info@clientvps.com.
This archive is provided as a free public service to the community.
Proudly supported by infrastructure from VPSPulse , RxServers , BuyNumber , UnitVPS , OffshoreName and secure payment technology by ArionPay.