predictprobs.m

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

function [pt,lhte]=predictprobs(X,y,alpha,Xtest,ytest,pars)

X=[ones(size(X,1),1) X];
Xtest=[ones(size(Xtest,1),1) Xtest];
K=length(unique(y));c=K-1;
[N,d]=size(X);

alphamat=reshape(alpha,c,N)'; % matrix notation
%p=alphaprobs(X,y,alphamat,pars,X);
pt=alphaprobs(Xtest,ytest,alphamat,pars,X);

%compute likelihood
%lhtr=lh(y,p,pars);
lhte=lh(ytest,pt,pars);

%p=[p ones(size(p,1),1)-sum(p,2)];
pt=[pt ones(size(pt,1),1)-sum(pt,2)];