function bayes_si_plot(data,si,a_si,xpt);

% function bayes_si_plot(data,si,a_si,xpt);

br = 1; %3

[nothing nsubj] = size(data);
%ax = .075*[-1 1 -1 1];
%tk = .075*[-1 0 1];
txtloc = [-.07,.06];
axlab = ['A_1,B_1';'A_1,B_2';'A_2,B_1';'A_2,B_2'];
msize_a = 12;
msize_b = 8;
lwid = 1;
cloud = .5*ones(1,3);
if strcmp(xpt,'pbr')
    xlab_b = '\pi_{ t}';
    ylab_b = '\pi_{ d}';
else
    xlab_b = '\pi_{ f}';
    ylab_b = '\pi_{ v}';
end
xlab_a = '\pi_{ p}';
ylab_a = '\pi_{ b}';

lenny = length(si);
rvec = randperm(lenny);
n_sub = floor(0.1*lenny);
si_sub = si(rvec(1:n_sub),:,:,:,:);
a_si_sub = a_si(rvec(1:n_sub),:,:,:);

si_x = [];
si_y = [];
for s = 1:nsubj
    si_x = [si_x; si(:,1,1,br,s)];
    si_y = [si_y; si(:,2,1,br,s)];
end
x_min = min(si_x);
x_max = max(si_x);
y_min = min(si_y);
y_max = max(si_y);
x_dif = abs(x_max - x_min);
y_dif = abs(y_max - y_min);
x_min = x_min - .05*x_dif;
x_max = x_max + .05*x_dif;
y_min = y_min - .05*y_dif;
y_max = y_max + .05*y_dif;
if x_dif > y_dif
    adj = (x_dif - y_dif)/2;
    y_min = y_min-adj;
    y_max = y_max+adj;
else
    adj = (y_dif - x_dif)/2;
    x_min = x_min-adj;
    x_max = x_max+adj;
end
ax = [x_min x_max y_min y_max];

figure
for s = 1:nsubj
    subplot(3,3,s)
    hold on
    axis square
    box on
    axis(ax)
    plot(si_sub(:,1,1,br,s),si_sub(:,2,1,br,s),'.',...
        'MarkerSize',msize_a,'MarkerEdgeColor',cloud)
    plot([-1 1],[0 0],'k:','LineWidth',lwid)
    plot([0 0],[-1 1],'k:','LineWidth',lwid)
    hpd_x = cilim(si(:,1,1,br,s));
    hpd_y = cilim(si(:,2,1,br,s));
    plot([hpd_x(1) hpd_x(1)],ax(3:4),'k')
    plot([hpd_x(2) hpd_x(2)],ax(3:4),'k')
    plot(ax(1:2),[hpd_y(1) hpd_y(1)],'k')
    plot(ax(1:2),[hpd_y(2) hpd_y(2)],'k')
    plot( ( data(s).pM(1,1,br) - ...
        ( ( data(s).pM(1,1,br) + data(s).pM(1,2,br) ) * ...
        ( data(s).pM(1,1,br) + data(s).pM(1,3,br) ) ) ), ...
        ( data(s).pM(2,2,br) - ...
        ( ( data(s).pM(2,1,br) + data(s).pM(2,2,br) ) * ...
        ( data(s).pM(2,2,br) + data(s).pM(2,4,br) ) ) ),'k+',...
        'MarkerSize',msize_b)
    plot( ( data(s).pM(1,1,br) - ...
        ( ( data(s).pM(1,1,br) + data(s).pM(1,2,br) ) * ...
        ( data(s).pM(1,1,br) + data(s).pM(1,3,br) ) ) ), ...
        ( data(s).pM(2,2,br) - ...
        ( ( data(s).pM(2,1,br) + data(s).pM(2,2,br) ) * ...
        ( data(s).pM(2,2,br) + data(s).pM(2,4,br) ) ) ),'ko',...
        'MarkerSize',msize_b)
    text(x_min+.1*x_dif,y_min+.1*y_dif,['S' num2str(s)])
    %text(txtloc(1),txtloc(2),['S' num2str(s)])
    %text(.58,-.75,axlab(1,:))
    %text(-.9,.75,axlab(2,:))
    if s == 1
        set(gca,'XTick',[],'YTick',[floor(y_max*100)/100])
        ylabel(ylab_a)
    elseif s == 4
        set(gca,'XTick',[],'YTick',[0])
        ylabel(ylab_a)
    elseif s == 8
        set(gca,'XTick',[0],'YTick',[])
        xlabel(xlab_a)
    elseif s == 7
        set(gca,'XTick',[ceil(x_min*100)/100],'YTick',[ceil(y_min*100)/100])
        ylabel(ylab_a)
        xlabel(xlab_a)
    elseif s == 2 || s == 3 || s == 5 || s == 6
        set(gca,'XTick',[],'YTick',[])
    end
end

subplot(3,3,9)
hold on
axis square
box on
axis(ax)
plot(a_si_sub(:,1,1,br),a_si_sub(:,2,1,br),'.',...
    'MarkerSize',msize_a,'MarkerEdgeColor',cloud)
plot([-1 1],[0 0],'k:','LineWidth',lwid)
plot([0 0],[-1 1],'k:','LineWidth',lwid)
hpd_x = cilim(a_si(:,1,1,br));
hpd_y = cilim(a_si(:,2,1,br));
plot([hpd_x(1) hpd_x(1)],ax(3:4),'k')
plot([hpd_x(2) hpd_x(2)],ax(3:4),'k')
plot(ax(1:2),[hpd_y(1) hpd_y(1)],'k')
plot(ax(1:2),[hpd_y(2) hpd_y(2)],'k')
%text(.58,-.75,axlab(1,:))
%text(-.9,.75,axlab(2,:))
%text(txtloc(1),txtloc(2),'\alpha')
set(gca,'XTick',[floor(x_max*100)/100],'YTick',[])
xlabel(xlab_a)
text(x_min+.1*x_dif,y_min+.1*y_dif,'\alpha')

axsize = .26;
axposx = [.15 .375 .6];
axposy = [.105 .405 .705];
for i = 1:9
   subplot(3,3,i)
   hax = gca;
   if i==1
       set(hax,'Position',[axposx(1) axposy(3) axsize axsize])
   elseif i==2
       set(hax,'Position',[axposx(2) axposy(3) axsize axsize])
   elseif i==3
       set(hax,'Position',[axposx(3) axposy(3) axsize axsize])
   elseif i==4
       set(hax,'Position',[axposx(1) axposy(2) axsize axsize])
   elseif i==5
       set(hax,'Position',[axposx(2) axposy(2) axsize axsize])
   elseif i==6
       set(hax,'Position',[axposx(3) axposy(2) axsize axsize])
   elseif i==7
       set(hax,'Position',[axposx(1) axposy(1) axsize axsize])
   elseif i==8
       set(hax,'Position',[axposx(2) axposy(1) axsize axsize])
   elseif i==9
       set(hax,'Position',[axposx(3) axposy(1) axsize axsize])
   end
end


si_x = [];
si_y = [];
for s = 1:nsubj
    si_x = [si_x; si(:,3,1,br,s)];
    si_y = [si_y; si(:,4,1,br,s)];
end
x_min = min(si_x);
x_max = max(si_x);
y_min = min(si_y);
y_max = max(si_y);
x_dif = abs(x_max - x_min);
y_dif = abs(y_max - y_min);
x_min = x_min - .05*x_dif;
x_max = x_max + .05*x_dif;
y_min = y_min - .05*y_dif;
y_max = y_max + .05*y_dif;
if x_dif > y_dif
    adj = (x_dif - y_dif)/2;
    y_min = y_min-adj;
    y_max = y_max+adj;
else
    adj = (y_dif - x_dif)/2;
    x_min = x_min-adj;
    x_max = x_max+adj;
end
ax = [x_min x_max y_min y_max];

figure
for s = 1:nsubj
    subplot(3,3,s)
    hold on
    axis square
    box on
    axis(ax)
    plot(si_sub(:,3,1,br,s),si_sub(:,4,1,br,s),'.',...
        'MarkerSize',msize_a,'MarkerEdgeColor',cloud)
    plot([-1 1],[0 0],'k:','LineWidth',lwid)
    plot([0 0],[-1 1],'k:','LineWidth',lwid)
    hpd_x = cilim(si(:,3,1,br,s));
    hpd_y = cilim(si(:,4,1,br,s));
    plot([hpd_x(1) hpd_x(1)],ax(3:4),'k')
    plot([hpd_x(2) hpd_x(2)],ax(3:4),'k')
    plot(ax(1:2),[hpd_y(1) hpd_y(1)],'k')
    plot(ax(1:2),[hpd_y(2) hpd_y(2)],'k')
    plot( ( data(s).pM(3,3,br) - ...
        ( ( data(s).pM(3,1,br) + data(s).pM(3,3,br) ) * ...
        ( data(s).pM(3,3,br) + data(s).pM(3,4,br) ) ) ), ...
        ( data(s).pM(4,4,br) - ...
        ( ( data(s).pM(4,4,br) + data(s).pM(4,2,br) ) * ...
        ( data(s).pM(4,4,br) + data(s).pM(4,3,br) ) ) ),'k+',...
        'MarkerSize',msize_b)
    plot( ( data(s).pM(3,3,br) - ...
        ( ( data(s).pM(3,1,br) + data(s).pM(3,3,br) ) * ...
        ( data(s).pM(3,3,br) + data(s).pM(3,4,br) ) ) ), ...
        ( data(s).pM(4,4,br) - ...
        ( ( data(s).pM(4,4,br) + data(s).pM(4,2,br) ) * ...
        ( data(s).pM(4,4,br) + data(s).pM(4,3,br) ) ) ),'ko',...
        'MarkerSize',msize_b)
    text(x_min+.1*x_dif,y_min+.1*y_dif,['S' num2str(s)])
    %text(txtloc(1),txtloc(2),['S' num2str(s)])
    %text(.58,-.75,axlab(1,:))
    %text(-.9,.75,axlab(2,:))    if s == 1 || s == 4
    if s == 1
        set(gca,'XTick',[],'YTick',[floor(y_max*100)/100])
        ylabel(ylab_b)
    elseif s == 4
        set(gca,'XTick',[],'YTick',[0])
        ylabel(ylab_b)
    elseif s == 8
        set(gca,'XTick',[0],'YTick',[])
        xlabel(xlab_b)
    elseif s == 7
        set(gca,'XTick',[ceil(x_min*100)/100],'YTick',[ceil(y_min*100)/100])
        ylabel(ylab_b)
        xlabel(xlab_b)
    elseif s == 2 || s == 3 || s == 5 || s == 6
        set(gca,'XTick',[],'YTick',[])
    end
end

subplot(3,3,9)
hold on
axis square
box on
axis(ax)
plot(a_si_sub(:,3,1,br),a_si_sub(:,4,1,br),'.',...
    'MarkerSize',msize_a,'MarkerEdgeColor',cloud)
plot([-1 1],[0 0],'k:','LineWidth',lwid)
plot([0 0],[-1 1],'k:','LineWidth',lwid)
hpd_x = cilim(a_si(:,3,1,br));
hpd_y = cilim(a_si(:,4,1,br));
plot([hpd_x(1) hpd_x(1)],ax(3:4),'k')
plot([hpd_x(2) hpd_x(2)],ax(3:4),'k')
plot(ax(1:2),[hpd_y(1) hpd_y(1)],'k')
plot(ax(1:2),[hpd_y(2) hpd_y(2)],'k')
text(x_min+.1*x_dif,y_min+.1*y_dif,'\alpha')
xlabel(xlab_b)
%text(.58,-.75,axlab(1,:))
%text(-.9,.75,axlab(2,:))
%text(txtloc(1),txtloc(2),'\alpha')
set(gca,'XTick',[floor(x_max*100)/100],'YTick',[])

axsize = .26;
axposx = [.15 .375 .6];
axposy = [.105 .405 .705];
for i = 1:9
   subplot(3,3,i)
   hax = gca;
   if i==1
       set(hax,'Position',[axposx(1) axposy(3) axsize axsize])
   elseif i==2
       set(hax,'Position',[axposx(2) axposy(3) axsize axsize])
   elseif i==3
       set(hax,'Position',[axposx(3) axposy(3) axsize axsize])
   elseif i==4
       set(hax,'Position',[axposx(1) axposy(2) axsize axsize])
   elseif i==5
       set(hax,'Position',[axposx(2) axposy(2) axsize axsize])
   elseif i==6
       set(hax,'Position',[axposx(3) axposy(2) axsize axsize])
   elseif i==7
       set(hax,'Position',[axposx(1) axposy(1) axsize axsize])
   elseif i==8
       set(hax,'Position',[axposx(2) axposy(1) axsize axsize])
   elseif i==9
       set(hax,'Position',[axposx(3) axposy(1) axsize axsize])
   end
end