tt> to allow automatic sampling and estimation of the error variance. The option <tt>qcov</tt> sets the initial covariance for the Gaussian proposal density of the MCMC sampler.</p><pre class="codeinput">options.nsimu = 4000;
options.updatesigma = 1;
options.qcov = tcov; <span class="comment">% covariance from the initial fit</span>
</pre><p>The actual MCMC simulation run is done using the function <tt>mcmcrun</tt>.</p><pre class="codeinput">[res,chain,s2chain] = mcmcrun(model,data,params,options);
</pre><p