indmat_Yes_S1
hi=[];
hi=xf(im_x1);
hi1=[reshape(hi,coldata.Nt(1),2)  1-sum(reshape(hi,coldata.Nt(1),2)')'];

hi=xf(im_x2);
hi2=[reshape(hi,coldata.Nt(2),2)  1-sum(reshape(hi,coldata.Nt(2),2)')'];

hi=[hi1;hi2];
subplot(2,1,1)
plot(hi)
hold on

hi=[];
hi=xf(im_x21);

hi1=[[xf(im_xt); 1-sum(xf(im_xt))]';reshape(hi,coldata.Nt(3),2)  1-sum(reshape(hi,coldata.Nt(3),2)')'];

hi=xf(im_x22);
hi2=[reshape(hi,coldata.Nt(4),2)  1-sum(reshape(hi,coldata.Nt(4),2)')'];

hiA=[hi1;hi2];
% plot(hiA)



hi=xf(im_x23);
hi1=[reshape(hi,coldata.Nt(5),2)  1-sum(reshape(hi,coldata.Nt(5),2)')'];

hi=xf(im_x24);
hi2=[[reshape(hi,coldata.Nt(6),2)  1-sum(reshape(hi,coldata.Nt(6),2)')'];[xf(im_xr); 1-sum(xf(im_xr))]'];

hiA=[hiA;hi1;hi2];
subplot(2,1,2)
plot(hiA)
hold on

% 
% % load Trial.mat
% indmat1
% % xf=xf
% hi=[];
% hi=xf(im_x1);
% hi1=[reshape(hi,coldata.Nt(1),2)  1-sum(reshape(hi,coldata.Nt(1),2)')'];
% 
% hi=xf(im_x2);
% hi2=[reshape(hi,coldata.Nt(2),2)  1-sum(reshape(hi,coldata.Nt(2),2)')'];
% 
% hi=[hi1;hi2];
% subplot(1,2,1)
% 
% plot(hi(end:-1:1,:),1:sum(coldata.Nt([1 2])))
% ylim([0 sum(coldata.Nt([1:2]))])
% 
% hold on
% 
% hi=[];
% hi=xf(im_x21);
% 
% hi1=[[xf(im_xt); 1-sum(xf(im_xt))]';reshape(hi,coldata.Nt(3),2)  1-sum(reshape(hi,coldata.Nt(3),2)')'];
% 
% hi=xf(im_x22);
% hi2=[reshape(hi,coldata.Nt(4),2)  1-sum(reshape(hi,coldata.Nt(4),2)')'];
% 
% hiA=[hi1;hi2];
% % plot(hiA)
% 
% 
% 
% hi=xf(im_x23);
% hi1=[reshape(hi,coldata.Nt(5),2)  1-sum(reshape(hi,coldata.Nt(5),2)')'];
% 
% hi=xf(im_x24);
% hi2=[[reshape(hi,coldata.Nt(6),2)  1-sum(reshape(hi,coldata.Nt(6),2)')'];[xf(im_xr); 1-sum(xf(im_xr))]'];
% 
% hiA=[hiA;hi1;hi2];
% subplot(1,2,2)
% plot(hiA(end:-1:1,:),1:sum(coldata.Nt([3:6]))+2)
% 
% ylim([0 sum(coldata.Nt([3:6]))+2])
% 
% 
% hold on
% 
% 
