PTE: How is the connectivity matrix stored?

Francois · October 14, 2020, 9:24am

Only the PTE data.
You can check how this is done directly in the code:

brainstorm-tools/brainstorm3/blob/master/toolbox/connectivity/bst_connectivity.m#L577


% Get frequency bands

nFreqBands = size(OPTIONS.Freqs, 1);

BandBounds = process_tf_bands('GetBounds', OPTIONS.Freqs);

% Intitialize returned matrix

R = zeros(size(sInputA.Data,1), size(sInputB.Data,1), nFreqBands);

% Loop on each frequency band

for iBand = 1:nFreqBands

    % Band-pass filter in one frequency band + Apply Hilbert transform

    DataAband = process_bandpass('Compute', sInputA.Data, sfreq, BandBounds(iBand,1), BandBounds(iBand,2), 'bst-hfilter-2019', OPTIONS.isMirror);

    % Compute PTE

    [dPTE, PTE] = PhaseTE_MF(permute(DataAband, [2 1]));

    if OPTIONS.isNormalized

        R(:,:,iBand) = dPTE;

        R(:,:,iBand) = R(:,:,iBand) - 0.5; % Center result around 0

    else

        R(:,:,iBand) = PTE;

    end

end

% We don't want to compute again the frequency bands

FreqBands = [];