HealthKit_data_plot.m

close all;
clear;

addpath('./utilities')
addpath('./utilities/isValidDate')

if ispc  % Use smaller font on Windows
    font_size = 11;
else
    font_size = 14;
end

fprintf('\nWelcome.\n');

% input_file_name = 'output.csv';

% Get data file path
if ~exist('input_file_name', 'var')
    fprintf('Please select a file (with a ''.csv'' extension) to process... ');
    [input_file_name, input_file_path] = uigetfile('*.csv');
    if input_file_name(1) == 0 || input_file_path(1) == 0
        fprintf('\nNo file selected.\nExiting.\n\n')
        return;
    end
    input_file_name = [input_file_path, input_file_name];
end

% Check if file exists, if not give error message and exit
if exist(input_file_name, 'file') ~= 2
    err_str = sprintf('ERROR: Specified input file ''%s'' could not be found.\nExiting.', input_file_name);
    msgbox(err_str)
    fprintf('\n%s\n\n', err_str);    
    return;
end

fprintf('Done.\nOpening file and determining record types... ');

% Read in data
data = readtable(input_file_name, 'Delimiter', ',', 'ReadVariableNames', true, 'HeaderLines', 0);
headers = data.Properties.VariableNames;
record_type_header = char(data.Properties.VariableNames(1));
date_header = char(data.Properties.VariableNames(2));
value_header = char(data.Properties.VariableNames(3));
unit_header = char(data.Properties.VariableNames(4));

% Create record and unit table
[record_types, ia, ~] = unique(data.(record_type_header)); % Find unique record types
record_types(:, 2) = data.(unit_header)(ia); % Find units associated with record types
record_types = cell2table(record_types, 'VariableNames', {record_type_header, unit_header});

% Have user select which record types to keep
fprintf('Done.\nPlease select record types to plot and plot type...  ');
select_record_types(record_types);
  % Returns check box booleans as 'selected_record_types'
uiwait(user_input_figure);  % Wait until user input window has been closed

% Exit if no record types were selected
if ~exist('selected_record_types', 'var') || sum(selected_record_types) == 0
    fprintf('\nNo record types selected.\nExiting.\n\n')
    return;
end

fprintf('Done.\nProcessing data... ');

% Keep only selected records
record_types = record_types(selected_record_types, :);
data = data(ismember(data.(record_type_header), record_types.(record_type_header)), :);

% Get dates as strings, remove UTC offsets, convert to numeric values
dates = data.(date_header);
regex = ' [-+±]\d\d\d\d';
dates = regexprep(dates, regex, ''); % Remove UTC offsets and convert to char
dates = datenum(dates);
data.(date_header) = dates;
clear dates;

% Store tables for each record type in individual cells in cell array
num_record_types = size(record_types, 1);
data_to_plot = cell(num_record_types, 1);
for record_type_ind = 1:num_record_types
    data_indices_bool = ismember(data.(record_type_header), record_types.(record_type_header)(record_type_ind));
    data_to_plot{record_type_ind} = data(data_indices_bool, :);
end

fprintf('Done.\n');

% Get initial min and max dates based on user selected record types
min_date_init = min(data_to_plot{:}.(date_header));
max_date_init = max(data_to_plot{:}.(date_header));

% Initial plot parameters
func_name = 'scatter';
x_label_interval = 'Auto';
aggregate_data_by = 'none';

plot_data