defaults
%set(0,'DefaultAxesFontSize', 8);
%set(0,'DefaultTextFontSize', 8);
nruns=10;
totfiring = zeros(nruns, 12, 2);
totactive = zeros(nruns, 12, 2);
coract = zeros(nruns, 12);
CUTOFF=5; % Hz
diffs = [-240, -180, -120, -90, -60, -30, 30, 60, 90, 120, 180, 240];
CONDITION='pairstrongL'
brsyns = zeros(nruns,npyrs*nbranches);
brsynsH = zeros(nruns,10);
nrnsyns = zeros(nruns,npyrs);
nrnbins = [0:40]
nrnsynsH = zeros(nruns, length(nrnbins));
brcors = zeros(nruns, 9);
nrncors = zeros(nruns, 9);
for run=1:nruns
for ncase=1:length(diffs)
diff = diffs(ncase)
if (diff < 0)
p1 = 0;
p2 = -diff;
else
p1 = diff;
p2 = 0;
end
sfn=sprintf('./data/%s_%d_%d_%d/spikesperpattern.dat', CONDITION, p2, p1,run-1)
spk = load( sfn);
if (diff>0)
spk = flipud(spk);
end
spk = spk(:, 1:npyrs)/2.0;
pop = spk(spk>=CUTOFF);
totfiring(run, ncase, 1) = mean(spk(1, spk(1,:)>2), 2);
totfiring(run, ncase, 2) = mean(spk(2, spk(2,:)>2), 2);
totactive(run, ncase, :) = sum(spk>5,2);
s = corrcoef(spk');
coract(run, ncase) = s(1,2);
[bw, bs, nw, ns] = getsynstate(sprintf('./data/%s_%d_%d_%d/synstate.dat', CONDITION, p2, p1,run-1));
brcors(run,ncase) = corr(bw(1,:)', bw(2,:)');
nrncors(run,ncase) = corr(nw(1,:)', nw(2,:)');
% brsyns(run,:) = bs(1,:);
% brsynsH(run,:) = hist(bs(1,:), [0:9]);
% brsynsH(run,:) = brsynsH(run,:)/sum(brsynsH(run,:));
% nrnsyns(run,:) = ns(1,:);
% nrnsynsH(run,:) = hist(ns(1,:), nrnbins);
% nrnsynsH(run,:) = nrnsynsH(run,:)/sum(nrnsynsH(run,:));
end
end
close all
mf = mean(totfiring,1);
sf = std(totfiring,0,1);
figure
barwitherr(sf(1,:,1), mf(1,:,1));
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Strong Memory')
%ylabel('Avg firing rate [Hz]')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,15])
xlim([0,13])
savefig(sprintf('./figs/%s_fstrong',CONDITION), 'pdf')
figure
barwitherr(sf(1,:,2), mf(1,:,2));
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Weak Memory')
%ylabel('Avg firing rate [Hz]')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,15])
xlim([0,13])
savefig(sprintf('./figs/%s_fweak',CONDITION), 'pdf')
figure
mact = mean(totactive,1)
sact = std(totactive,0,1)
barwitherr(sact(1,:,1), mact(1,:,1));
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Strong Memory')
%ylabel('% Active neurons')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,70])
xlim([0,13])
savefig(sprintf('./figs/%s_actstrong',CONDITION), 'pdf')
figure
barwitherr(sact(1,:,2), mact(1,:,2)); % weak
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Weak Memory')
%ylabel('% Active neurons')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,70])
xlim([0,13])
savefig(sprintf('./figs/%s_actweak',CONDITION), 'pdf')
figure
barwitherr(std(coract, 0,1), mean(coract));
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Correlation between population firing rates')
%ylabel('Correlation coefficient')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,1.0])
xlim([0,13])
savefig(sprintf('./figs/%s_corr',CONDITION), 'pdf')
figure
barwitherr(std(brcors, 0,1), mean(brcors));
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Correlated synapses per branch')
%ylabel('Correlation')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,1.0])
xlim([0,13])
savefig(sprintf('./figs/%s_brcorr',CONDITION), 'pdf')
figure
barwitherr(std(nrncors, 0,1), mean(nrncors));
set(gca, 'XTick', [1:length(diffs)])
set(gca, 'XTickLabel', diffs)
title('Correlated synapses per neuron')
%ylabel('Correlation')
%xlabel('Weak-Strong Interval [minutes]')
ylim([0,1.0])
xlim([0,13])
savefig(sprintf('./figs/%s_nrnrcorr',CONDITION), 'pdf')
a = totactive(:, [1 2 3 10 11 12],2)
b = totactive(:, [4 5 6 7 8 9],2)
[h,p] = ttest2(a(:), b(:))
a = totfiring(:, [1 2 3 10 11 12],2)
b = totfiring(:, [4 5 6 7 8 9],2)
[h,p] = ttest2(a(:), b(:))