loglikelihood.m.txt

这是马尔可夫-蒙特卡罗算法的MATLAB源程序.

function llv = logLikelihood(p,d)
% llv = logLikelihood(p,d)
% llv = logLikelihood(G,d)
% llv = logLikelihood(Data,d)
% Calculate log likelihood for state p, Green's function G, or data Data, given data structure d

% Colin Fox, 21 Jan 2003.

if isstruct(p)
    pData = p.G(d.ElecNodes,:);
else
    if all(size(p) == size(d.Data))
        % p is mock data
        pData = p;
    else
        % p is a set of Greens functions
        pData = p(d.ElecNodes,:);
    end
end

llv = -sum(sum((pData-d.Data).^2))/(2*d.StdDev.^2);