lhlambdaalpha.m

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

function lh=lhlambdaAlpha(y,alpha,pall,ff,pars)

K=length(unique(y));c=K-1;

pref=1-sum(pall,2);
N=length(y);

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))));
normalpha=0;
for cc=1:c,
  normalpha=normalpha+alpha(:,cc)'*ff(:,cc);
end

lh=lh - 0.5*N*pars.lambda*normalpha;