lh.m

The goal of SPID is to provide the user with tools capable to simulate, preprocess, process and clas

function lh=lh(y,pall,pars)

c=size(pall,2);
pref=1-sum(pall,2);

lh=0;
for i=1:c
    index=find(y==i);
    lh=lh + sum(log(max(eps,pall(index,i))));
end
index=find(y==c+1);
lh=lh + sum(log(max(eps,pref(index))));
%normalise
lh=(1/length(y))*lh;