Measuring brain synchronization using PLV calculation of 2 subjects EEG recordnigs

Francois · July 23, 2019, 2:07pm

You would need first to learn how to use Brainstorm in general. Follow all the introduction tutorials at least until #19, using the example dataset provided:
https://neuroimage.usc.edu/brainstorm/Tutorials

Then you could read this advanced tutorial that deals more specifically with EEG:
https://neuroimage.usc.edu/brainstorm/Tutorials/Epilepsy

Finally, you will need to find a way to import your recordings into Brainstorm. I don't know what file formats you can export your data to. The list of supported file formats is available here:
https://neuroimage.usc.edu/brainstorm/Introduction#Supported_file_formats

When you have your recordings in the database, it will be very straightforward to compute PLV between pairs of signals.

If you don't want to go through all this training, you could try to reuse directly parts of the Matlab scripts for the PLV computation:

github.com

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


    FreqBands = [];            

    

% ==== PLV ====

case 'plv'

    bst_progress('text', sprintf('Calculating: PLV [%dx%d]...', size(sInputA.Data,1), size(sInputB.Data,1)));

    Comment = 'PLV: ';

    % Get frequency bands

    nFreqBands = size(OPTIONS.Freqs, 1);

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

    

    % ===== IMPLEMENTATION G.DUMAS =====

    % 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

        if isConnNN

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

            HA = hilbert_fcn(DataAband')';

            HB = HA;

        else