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