train,ytrain,'.b', 'markersize', 50);
    hold on;
    plot(xtest, ypredTest, 'k', 'linewidth', 3, 'markersize', 20);
    hold off
    title(sprintf('degree %d', deg))
    set(gca,'ylim',[-10 15]);
    set(gca,'xlim',[-1 21]);
    printPmtkFigure(sprintf('polyfitDemo%d', deg))
end

%% Compute log evidence for each model
for m=1:length(degs)
    deg = degs(m);
    pp = preprocessorCreate('rescaleX', true, 'poly', deg, 'addOnes', true);
    [modelEB, logev(m)] = linregFitBayes(xtrain, ytrain, 'preproc', pp, .