function sc3tv(X,Z,itr,ival,j)


k0=itr(Z(itr)==0); k1=itr(Z(itr)==1); 
 
scatter3(X(j(1),k0),X(j(2),k0),X(j(3),k0),'b','filled')
hold on
scatter3(X(j(1),k1),X(j(2),k1),X(j(3),k1),'r','filled')

k0=ival(Z(ival)==0); k1=ival(Z(ival)==1); 
scatter3(X(j(1),k0),X(j(2),k0),X(j(3),k0),'b')
scatter3(X(j(1),k1),X(j(2),k1),X(j(3),k1),'r')

xlabel(['Metagene ',int2str(j(1))])
ylabel(['Metagene ',int2str(j(2))])
zlabel(['Metagene ',int2str(j(3))])
hold off