mpib_memoreye/behavioural_modeling/quickplotcompBerrors.m

function quickplotcompBerrors(theta,uno,dos,titulo,legyes,Aave);
    
colores=(colormap((summer(3))));
colores(1,:)=[0.9 0 0.7];


theta=[theta theta(1)];
uno=[uno; uno(1,:)];
dos=[dos; dos(1,:)];

dat{1}=uno;
dat{2}=dos;
    legcoun=2;
        for i = 1:2

        ap=nanmean(dat{i},2);
        chance=ap;
    chance(:)=0.5;
        sem=nanstd(dat{i},1,2)/sqrt(size(dat{i},2));
        uppererror=ap+sem;
        lowerror=ap-sem;

        pl{1}=polarplot(theta,chance,'k--','LineWidth',2);hold on;
        polarplot(theta,ap,'-','Color',colores(1+i-1,:),'LineWidth',2.5);hold on;

        pl{legcoun}=polarplot(theta,ap,'o','Color',colores(1+i-1,:),'MarkerFaceColor',colores(1+i-1,:),'MarkerSize',8); hold on;
        legcoun=legcoun+1;
        if legcoun~=3
        pl{legcoun}=polarplot(theta,uppererror,'--','Color',colores(1+i-1,:),'LineWidth',2);hold on;

        polarplot(theta,lowerror,'--','Color',colores(1+i-1,:),'LineWidth',2);hold on;
        end
        legcoun=legcoun+1;
        end
    
    if legyes
    legend([pl{1} pl{2} pl{4} pl{5}],'unbiased (B=0)',['model (B=' num2str(Aave) ')'],'mean behaviour','SEM behaviour','Location','SouthOutside')
    end
        
    rlim([0 0.7]);
    title(titulo); %legend('model','behav','Chance')
    
end