% _________________________________________________________________________
% จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ
% PREDICTION OF ACUTE HYPOTENSIVE EPISODES USING NEURAL NETWORK MULTIMODELS
% _________________________________________________________________________
% จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ
%        Copyright (C) 2009
%          : Jorge Henriques ...  <jh@dei.uc.pt>,
%          : Teresa Rocha    ...  <teresa@sun.isec.pt>
%        This software is released under the terms of the GNU
%        General Public License (http://www.gnu.org/copyleft/gpl.html)
% จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ
%  mainShowAHE
%__________________________________________________________________________
% จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ
%  Shoe Results from ..\Results
%__________________________________________________________________________
% จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ

clc; close all
disp('___________________________________________________________________')
disp('จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ')
disp('                                mainShowResults                    ')
disp('___________________________________________________________________')
disp('จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ')
pause(1)

%จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ
FORECAST=60;
mmHG   =60;
pathResult='..\Results\';

sinal = 'BP';
nS    = '1';
option=menu(' -------- DATASET -------- ', ...
    ' .  A ', ' .  B ');
switch option
    case 1; dS='A';
    case 2; dS='B';
end

if dS=='A'
    load ..\Results\episodeAHE_A
    load ..\Results\predictAHE_A
    load ..\Results\positionAHE_A
end
if dS=='B'
    load ..\Results\episodeAHE_B
    load ..\Results\predictAHE_B
    load ..\Results\positionAHE_B
end


disp('บบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบ')
disp(['  RESULTS dataSet : '  char(dS) ]);
disp('บบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบบ')
disp(' ')
disp(['      AHE : ' num2str(episodeAHE) ])
disp([' positions: ' num2str(positionAHE(episodeAHE)) ])
disp(' ')
pause

%--------------------------------------------------
DIM=size(predictAHE,1);
for iD=1:size(predictAHE,2)
    ABP=predictAHE(:,iD);
    indm=positionAHE(iD);

    mmMaximize
    id0=DIM-FORECAST;
    figure(1)
    plot(1:id0,ABP(1:id0),'g','LineWidth',10)
    hold on
    plot(1:DIM,ABP,'r','LineWidth',3)
    minV=min(mmHG-10,min(ABP));
    maxV=max(ABP);
    plot(1:DIM,mmHG*ones(DIM,1),'k','LineWidth',6)
    plot([id0 id0],[minV maxV ],'k:','LineWidth',3 );
    plot([id0+FORECAST-1 id0+FORECAST-1],...
        [minV maxV ],'k:','LineWidth',3 );
    if indm<FORECAST
        plot(indm+id0,ABP(indm+id0),'bs','MarkerSize',9,'LineWidth',8  );
        plot(indm+id0,ABP(indm+id0),'y+','MarkerSize',9,'LineWidth',5  );
    end

    axis([0 DIM  minV maxV  ])
    legend('ABP Signal','Model+Prediction','Limit','Orientation','horizontal') 
    strTitle=[  char(dS)  '-'  num2str(iD) ];
    mmLabel(strTitle,'minutes','mmHG')

    %จจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจจ
    pause
    close all

end
option=1;