function [likelihood] = evaluateSwatch(basis,imagename,parameters,mndata,method)

sz = size(basis);

region = getSwatch(imagename);
region = region(1:sz(1));

if(mndata(1,1) ~= 0)
66
region = region - mndata';
end

projection = region'*basis;

likelihood = marginalProbability(projection,parameters,method);