% create new VMP
vmp = xff('new:vmp');

% generate random data
vmp.Map.VMPData = single(15.*smoothdata3(randn(size(vmp.Map.VMPData)), [2,2,2]));

% mask with Colin's brain (we don't want stuff that couldn't be good)
cb = neuroelf_file('c', 'colin_brain.vmr');
vmp.MaskWithVMR(cb, 11);

% produce cluster table outputs
[c,t,v,vo] = vmp.ClusterTable(1, -sdist('tinv', 1e-6, vmp.Map.DF1), struct( ...
    'minsize', 40, 'mni2tal', true, 'tdclient', true));