wolffd@0: function [classHatTest, probPos] = mixgauss_classifier_apply(mixgauss, testFeatures)
wolffd@0: 
wolffd@0: Bpos = mixgauss_prob(testFeatures, mixgauss.pos.mu, mixgauss.pos.Sigma, mixgauss.pos.prior);
wolffd@0: Bneg = mixgauss_prob(testFeatures, mixgauss.neg.mu, mixgauss.neg.Sigma,  mixgauss.neg.prior);
wolffd@0: prior_pos = mixgauss.priorC(1);
wolffd@0: prior_neg = mixgauss.priorC(2);
wolffd@0: post = normalize([Bpos * prior_pos; Bneg * prior_neg], 1);
wolffd@0: probPos = post(1,:)';
wolffd@0: [junk, classHatTest] = max(post);
wolffd@0: classHatTest(find(classHatTest==2))=0;
wolffd@0: classHatTest = classHatTest(:);