disp report

subplot(2, 3, 1);
disp mu;
disp(showface(mu));
subplot(2, 3, 4);
disp diag(psi);
disp(showface(diag(psi)));
subplot(2, 3, 2);
disp lambda_init_1;
disp(showface(lambda_init(:,1)));
subplot(2, 3, 5);
disp lambda_1;
disp(showface(lambda(:,1)));
subplot(2, 3, 3);
disp lambda_init_2;
disp(showface(lambda_init(:,2)));
subplot(2, 3, 6);
disp lambda_2;
disp(showface(lambda(:,2)));
input('faces 1');
subplot(1,1,1);

plot(likhist);
xlabel('iteration');
ylabel('shifted mean log likelihood');
input('trajectory');

hist(lik_train, bins);
xlabel('log likelihood on training set');
input('log likelihood on training set');

hist(lik_test, bins);
xlabel('log likelihood on test set');
input('log likelihood on test set');

plot(z_train(1,:), z_train(2,:), '.');
xlabel('z1');
ylabel('z2');
input('z2 over z1 for train');

[junk, i] = sort(z_train(2,:));
z2_index = i;
subplot(3, 4, 3);
showface(xx_train(:,i(n_train)));
subplot(3, 4, 7);
showface(xx_train(:,i(round(n_train/2))));
subplot(3, 4, 11);
showface(xx_train(:,i(1)));
[junk, i] = sort(z_train(1,:));
z1_index = i;
subplot(3, 4, 8);
showface(xx_train(:,i(n_train)));
subplot(3, 4, 6);
showface(xx_train(:,i(round(n_train/2))));
subplot(3, 4, 5);
showface(xx_train(:,i(1)));
input('faces 2');
subplot(1,1,1);

disp naive_threshold;
disp(naive_threshold);
disp naive_test_error;
disp(naive_error);
disp my_threshold;
disp(lik_threshold);
disp my_test_error;
disp(test_error);
disp my_train_error;
disp(train_error);
input('classification error');

input('images sorted over z1');
for i = 1:group_size:n_train
    showface(xx_train(:,z1_index(i)));
    drawnow;
end

input('images sorted over z2');
for i = 1:group_size:n_train
    showface(xx_train(:,z2_index(i)));
    drawnow;
end