function patient_data = load_patient_data_from_filename(filename,patient_data,total_time)
% This function takes raw patient data from spreadsheet and returns 
% only useful data

    % only take relevant data (time, pain, LA drug, SA drug, Non-opiod drug)
    patient_data = [patient_data(:,1),patient_data(:,10),patient_data(:,30),patient_data(:,31),patient_data(:,32)];
    % for patients with initial weird report, skip first report
    if (strcmp(filename,'PatientA0004') || strcmp(filename,'PatientA0009')) || (strcmp(filename,'PatientB1009') || strcmp(filename,'PatientB1010'))
        % FOR PATIENTA0004/9, PATIENTB1009/10 ONLY!!! Delete first report time (outlier)
        patient_data(:,1) = patient_data(:,1) - patient_data(2,1);
        patient_data(1,:) = [];
    end
    % only look at small period of time (parameters may shift over time)
    patient_data = patient_data(patient_data(:,1)<total_time,:);

end