%Written by Martin Macas 2012
%Implementation tested for MATLAB 2011b 7.13.0.564
%needs Statistics toolbox
function [risk,prediction]=physionet2012(time,param,value)
Pw2=0.28;%POZOR
beta=0.625;%risk correction

if (~iscell(value))
value=cellstr(num2str(value));%This works for genresults,
%but does not correspond to 
%physionet2012.m, where "value" should be cell array. 
end

%Feature extraction uses extraction scripts f_******.mat
x=extractFeatures(time,param,value);
load load_classifier;
[y,p]=predict(cl, x);

Pw2x=p(:,2);

%Risk is proportional to aposterior probability
risk=beta*Pw2x;%POZOR


%Prediction is based on MAP principle
prediction=y;