Critical_reorganization_of_adolescent_prefrontal_circuitry/code/LFP/PAC/MainSinglePAC.m

function Analysis= MainSinglePAC(rangePhase,rangeAmplitude,measure,Merge_num,Signal,fs)

% This is a helperfunction to calculate a matrix of PAC values using either the ESC, MI(MVL) 
% or CFC measure. 

% INPUTS:
%  fs - sampling frequency;
%  Signal - vector of 0.5-100 Hz filtered LFP signal
%  measure - measure to be used - it should be: 'esc', 'mi' or 'cfc';
%  rangePhase - range of frequencies for low signal; for example 1:1:20
%  rangeAmplitude - range of frequencies for high signal; for example 30:2:200
%  Merge_num - used 50

% Output
%PAC results for each analyzed Signal window, including:
    %pac.pacmat=pacmat;
    %pac.freqvec_ph---the phase frequency band;
    %pac.freqvec_amp---the amplitude frequency band;
    %pac.shf_data_mean---shuffle mean;
    %pac.shf_data_std----shuffle std;
    %pac.relat_mi-----Zscore of the pac.pacmat

% Author: Angela Onslow, May 2010
%% calc PAC
segs=size(Signal,1);
Merge_segs=floor(2*segs/Merge_num-1);
for Merge_seg=1:Merge_segs
    Tem=floor(Merge_seg-1)*Merge_num/2+1:(Merge_seg+1)*Merge_num/2;
    X1=Signal(Tem,:);
    X2=Signal(Tem,:);
    
    [PAC]=find_pac_shf_several_CutGlue_segs(X2,fs,measure,X1,rangePhase,rangeAmplitude);
    
    
    pacwin{Merge_seg}=PAC;
    clear PAC
end
Analysis.pacwin=pacwin;
clear Signal pacwin