%Written by Martin Macas 2012
%Function needs statistical toolbox
function [risk,prediction]=physionet2012(time,param,value)

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);

%Computation of posterior probabilities using Bayes formula
[Pw1x, Pw2x]=computePoster(x);

%Correction of posterior probabilities
if Pw2x>1,Pw2x=1;end
if Pw2x<0,Pw2x=0;end

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

%Prediction is based on MAP principle
prediction=Pw2x>Pw1x;