modelEB, logev(m)] = linregFitBayes(xtrain, ytrain, <span class="string">'preproc'</span>, pp, <span class="keyword">...</span>
        <span class="string">'prior'</span>, <span class="string">'eb'</span>);
<span class="keyword">end</span>

figure;
probs = exp(normalizeLogspace(logev));
plot(degs, logev ,<span class="string">'ko-'</span>, <span class="string">'linewidth'</span>, 2, <span class="string">'markersize'</span>, 24);
xlabel(<span class="string">'degree'</span>); ylabel(<span class="string">'log 