str(c) <span class="string">'   k = '</span> num2str(kk) <span class="string">'  :  posterior mass on all C(n,k) modes'</span>],<span class="string">'fontsize'</span>,14)
       xlabel(<span class="string">'mode #'</span>), axis([0.5 Ns+.5 0 1]); drawnow

       <span class="comment">% what if you had equal-weighting?</span>
       <span class="comment">%Pmode = ones(Ns,1)/Ns;   % but this is very unrealistic!</span>

       <span class="comment">% accumulate marginal histogram (x_i are exchangeable)</span>