lssvroutput.m

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

function [Yhat,ROCarea,se] = lssvroutput(alpha,bias,sig2,trnX,trnY, tstX,tstY)
% Calculate lssSVR Output
%  [Yhat,ROCarea,se] = lssvroutput(alfa,bias,sig2,trnX,tstX,tstY)
%  Parameters: trnX   - Training inputs
%              tstX   - Test inputs
%              alfa   - Difference of Lagrange Multipliers
%              bias   - bias
%					sig2 = sig^2 for RBF kernel

  if (nargin < 6)|(nargout>1&nargin<7) % check correct number of arguments
    help svroutput,
else
    n = size(trnX,1);
    m = size(tstX,1);
    H = zeros(m,n);  
    for i=1:m
       for j=1:n
        u=tstX(i,:); v=trnX(j,:);  
%        K(i,j) = exp(-sqrt((u-v)*(u-v)')/(2*sig2)) 
%        K(i,j)=exp(-(norm(tstX(i,:)-trnX(j,:))^2)/sig2);
        K(i,j)=exp(-(norm(u-v)^2)/sig2);       
    end
  end
 
   Yhat = (K*alpha +bias);

   if nargout>1,
   
   Y0=[];
   Y1=[];
   for i=1:length(tstY)
     if sign(tstY(i))<=0,
       Y0=[Y0;Yhat(i)];
     else
       Y1=[Y1;Yhat(i)];
     end
   end
   
   [TN,TP,FP,FN,ROCarea,se]=roc_tf(Y0,Y1);  
end,
end,