LogLik(muSigma,measurements,resolution)

function logLik = LogLik(muSigma,measurements,resolution)
% LogLik evaluates the log-likelihood for discrete (digitized)
% measurements

%Viktor Witkovsky ([e-mail])
%Revised: 20-May-2008 11:15:52

measurements = sort(measurements(:));
minMeasurements = min(measurements);
maxMeasurements = max(measurements);
scale =(minMeasurements - 2 * resolution : resolution : maxMeasurements)';
edges = scale + resolution/2;
counts = histc(measurements,edges);
scale = scale + resolution;
nScale = numel(scale);

rounding = 0.5;
upp = (1 - rounding) * resolution;
low = -rounding * resolution;
sampleSize = size(muSigma,1);
logLik = zeros(sampleSize,1);
for i = 1:nScale
    if counts(i) ~=  0
        cdif = normcdf((scale(i) + upp - muSigma(:,1)) ./ (muSigma(:,2)+eps)) ...
            - normcdf((scale(i) + low - muSigma(:,1)) ./ (muSigma(:,2)+eps));
        logLik = logLik + counts(i) * log(cdif + eps);
    end
end