uce the example on p64 of Marin and Robert, 2007 This computes logBF(i) = log p(D|X(:,:)) / p(D|X(:,-j)) for each feature j using a g-prior</p><!--/introduction--><pre class="codeinput"><span class="comment">% This file is from pmtk3.googlecode.com</span>

X = loadData(<span class="string">'caterpillar'</span>);
<span class="comment">% from http://www.ceremade.dauphine.fr/~xian/BCS/caterpillar</span>
y = log(X(:,11)); <span class="comment">% log number of nests</span>
X = X(:,1:10);
[n,d] = size(X);
X1 = [o