>%</span>
<span class="comment">% Written by Emtiyaz, CS, UBC</span>
<span class="comment">% Modified on Jan 29, 2010</span>


<span class="comment">% hyperparameters</span>
mean0 = params.mean0;
precMat0 = inv(params.covMat0);
nu0 = params.nu0;
S0 = params.S0;
[d,n] = size(data);

<span class="comment">% Initialize chain by setting mu= sample mean</span>
xbar = mean(data,2);
mu = xbar;
<span class="keyword">for</span> i = 1:nSamples
  <span class="comment">% sample Sigma given mu</span>
  nu = nu0 + n;
  d