Critical_reorganization_of_adolescent_prefrontal_circuitry/code/LFP/PAC/subfunctions/PACmeanALL.m

function []=PACmeanALL(Param,Path,Experiment,exp_num)

for inum= 1: length(exp_num)
    iexperiment=exp_num(inum);
    
filename1=Experiment(iexperiment).name;
filename2=Experiment(iexperiment).name2;


    
    
load( strcat(Path.output,filesep,'Mi_PAC OB Phase LEC Amplitude 1s no zeropad filter 1',filesep,filename1,'.mat'));
Phases(inum,:)=pacwin{1}.freqvec_ph;
Amplitudes(inum,:)=pacwin{1}.freqvec_amp;
for iseg=1:length(pacwin)
Miseg(iseg,:,:)=pacwin{iseg}.relat_mi;
end
Mi(inum,:,:)=mean(Miseg,1);

clear pacwin
load( strcat(Path.output,filesep,'Mi_PAC OB Phase LEC Amplitude 1s no zeropad filter 1',filesep,filename2,'.mat'));
Phases_nO(inum,:)=pacwin{1}.freqvec_ph;
Amplitudes_nO(inum,:)=pacwin{1}.freqvec_amp;
for iseg=1:length(pacwin)
Miseg_nO(iseg,:,:)=pacwin{iseg}.relat_mi;
end
size(Miseg_nO)
Mi_nO(inum,:,:)=mean(Miseg_nO,1);
size(Mi_nO)

Mi_sub(inum,:,:)=Mi(inum,:,:)-Mi_nO(inum,:,:);
Mi_div(inum,:,:)=Mi(inum,:,:)./Mi_nO(inum,:,:)

figure
subplot(2,2,1)
imagesc(Phases(inum,:),Amplitudes(inum,:),squeeze(Mi(inum,:,:)));axis xy;  
title('Baseline PAC OB-LEC')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
subplot(2,2,2)
imagesc(Phases_nO(inum,:),Amplitudes_nO(inum,:),squeeze(Mi_nO(inum,:,:)));axis xy;  
title('naris occlusion PAC OB-LEC')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
colormap(jet)
subplot(2,2,3)
imagesc(Phases(inum,:),Amplitudes(inum,:),squeeze(Mi_sub(inum,:,:)));axis xy;  
title(' PAC OB-LEC baseline - nO')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
subplot(2,2,4)
imagesc(Phases_nO(inum,:),Amplitudes_nO(inum,:),squeeze(Mi_div(inum,:,:)));axis xy;  
title('PAC OB-LEC baseline/nO')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
colormap(jet)


end

MeanMi=mean(Mi,1);
MeanMi_nO=mean(Mi_nO,1);
MeanMi_sub=mean(Mi_sub,1);
MeanMi_div=mean(Mi_div,1);

%permutation testing
c = Mi;
d = Mi_nO;


npermutes = 1000;

diff = mean(c-d);
diff_null = zeros(npermutes,1);
for permi = 1:npermutes
    anull = [c; d];
    ind = randperm(length(anull),length(anull)/2);
    bnull = anull(ind);
    anull(ind) = [];
    
    diff_null(permi) = mean(anull-bnull);
end

pMi = 1-sum(abs(diff)>=abs(diff_null))/npermutes



figure
subplot(2,2,1)
imagesc(Phases(1,:),Amplitudes(1,:),squeeze(MeanMi));axis xy;  
title('mean Baseline PAC OB-LEC')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
subplot(2,2,2)
imagesc(Phases_nO(1,:),Amplitudes_nO(1,:),squeeze(MeanMi_nO));axis xy; 
title('naris occlusion PAC OB-LEC')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
colormap(jet)
subplot(2,2,3)
imagesc(Phases(1,:),Amplitudes(1,:),squeeze(MeanMi_sub));axis xy;  
title('mean PAC OB-LEC baselin-nO')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
subplot(2,2,4)
imagesc(Phases_nO(1,:),Amplitudes_nO(1,:),squeeze(MeanMi_div));axis xy; 
title('mean PAC OB-LEC baseline/nO')
xlabel('OB phase frequency [Hz]')
ylabel('LEC amplitude frequency [Hz]')
colormap(jet)
end