iter2G_2Ind.m

Z5 = SampleZ5(mu1(:,:,tt-1)',Sigma1(:,:,:,tt-1),mu2(:,:,tt-1)',Sigma2(:,:,:,tt-1),W1(:,tt-1),WWk(:,:,tt-1));

[Z1 Z2 ] = updateZ1Z2_2Ind(W1(:,tt-1),mu1(:,:,tt-1),Sigma1(:,:,:,tt-1),...
     WWk(:,:,tt-1),mu2(:,:,tt-1),Sigma2(:,:,:,tt-1),Z5);
 

%[U1(:,tt) W1(:,tt)] = updateUW_newG(Z1,Z4,J,alpha1(tt-1));
[U1(:,tt) W1(:,tt) accept_U1(:,tt)] = updateUW_2Ind(Z1,Z2,Z5,J,K,alpha1(tt-1),U1(:,tt-1),W1(:,tt-1), mu1(:,:,tt-1),Sigma1(:,:,:,tt-1),SPU,...
    mu2(:,:,tt-1),Sigma2(:,:,:,tt-1),WWk(:,:,tt-1));



alpha1(tt) = gamrnd(J+e1-1, 1/(f1-sum(log(1-U1(:,tt)))));

%[mu1(:,:,tt) Sigma1(:,:,:,tt)]=updatemu1Sigma1_newG(J,Z1,Z4,X1,...
 %   lambda1,m1,delta1,Phi1, Sigma1(:,:,:,tt-1),p1);

[mu1(:,:,tt) Sigma1(:,:,:,tt) accept_muSig1(:,tt)]=updatemu1Sigma1_2Ind(J,K,Z1,Z2,Z5,X1,lambda1,m1,delta1,Phi1,...
    mu1(:,:,tt-1),Sigma1(:,:,:,tt-1),SP,W1(:,tt),mu2(:,:,tt-1),Sigma2(:,:,:,tt-1),WWk(:,:,tt-1));

[V(:,tt) WW(:,tt) accept_V(:,tt)] = updateVWWG_2Ind(K,J,gamma2(tt-1),Vk(:,:,tt-1),alpha2, V(:,tt-1),WW(:,tt-1),SPV,mu2(:,:,tt-1),...
    Sigma2(:,:,:,tt-1));

gamma2(tt) = gamrnd(K+ee-1, 1/(ff-sum(log(1-V(:,tt)))));

[Vk(:,:,tt) WWk(:,:,tt) accept_V1(:,tt)] = updateVkWWk1_2Ind(Z5,Z2,WW(:,tt),alpha2,J,K,Vk(:,:,tt-1),mu1(:,:,tt),Sigma1(:,:,:,tt),W1(:,tt), mu2(:,:,tt-1),...
    Sigma2(:,:,:,tt-1), WWk(:,:,tt-1));

[mu2(:,:,tt) Sigma2(:,:,:,tt),C2] = updatemu2Sigma2_2Ind(Z5,x2,mu2(:,:,tt-1),Sigma2(:,:,:,tt-1),m2,PPII,phis,Sigma_mu_inv,delta2,Phi2,K);