autocov_to_spwcgc
Calculate pairwise-conditional frequency-domain MVGCs (spectral multivariate Granger causalites)
Syntax
[f,fres] = autocov_to_spwcgc(G,fres,useFFT)
Arguments
See also Common variable names and data structures.
input
G autocovariance sequence fres frequency resolution useFFT use FFT method for autocovariance transform (default: as for autocov_xform)
output
f spectral Granger causality matrix
Description
Returns the matrix f of pairwise-conditional frequency-domain (spectral) MVGCs
(where [ij] denotes omission of the ij-th variables) between all pairs of variables represented in G, for a stationary VAR process with autocovariance sequence G. The first index i of f is the target (causee) variable, the second j the source (causal) variable and the third indexes the frequency. See ref. [1] for details.
Spectral causality is calculated up to the Nyqvist frequency at a resolution fres. If fres is not supplied it is calculated optimally as the number of autocovariance lags. Call freqs = sfreqs(fres,fs), where fs is the sampling rate, to get a corresponding vector freqs of frequencies on [0,fs/2].
The useFFT flag specifies the algorithm used to transform the autocovariance sequence; see autocov_xform for details.
The caller should take note of any warnings issued by this function and test results with a call isbad|(f,false)|.
For details of the algorithm, see autocov_to_smvgc and [1].
References
[1] L. Barnett and A. K. Seth, The MVGC Multivariate Granger Causality Toolbox: A New Approach to Granger-causal Inference, J. Neurosci. Methods 223, 2014 [ preprint ].
See also
autocov_to_smvgc | autocov_to_pwcgc | autocov_to_var | var2trfun | autocov_xform | sfreqs | isbad