', 'fontsize', fs)
xlabel('w^{MLE}', 'fontsize', fs)
ylabel('w^{MAP}', 'fontsize', fs)
printPmtkFigure('NGthreshold')
%}
%%
%{
    figure();
    hold on
    for i=1:length(deltas)
        plot(z,x(outNIG{i}),styles{i}, 'linewidth', 3);
        str{i} = sprintf('%s = %5.3f, c=1', '\delta', deltas(i));
    end
    legend(str, 'location', 'southeast')
    plot(z,z,'r:', 'linewidth', 3)
    title('normalInvGauss')
    printPmtkFigure('NIGthreshold')
end
%}
%%
%{
figure
hold on
for i=1:length(outNEG)
    plot(z,