clc fiid = 75; bn = 100; fgl = length(fgtype); ocmat = cell(1,fgl); clrs = getstateclr; cx = -.01:.005:.115; hrn = [.02 .04]; pmmt = nan(bn,fgl); for fgi = 1:fgl gtype = fgtype{fgi}; load([savpath gtype '_alldata.mat']) fids = find(Fdx>0); cdata = Cdat(fids,:); cdata(cdata==0) = nan; fdata = Fdx(fids); vdata = Vdat(fids)/20000; bdata = Bst2(fids); if fgi==1; bdata = 2-bdata; end cln = size(cdata,1); for bi = 1:bn pid = datasample(1:cln,round(cln/4)); pdt = abs(vdata(pid)); pmd(bi,fgi) = sum(pdt>=hrn(1)&pdt<=hrn(2))/length(pdt); end end pmci = prctile(pmd,[2.5 97.5]); pmn = mean(pmd); [fdrp,rp] = multicomp_bt(pmd(:,2:end),pmd(:,1)) %% plotting figure(fiid);clf;hold all bw = []; bx = []; bbol = 1; msiz = 10; mtp = []; pclr = [.24*ones(1,3);.9 .5 .1;1 .8 .3;.1 0 .9;.5 .1 .56;0.8 .6 .8]; plot_bcibar(bx,pmci,pmn,pclr,msiz,mtp,bw,bbol) plotstandard set(gca,'xlim',[.3 6.7],'ylim',[0 .5],'ytick',0:.25:1,'yticklabel','') setfigsiz([-10.2000 520.6000 306.6000 280.4000])