redictive'</span>)
printPmtkFigure(<span class="string">'gibbsGaussParamsPostPred'</span>)

<span class="keyword">end</span>

<span class="keyword">function</span> [samples] = gibbsSampler(data, params, nSamples)
<span class="comment">% Sample mu and Sigma given a fully observed MVN</span>
<span class="comment">% samples.mu(:,s), samples.Sigma(:,:,s)</span>
<span class="comment">% We assume a factored Normal * IW prior</span>
<span class="comment">% See Peter Hoff's book sec 7.4</span>
<span class="comment"