function likelihood = evaluateSample(basis,imageVector,parameters,mndata,method)


sz = size(basis);

region = imageVector(1:sz(1));

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

projection = region'*basis;

likelihood = marginalProbability(projection,parameters,method);