simulateVAR(N = 100, p = 1, nobs = 250, rho = 0.5, sparsity = 0.05,
method = "normal", covariance = "toeplitz")"normal" or "bimodal"."toeplitz", "block1", "block2" or simply "diagonal".data a list with two elements: series the multivariate time series and
noises the time series of errors
S the variance/covariance matrix of the process