% Copyright (C) 2016 William Craelius (craelius@rci.rutgers.edu)

% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% any later version.

% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <http://www.gnu.org/licenses/>.

function varargout = CHARIS_GUI(varargin)

% CHARIS_GUI Version 1.0 

% CHARIS_GUI MATLAB code for CHARIS_GUI.fig
%      CHARIS_GUI, by itself, creates a new CHARIS_GUI or raises the existing
%      singleton*.
%
%      H = CHARIS_GUI returns the handle to a new CHARIS_GUI or the handle to
%      the existing singleton*.
%
%      CHARIS_GUI('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in CHARIS_GUI.M with the given input arguments.
%
%      CHARIS_GUI('Property','Value',...) creates a new CHARIS_GUI or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before CHARIS_GUI_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to CHARIS_GUI_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help CHARIS_GUI

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
    'gui_Singleton',  gui_Singleton, ...
    'gui_OpeningFcn', @CHARIS_GUI_OpeningFcn, ...
    'gui_OutputFcn',  @CHARIS_GUI_OutputFcn, ...
    'gui_LayoutFcn',  [] , ...
    'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT

% --- Executes just before CHARIS_GUI is made visible.
function CHARIS_GUI_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to CHARIS_GUI (see VARARGIN)

% Choose default command line output for CHARIS_GUI
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% --- Outputs from this function are returned to the command line.
function varargout = CHARIS_GUI_OutputFcn(hObject, eventdata, handles)
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;

% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global ICP
global ABP
global ECG
global Time
global OrigTable
global count
global PatientData
global eventTimes

PatientData=[];
count=0;
ICP=[];
ABP=[];
ECG=[];
Time=[];
eventTimes=[];


set(handles.text42,'String',' ')
% Get PRx Data
edit8=str2num(char(get(handles.edit8,'String')));edit9=str2num(char(get(handles.edit9,'String')));edit10=str2num(char(get(handles.edit10,'String')));
dirOriginal=pwd;

% Load Source Files
set(handles.text42,'String','Data can be loaded from LVMs, source files, or from CHARIS Data.mat');
choice=questdlg('Has data for this patient already been extracted from the LVMs into source files (.mat) or a CHARIS Data.mat?','Find Source Files','Yes','No','Yes');
switch choice
    case 'Yes'
        SaveFolder=1;
    case 'No'
        SaveFolder=0;
    case 'Cancel'
      cd(dirOriginal)
end

if SaveFolder==1
    set(handles.text42,'String','Select folder containing source files (.mat) or the CHARIS Data.mat');
    fprintf('Select folder with the desired source files (.mat)\n ');
    dirRawLocation = uigetdir('Select folder with source files (.mat)\n C:\...\PatientFolder\Source Channels \n');
    cd(dirRawLocation)
    found=0;
    found2=0;
    count1=0;
    count2=0;
    
    % Count files for CHARIS Data
    files = dir('*.mat');
    for i=1:length(files)
        if strncmpi(files(i).name,'CHARIS',5)==1
            count1=count1+1;
        end
    end
    
    % Multiple CHARIS found
    if count1>1
        str=sprintf('CHARIS parts 1-%d detected. Type a number to load desired part.\nCHARIS files are ordered as they are ordered in the folder.',count1);
        whichPart1=inputdlg(str);
        
        partLoad1=str2num(whichPart1{1,1});
        for i=1:length(files)
            if strncmpi(files(i).name,'CHARIS',5)==1
                count2=count2+1;
                if count2==partLoad1
                    load(files(i).name);
                    cd(dirOriginal)
                    found=1;
                end
            end
        end
        
        
    else
        %Only 1 CHARIS found
        for i=1:length(files)
            if strncmpi(files(i).name,'CHARIS',5)==1
                load(files(i).name);
                cd(dirOriginal)
                found=1;
                
            end
        end
        
    end

    if found==0
        for i=1:length(files)
            if strncmpi(files(i).name,'Time Part',9)==1                
               found2=1;
               count2=count2+1;
            end
        end
    end

    % No CHARIS Data found. Source Parts found (.mat)
    if found2==1
        str1=sprintf('Source Parts 1-%d detected. Type a number to load desired part.\nSource Parts are ordered as they are ordered in the folder.',count2);
        whichPart=inputdlg(str1);
        partLoad=str2num(whichPart{1,1});
        for i=1:length(files)
            strT=sprintf('Time Part %d',partLoad);
            strA=sprintf('ABP Part %d',partLoad);
            strI=sprintf('ICP Part %d',partLoad);
            strE=sprintf('ECG Part %d',partLoad);
            
            if strncmpi(files(i).name,strA,10)==1
                ABP=importdata(files(i).name);
            elseif strncmpi(files(i).name,strI,10)==1
                ICP=importdata(files(i).name);
            elseif strncmpi(files(i).name,strE,10)==1
                ECG=importdata(files(i).name);
            elseif strncmpi(files(i).name,strT,11)==1
                Time=importdata(files(i).name);
            end
        end
         cd(dirOriginal)
    end
    
    
    % No CHARIS data found and no source files found
    if found==0&&found2==0
        for i=1:length(files)
            if strncmpi(files(i).name,'ABP',3)==1
                ABP=importdata(files(i).name);
            elseif strncmpi(files(i).name,'ICP',3)==1
                ICP=importdata(files(i).name);
            elseif strncmpi(files(i).name,'ECG',3)==1
                ECG=importdata(files(i).name);
            elseif strncmpi(files(i).name,'Time',4)==1
                Time=importdata(files(i).name);
            end
        end
    end
    set(handles.text42,'String','Finished Loading Source Files')
else
    set(handles.text42,'String','Select folder containing LVM files.');
    dirName = uigetdir('Select folder containing LVM files.');
    set(handles.text42,'String','Select or create a folder for the source files.');
    dirRawSave = uigetdir('Select or create a folder for the source files.');
    set(handles.text42,'String','Compiling Data');
    
    cd(dirName)
    lvmFiles = dir('*.lvm');
    numfiles = length(lvmFiles);
    roundUp=ceil(numfiles/100);
    
    %Make sure that this is the correct column
    tempMatrix = importLVMfile(lvmFiles(1).name, 25, 180024);
    SaveABP=0;
    SaveICP=0;
    SaveECG=0;
    
    while SaveABP==0||SaveICP==0||SaveECG==0
        if SaveABP==0
        
        figure
        plot(tempMatrix(:,edit8));
        set(handles.text42,'String','Examine the figure for ABP data and then close figure to continue.');
        uiwait(gcf)
        
        figure
        plot(tempMatrix(:,edit8));
        choice=questdlg('Is this ABP data?','Find ABP Column','Yes','No','Yes');        
        switch choice
            case 'Yes'
                SaveABP=1;
            case 'No'
                SaveABP=0;
                edit8=inputdlg('Set new ABP Channel (Channels are 2-9)');
                edit9=str2num(edit9{1,1});
                set(handles.edit8,'String',edit8)
            case 'Cancel'
                break
        end
        
        close
        end
        if SaveICP==0
        figure
        plot(tempMatrix(:,edit9));
        set(handles.text42,'String','Examine figure for ICP data and then close figure to continue.');
        uiwait(gcf)
        
        figure     
        plot(tempMatrix(:,edit9));
        choice=questdlg('Is this ICP data?','Find ICP Column','Yes','No','Yes');
        switch choice
            case 'Yes'
                SaveICP=1;
            case 'No'
                SaveICP=0;
                edit9=inputdlg('Set new ICP Channel (Channels are 2-9)');
                edit9=str2num(edit9{1,1});
                set(handles.edit9,'String',edit9)
            case 'Cancel'
                break
        end
        close
        end
        
        if SaveECG==0
        figure
        plot(tempMatrix(:,edit10));
        set(handles.text42,'String','Examine figure for ECG and then close figure to continue.');
        uiwait(gcf)
        
        figure 
        plot(tempMatrix(:,edit10));
        choice=questdlg('Is this ECG data?','Find ECG Column','Yes','No','Yes');
        switch choice
            case 'Yes'
                SaveECG=1;
            case 'No'
                SaveECG=0;
                edit10=inputdlg('Set new ECG Channel (Channels are 2-9)');
                edit10=str2num(edit10{1,1});
                set(handles.edit10,'String',edit10)
            case 'Cancel'
                break
        end
        close
        end
    end
    cd(dirOriginal)
    
    if numfiles>=100
        set(handles.text42,'String','Over 100 files detected. Source files will be split into 100-file parts to improve performance');
        h=waitbar(0,'Please Wait (Files exceed 100)');
        
        for i=1:roundUp
            if i==roundUp;
                indTop=numfiles;
                indBot=roundUp*100-100+1;
            else
                indTop=i*100;
                indBot=indTop-100+1;
            end
            
            [Time1] = CHARISGUIcollectData(indBot,indTop,1,dirName);
            waitbar(((i-1)*4+1)/(roundUp*4))
            [ABP1] = CHARISGUIcollectData(indBot,indTop,edit8,dirName);
            waitbar(((i-1)*4+2)/(roundUp*4))
            [ICP1] = CHARISGUIcollectData(indBot,indTop,edit9,dirName);
            waitbar(((i-1)*4+3)/(roundUp*4))
            [ECG1] = CHARISGUIcollectData(indBot,indTop,edit10,dirName);
            waitbar(((i-1)*4+4)/(roundUp*4))
            
            
            cd(dirRawSave)
            ECG1=cell2mat(ECG1);
            ICP1=cell2mat(ICP1).*100;
            ABP1=cell2mat(ABP1);
            Time1=cell2mat(Time1);
            str1=sprintf('ECG Part %d.mat',i);
            str2=sprintf('ICP Part %d.mat',i);
            str3=sprintf('ABP Part %d.mat',i);
            str4=sprintf('Time Part %d.mat',i);
            save(str1,'ECG1');
            save(str2,'ICP1');
            save(str3,'ABP1');
            save(str4,'Time1');
            
            if i==1
                Time=Time1;
                ABP=ABP1;
                ICP=ICP1;
                ECG=ECG1;
            end
            cd(dirOriginal)
        end
        close(h)
        
        
    else
        h=waitbar(0,'Please wait...');
        [Time] = CHARISGUIcollectData(1,numfiles,1,dirName);
        waitbar(1/4)
        [ABP] = CHARISGUIcollectData(1,numfiles,edit8,dirName);
        waitbar(2/4)
        [ICP] = CHARISGUIcollectData(1,numfiles,edit9,dirName);
        waitbar(3/4)
        [ECG] = CHARISGUIcollectData(1,numfiles,edit10,dirName);
        waitbar(4/4)
        close(h)
        
        cd(dirRawSave)
        ECG=cell2mat(ECG);
        ICP=cell2mat(ICP).*100;
        ABP=cell2mat(ABP);
        Time=cell2mat(Time);
        
        save('ECG.mat','ECG');
        save('ICP.mat','ICP');
        save('ABP.mat','ABP');
        save('Time.mat','Time');
        set(handles.text42,'String','Finished Compiling');
    end
    
end
if isequal(length(Time),length(ABP))==0
    if length(Time)>length(ABP)
        fix=zeros(length(Time)-length(ABP),1);
        ABP=[ABP;fix];
    elseif length(ABP)>length(Time)
        ABP=ABP(1:length(Time));
    end
end
if isequal(length(Time),length(ICP))==0
    if length(Time)>length(ICP)
        fix=zeros(length(Time)-length(ICP),1);
        ICP=[ICP;fix];
    elseif length(ICP)>length(Time)
        ICP=ICP(1:length(Time));
    end
end
if isequal(length(Time),length(ECG))==0
    if length(Time)>length(ECG)
        fix=zeros(length(Time)-length(ECG),1);
        ECG=[ECG;fix];
    elseif length(ECG)>length(Time)
        ECG=ECG(1:length(Time));
    end
end
cd(dirOriginal);
cla(handles.axes1)
cla(handles.axes2)
cla(handles.axes3)
cla(handles.axes4)
set(handles.axes1, 'Xdir', 'normal')
set(handles.axes2, 'Xdir', 'normal')
set(handles.axes3, 'Xdir', 'normal')
set(handles.axes4, 'Xdir', 'normal')
if length(Time)==length(ICP)
    plot(handles.axes1,Time,ICP);
end
set(handles.text1, 'String', 'Raw ICP');
set(handles.text2, 'String', 'Time in seconds');
set(handles.text3, 'String', 'mmHg');
if length(Time)==length(ABP)
    plot(handles.axes2,Time,ABP);
end
set(handles.text4, 'String', 'Raw ABP');
set(handles.text10, 'String', 'Time in seconds');
set(handles.text11, 'String', 'mmHg');
if length(Time)==length(ECG)
    plot(handles.axes3,Time,ECG);
end
set(handles.text12, 'String', 'Raw ECG');
set(handles.text13, 'String', 'Time in seconds');
set(handles.text14, 'String', 'mV');
set(handles.uitable1,'Data',[Time,ABP,ICP,ECG])
OrigTable=[Time,ABP,ICP,ECG];

% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global PatientData
global PRx
global PacketAvg_Time
global ICP
global ABP
global ECG
global Time
global eventTimes
global critStr
global ICPevent_index

ICPevent_index=[];
tableData=get(handles.uitable1,'Data');
edit1 = get(handles.edit1,'String');
edit2=str2num(char(get(handles.edit2,'String')));
edit3=str2num(char(get(handles.edit3,'String')));
edit4=str2num(char(get(handles.edit4,'String')));
edit5=str2num(char(get(handles.edit5,'String')));
set(handles.edit17,'String',edit1);
critStr1=get(handles.edit2,'String');
critStr2=get(handles.edit3,'String');
critStr3=get(handles.edit4,'String');
critStr4=get(handles.edit5,'String');


edit6=edit5*12;
set(handles.text42,'String','Finding Events')
critStr=sprintf('Patient ID: %s \n ICP Threshold: %s mmHg \n Event Length: %s minutes \n Pre-Event: %s minutes \n Examined: %s minutes', edit1,critStr1,critStr2,critStr3,critStr4);

[PatientData,PRxEventTimes,PTM]=CHARISGUIReportPRxEvents(edit1,edit2,edit3,edit4,edit6,tableData(:,3), tableData(:,1),PRx,PacketAvg_Time);
eventTimes=PatientData.eventTimes;
[row,col]=size(eventTimes);
if row==0
    str=sprintf('Events: %d',row);
    set(handles.text42, 'String', str);
    
else
    cla(handles.axes1)
    cla(handles.axes2)
    cla(handles.axes3)
    cla(handles.axes4)
    set(handles.axes1, 'Xdir', 'normal')
    set(handles.axes2, 'Xdir', 'normal')
    set(handles.axes3, 'Xdir', 'normal')
    set(handles.axes4, 'Xdir', 'normal')
    str=sprintf('Events: %d',row);
    set(handles.text42, 'String', str);
    cla(handles.axes1)
    plot(handles.axes1,Time,ICP);
    set(handles.text1, 'String', 'ICP');
    set(handles.text2, 'String', 'Time in seconds');
    set(handles.text3, 'String', 'mmHg');
    
    for i=1:length(eventTimes)
        if isempty(find(Time==eventTimes(i)))==0
            TimeIndex(i,1)=find(Time==eventTimes(i));
        else
            Time=round(Time(:,1)*(10^2))/(10^2);
            TimeIndex(i,1)=find(Time==eventTimes(i));
            
        end
        ICPevent_index(:,i)=ICP(TimeIndex(i,1)-edit5*60/0.02:TimeIndex(i,1));
    end
    
    ICPindex=ICP(TimeIndex);
    
    hold(handles.axes1,'on')
    h=scatter(handles.axes1,eventTimes,ICPindex,'ro');
    set(h,'MarkerEdgeColor','r','MarkerFaceColor','r')
end
cla(handles.axes2)
plot(handles.axes2,Time,ABP);
set(handles.text4, 'String', 'ABP');
set(handles.text10, 'String', 'Time in seconds');
set(handles.text11, 'String', 'mmHg');
if exist('eventTimes')==1
    for i=1:length(eventTimes)
        TimeIndex(i,1)=find(Time==eventTimes(i));
    end
    ABPindex=ABP(TimeIndex);
    hold(handles.axes2,'on')
    h=scatter(handles.axes2,eventTimes,ABPindex,'ro');
    set(h,'MarkerEdgeColor','r','MarkerFaceColor','r')
    hold(handles.axes2,'off')
end
cla(handles.axes3)
plot(handles.axes3,Time,ECG);
set(handles.text12, 'String', 'ECG');
set(handles.text13, 'String', 'Time in seconds');
set(handles.text14, 'String', 'mV');
if exist('eventTimes')==1
    for i=1:length(eventTimes)
        TimeIndex(i,1)=find(Time==eventTimes(i));
    end
    ECGindex=ECG(TimeIndex);
    hold(handles.axes3,'on')
    h=scatter(handles.axes3,eventTimes,ECGindex,'ro');
    set(h,'MarkerEdgeColor','r','MarkerFaceColor','r')
    hold(handles.axes3,'off')
end
set(handles.text42,'String','Finished Finding Events')
cla(handles.axes4)

function edit1_Callback(hObject, eventdata, handles)
% hObject    handle to edit1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
handles.edit1 = get(hObject,'String');

% --- Executes during object creation, after setting all properties.
function edit1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

function edit2_Callback(hObject, eventdata, handles)
% hObject    handle to edit2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
input = str2double(get(hObject,'string'));
if isnan(input)
    errordlg('You must enter a numeric value','Invalid Input','modal')
    uicontrol(hObject)
    return
else
    display(input);
end

% --- Executes during object creation, after setting all properties.
function edit2_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

function edit3_Callback(hObject, eventdata, handles)
% hObject    handle to edit3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
str2double(get(hObject,'String'));

% --- Executes during object creation, after setting all properties.
function edit3_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

function edit4_Callback(hObject, eventdata, handles)
% hObject    handle to edit4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
str2double(get(hObject,'String'))

% --- Executes during object creation, after setting all properties.
function edit4_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

% --- Executes on button press in pushbutton3.
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global ICP
global ABP
global ECG
global Time
global PacketAvg_ICP
global PacketAvg_ABP
global PRx
global count
global eventTimes
global PacketAvg_Time

count=count+1;

tableData=get(handles.uitable1,'Data');

listed=get(handles.listbox2,'String');
selected=get(handles.listbox2,'Value');

graphed=get(handles.listbox3,'Value');
if graphed==1
    graph=handles.axes1;
elseif graphed==2
    graph=handles.axes2;
elseif graphed==3
    graph=handles.axes3;
elseif graphed==4
    graph=handles.axes4;
end

    
if strcmp(listed(selected),'Raw ICP')==1
    
    cla(graph)
    set(graph, 'Xdir', 'normal')
    
    plot(graph,Time,ICP)
    if graphed==1
        %Title
        set(handles.text1, 'String', 'Raw ICP');
        %X values
        set(handles.text2, 'String', 'Time in seconds');
        %Y values
        set(handles.text3, 'String', 'mmHg');
        
    elseif graphed==2
        set(handles.text4, 'String', 'Raw ICP');
        set(handles.text10, 'String', 'Time in seconds');
        set(handles.text11, 'String', 'mmHg');
        
    elseif graphed==3
        set(handles.text12, 'String', 'Raw ICP');
        set(handles.text13, 'String', 'Time in seconds');
        set(handles.text14, 'String', 'mmHg');
        
    elseif graphed==4
        set(handles.text15, 'String', 'Raw ICP');
        set(handles.text16, 'String', 'Time in seconds');
        set(handles.text17, 'String', 'mmHg');
    end
    
    if exist('eventTimes')==1
        for i=1:length(eventTimes)
            TimeIndex(i,1)=find(Time==eventTimes(i));
        end
        ICPindex=ICP(TimeIndex);
        hold(graph,'on')
        h=scatter(graph,eventTimes,ICPindex,'ro');
        set(h,'MarkerEdgeColor','r','MarkerFaceColor','r')
        
    end
    
elseif strcmp(listed(selected),'Raw ABP')==1
    
    cla(graph)
    set(graph, 'Xdir', 'normal')
    
    plot(graph,Time,ABP)
    if graphed==1
        %Title
        set(handles.text1, 'String', 'Raw ABP');
        %X values
        set(handles.text2, 'String', 'Time in seconds');
        %Y values
        set(handles.text3, 'String', 'mmHg');
        
    elseif graphed==2
        set(handles.text4, 'String', 'Raw ABP');
        set(handles.text10, 'String', 'Time in seconds');
        set(handles.text11, 'String', 'mmHg');
        
    elseif graphed==3
        set(handles.text12, 'String', 'Raw ABP');
        set(handles.text13, 'String', 'Time in seconds');
        set(handles.text14, 'String', 'mmHg');
        
    elseif graphed==4
        set(handles.text15, 'String', 'Raw ABP');
        set(handles.text16, 'String', 'Time in seconds');
        set(handles.text17, 'String', 'mmHg');
    end
    
    if exist('eventTimes')==1
        for i=1:length(eventTimes)
            TimeIndex(i,1)=find(Time==eventTimes(i));
        end
        ABPindex=ABP(TimeIndex);
        hold(graph,'on')
        h=scatter(graph,eventTimes,ABPindex,'ro');
        set(h,'MarkerEdgeColor','r','MarkerFaceColor','r')
    end
    
elseif strcmp(listed(selected),'Raw ECG')==1
    
    cla(graph)
    
    set(graph, 'Xdir', 'normal')
    
    plot(graph,Time,ECG)
    if graphed==1
        %Title
        set(handles.text1, 'String', 'Raw ECG');
        %X values
        set(handles.text2, 'String', 'Time in seconds');
        %Y values
        set(handles.text3, 'String', 'mV');
        
    elseif graphed==2
        set(handles.text4, 'String', 'Raw ECG');
        set(handles.text10, 'String', 'Time in seconds');
        set(handles.text11, 'String', 'mV');
        
    elseif graphed==3
        set(handles.text12, 'String', 'Raw ECG');
        set(handles.text13, 'String', 'Time in seconds');
        set(handles.text14, 'String', 'mV');
        
    elseif graphed==4
        set(handles.text15, 'String', 'Raw ECG');
        set(handles.text16, 'String', 'Time in seconds');
        set(handles.text17, 'String', 'mV');
    end
    
    for i=1:length(eventTimes)
        TimeIndex(i,1)=find(Time==eventTimes(i));
    end
    if exist('eventTimes')==1
        ECGindex=ECG(TimeIndex);
        hold(graph,'on')
        h=scatter(graph,eventTimes,ECGindex,'ro');
        set(h,'MarkerEdgeColor','r','MarkerFaceColor','r')
    end
    
    
elseif strcmp(listed(selected),'Packeted ABP')==1
    cla(graph)
    
    set(graph, 'Xdir', 'normal')
    % ABP packeted data plot.
    plot(graph, PacketAvg_ABP.*100,'r.-');

    
    if graphed==1
        
        %Title
        set(handles.text1, 'String', 'Packeted ABP');
        %X values
        set(handles.text2, 'String', 'Seconds');
        %Y values
        set(handles.text3, 'String', '');
        
    elseif graphed==2
        set(handles.text4, 'String', 'Packeted ABP');
        set(handles.text10, 'String', 'Seconds');
        set(handles.text11, 'String', '');
        
    elseif graphed==3
        set(handles.text12, 'String', 'Packeted ABP');
        set(handles.text13, 'String', 'Seconds');
        set(handles.text14, 'String', '');
        
    elseif graphed==4
        set(handles.text15, 'String', 'Packeted ABP');
        set(handles.text16, 'String', 'Seconds');
        set(handles.text17, 'String', '');
    end
    
elseif strcmp(listed(selected),'Packeted ICP')==1
    cla(graph)
    
    set(graph, 'Xdir', 'normal')
    % ICP packeted data plot.
    plot(graph,PacketAvg_ICP.*100,'k.-');

    if graphed==1
        
        %Title
        set(handles.text1, 'String', 'Packeted ICP');
        %X values
        set(handles.text2, 'String', 'Seconds');
        %Y values
        set(handles.text3, 'String', 'mmHg');
        
    elseif graphed==2
        set(handles.text4, 'String', 'Packeted ICP');
        set(handles.text10, 'String', 'Seconds');
        set(handles.text11, 'String', 'mmHg');
        
    elseif graphed==3
        set(handles.text12, 'String', 'Packeted ICP');
        set(handles.text13, 'String', 'Seconds');
        set(handles.text14, 'String', 'mmHg');
        
    elseif graphed==4
        set(handles.text15, 'String', 'Packeted ICP');
        set(handles.text16, 'String', 'Seconds');
        set(handles.text17, 'String', 'mmHg');
    end
    
elseif strcmp(listed(selected),'PRx 5 Minute Packet')==1
    cla(graph)
    
    set(graph, 'Xdir', 'normal')

    plot(graph,PacketAvg_Time,PRx(:,1), 'b.-'); hold on;

    hold off;
    if graphed==1
        
        %Title
        set(handles.text1, 'String', 'PRx 5 Minute Packet');
        %X values
        set(handles.text2, 'String', '');
        %Y values
        set(handles.text3, 'String', '');
        
    elseif graphed==2
        set(handles.text4, 'String', 'PRx 5 Minute Packet');
        set(handles.text10, 'String', '');
        set(handles.text11, 'String', '');
        
    elseif graphed==3
        set(handles.text12, 'String', 'PRx 5 Minute Packet');
        set(handles.text13, 'String', '');
        set(handles.text14, 'String', '');
        
    elseif graphed==4
        set(handles.text15, 'String', 'PRx 5 Minute Packet');
        set(handles.text16, 'String', '');
        set(handles.text17, 'String', '');
    end
elseif strcmp(listed(selected),'PRx 10 Minute Packet')==1
    cla(graph)
    
    set(graph, 'Xdir', 'normal')

    plot(graph,PacketAvg_Time,PRx(:,2), 'b.-'); hold on;

    hold off;
    if graphed==1
        
        %Title
        set(handles.text1, 'String', 'PRx 10 Minute Packet');
        %X values
        set(handles.text2, 'String', '');
        %Y values
        set(handles.text3, 'String', '');
        
    elseif graphed==2
        set(handles.text4, 'String', 'PRx 10 Minute Packet');
        set(handles.text10, 'String', '');
        set(handles.text11, 'String', '');
        
    elseif graphed==3
        set(handles.text12, 'String', 'PRx 10 Minute Packet');
        set(handles.text13, 'String', '');
        set(handles.text14, 'String', '');
        
    elseif graphed==4
        set(handles.text15, 'String', 'PRx 10 Minute Packet');
        set(handles.text16, 'String', '');
        set(handles.text17, 'String', '');
    end
elseif strcmp(listed(selected),'PRx 20 Minute Packet')==1
    cla(graph)
    
    set(graph, 'Xdir', 'normal')

    
    plot(graph,PRx(:,3), 'b.-');  hold on;

    hold off;
    if graphed==1
        
        %Title
        set(handles.text1, 'String', 'PRx 20 Minute Packet');
        %X values
        set(handles.text2, 'String', '');
        %Y values
        set(handles.text3, 'String', '');
        
    elseif graphed==2
        set(handles.text4, 'String', 'PRx 20 Minute Packet');
        set(handles.text10, 'String', '');
        set(handles.text11, 'String', '');
        
    elseif graphed==3
        set(handles.text12, 'String', 'PRx 20 Minute Packet');
        set(handles.text13, 'String', '');
        set(handles.text14, 'String', '');
        
    elseif graphed==4
        set(handles.text15, 'String', 'PRx 20 Minute Packet');
        set(handles.text16, 'String', '');
        set(handles.text17, 'String', '');
    end
    
elseif strcmp(listed(selected),'Event Finder')==1
    cla(graph)
    
    set(graph, 'Xdir', 'normal')

    
    plot(graph,PRx(:,3), 'b.-');  hold on;

    hold off;
    if graphed==1
        
        %Title
        set(handles.text1, 'String', 'PRx 20 Minute Packet');
        %X values
        set(handles.text2, 'String', '');
        %Y values
        set(handles.text3, 'String', '');
        
    elseif graphed==2
        set(handles.text4, 'String', 'PRx 20 Minute Packet');
        set(handles.text10, 'String', '');
        set(handles.text11, 'String', '');
        
    elseif graphed==3
        set(handles.text12, 'String', 'PRx 20 Minute Packet');
        set(handles.text13, 'String', '');
        set(handles.text14, 'String', '');
        
    elseif graphed==4
        set(handles.text15, 'String', 'PRx 20 Minute Packet');
        set(handles.text16, 'String', '');
        set(handles.text17, 'String', '');
    end
end

function edit5_Callback(hObject, eventdata, handles)
% hObject    handle to edit5 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% --- Executes during object creation, after setting all properties.
function edit5_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit5 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

% --- Executes on selection change in listbox2.
function listbox2_Callback(hObject, eventdata, handles)
% hObject    handle to listbox2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% --- Executes during object creation, after setting all properties.
function listbox2_CreateFcn(hObject, eventdata, handles)
% hObject    handle to listbox2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: listbox controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

% --- Executes on selection change in listbox3.
function listbox3_Callback(hObject, eventdata, handles)
% hObject    handle to listbox3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% --- Executes during object creation, after setting all properties.
function listbox3_CreateFcn(hObject, eventdata, handles)
% hObject    handle to listbox3 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: listbox controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

% --- Executes on button press in pushbutton6.
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

cla(handles.axes1)
cla(handles.axes2)
cla(handles.axes3)
cla(handles.axes4)
set(handles.text1, 'String', '');
set(handles.text2, 'String', '');
set(handles.text3, 'String', ' ');
set(handles.text4, 'String', '');
set(handles.text10, 'String', '');
set(handles.text11, 'String', ' ');
set(handles.text12, 'String', '');
set(handles.text13, 'String', '');
set(handles.text14, 'String', ' ');
set(handles.text15, 'String', '');
set(handles.text16, 'String', '');
set(handles.text17, 'String', ' ');

function edit8_Callback(hObject, eventdata, handles)
% hObject    handle to edit8 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% --- Executes during object creation, after setting all properties.
function edit8_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit8 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

function edit9_Callback(hObject, eventdata, handles)
% hObject    handle to edit9 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% --- Executes during object creation, after setting all properties.
function edit9_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit9 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

function edit10_Callback(hObject, eventdata, handles)
% hObject    handle to edit10 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% --- Executes during object creation, after setting all properties.
function edit10_CreateFcn(hObject, eventdata, handles)
% hObject    handle to edit10 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end

% --- Executes on button press in pushbutton8.
function pushbutton8_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton8 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global PRx
global PacketAvg_Time
global PacketAvg_ABP
global PacketAvg_ICP
global ICP
global ABP
global Time

h=waitbar(0,'Please Wait...');

edit8=str2num(char(get(handles.edit8,'String')));edit9=str2num(char(get(handles.edit9,'String')));
SelectedHours=100000000000000000;
set(handles.text42,'String','Please select LVM folder')
[PRx,PacketAvg_Time,PacketAvg_ABP,PacketAvg_ICP]=CHARISGUIpacketAve2(1,SelectedHours,ABP,ICP,Time,1);
waitbar(1/2)
cla(handles.axes1);
cla(handles.axes2);
cla(handles.axes3);
cla(handles.axes4);

% ABP packeted data plot.
plot(handles.axes1,PacketAvg_Time, PacketAvg_ABP.*100,'r.-');

set(handles.text3,'String','Seconds')
set(handles.text1,'String','Packeted ABP')

% ICP packeted data plot.
plot(handles.axes2,PacketAvg_Time,PacketAvg_ICP.*100,'k.-');

set(handles.text4,'String','Packeted ICP')
set(handles.text11,'String','Seconds')

% 5 minutes 
plot(handles.axes3,PacketAvg_Time,PRx(:,1), 'b.-'); hold on;

hold off
set(handles.text12,'String','PRx 5 Minute Packet')
set(handles.text13,'String','')

% 10 minutes 
plot(handles.axes4,PacketAvg_Time,PRx(:,2), 'b.-'); hold on;

hold off
set(handles.text15,'String','PRx 10 Minute Packet');
set(handles.text16,'String','')
set(handles.text42,'String','Finished Packet Averaging All')
waitbar(1)
close(h)

% --- Executes on button press in pushbutton12.
function pushbutton12_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton12 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global PRx
global PacketAvg_Time
global PacketAvg_ABP
global PacketAvg_ICP
global ICP
global ABP
global ECG
global Time
global PatientData
global eventTimes
global critStr
global ICPevent_index

set(handles.text42, 'String', 'Please select a folder to save CHARIS Data.mat');
edit17=get(handles.edit17,'String');
dirNewSave = uigetdir('Select a folder to save CHARIS Data.mat');
olddir=pwd;
cd(dirNewSave)
str1=sprintf('CHARIS Data %s (All Data).mat',edit17);
save(str1,'PRx','PacketAvg_Time','PacketAvg_ABP','PacketAvg_ICP','ICP','ABP','ECG','Time','ICPevent_index','eventTimes');

str3=sprintf('Patient Object %s.mat',edit17);
save(str3,'PatientData');

fileStr=sprintf('Patient %s Criteria.txt',edit17);
fid=fopen(fileStr,'wt');
fprintf(fid,critStr);
fclose(fid);

cd(olddir)
set(handles.text42,'String', 'Finished Saving Data')

function SICUNI62101405051824 = importLVMfile(filename, startRow, endRow)

%% Initialize variables.
delimiter = '\t';
if nargin<=2
    startRow = 25;
    endRow = inf;
end

%% Format string for each line of text:
%   column1: double (%f)
%	column2: double (%f)
%   column3: double (%f)
%	column4: double (%f)
%   column5: double (%f)
%	column6: double (%f)
%   column7: double (%f)
%	column8: double (%f)
%   column9: double (%f)
% For more information, see the TEXTSCAN documentation.
formatSpec = '%f%f%f%f%f%f%f%f%f%*s%*s%*s%*s%[^\n\r]';

%% Open the text file.
fileID = fopen(filename,'r');

%% Read columns of data according to format string.
% This call is based on the structure of the file used to generate this
% code. If an error occurs for a different file, try regenerating the code
% from the Import Tool.
dataArray = textscan(fileID, formatSpec, endRow(1)-startRow(1)+1, 'Delimiter', delimiter, 'EmptyValue' ,NaN,'HeaderLines', startRow(1)-1, 'ReturnOnError', false);
for block=2:length(startRow)
    frewind(fileID);
    dataArrayBlock = textscan(fileID, formatSpec, endRow(block)-startRow(block)+1, 'Delimiter', delimiter, 'EmptyValue' ,NaN,'HeaderLines', startRow(block)-1, 'ReturnOnError', false);
    for col=1:length(dataArray)
        dataArray{col} = [dataArray{col};dataArrayBlock{col}];
    end
end

%% Close the text file.
fclose(fileID);


%% Create output variable
SICUNI62101405051824 = [dataArray{1:end-1}];