ly selected subset of the chain and calculate the predictive envelope of the model. The grey areas in the plot correspond to 50%, 90%, 95%, and 99% posterior regions.</p><pre class="codeinput">figure(5); clf
out = mcmcpred(res,chain,[],x,modelfun);
mcmcpredplot(out);
hold <span class="string">on</span>
plot(data.xdata,data.ydata,<span class="string">'s'</span>); <span class="comment">% add data points to the plot</span>
xlabel(<span class="string">'x [mg/L COD]'</span>); ylabel(<span class="string">'y [1/h]