% For each unique electrode pair, derive a measure of coherence
        % between the recorded signals
        % Build up a model of statistical relationships between signals
        % across all channels
        for elecone = 1:8
            for electwo = 1:8
                
                [cohere_allpairs(i, j, elecone, electwo, :), cohere_freqs] = mscohere(epoch(elecone, :), epoch(electwo, :), [], [], [], 128);

            end
            
        end