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, 200);
xlabel('log likelihood on training set');
input('log likelihood on training set');

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

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

[junk, i] = sort(z_train(2,:));
subplot(4, 4, 3);
showface(xx_train(:,i(n_train)));
subplot(4, 4, 7);
showface(xx_train(:,i(n_train)));
subplot(4, 4, 15);
showface(xx_train(:,i(1)));

subplot(4, 4, 12);
[junk, i] = max(z_train(1,:));
showface(xx_train(:,i));
subplot(4, 4, 9);
[junk, i] = min(z_train(1,:));
showface(xx_train(:,i));