function RemovedWindowsStats(Windows, AFWindows, HRVparams, sub_id)
% RemovedWindowsStats(Windows, AFWindows, HRVparams, sub_id)
% OVERVIEW:
% Returns a file with the total %ages of windows removed, the %ages
% removed due to AF and %ages low quality signal
%
% INPUT:
% Windows : vector containing windows indexes corresponding
% to RR interval segmentation for HRV analysis
% AFWindows : vector containing windows indexes corresponding
% to RR interval segmentation for AF analysiss
% HRVparams : struct containing the parameters for HRV analysis
% subjectID : string with subject name
%
% DEPENDENCIES & LIBRARIES:
% PhysioNet Cardiovascular Signal Toolbox
% https://github.com/cliffordlab/PhysioNet-Cardiovascular-Signal-Toolbox
%
% REFERENCE:
% Vest et al. "An Open Source Benchmarked HRV Toolbox for Cardiovascular
% Waveform and Interval Analysis" Physiological Measurement (In Press), 2018.
%
% REPO:
% https://github.com/cliffordlab/PhysioNet-Cardiovascular-Signal-Toolbox% ORIGINAL SOURCE AND AUTHORS:
% Written by Giulia Da Poian (giulia.dap@gmail.com) on Sep 6, 2017.
% Dependent scripts written by various authors
% (see functions for details)
% COPYRIGHT (C) 2018
% LICENSE:
% This software is offered freely and without warranty under
% the GNU (v3 or later) public license. See license file for
% more information
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TotWind = length(Windows);
WindNotAnalyzed = numel(find(isnan(Windows)));
AFWind = numel(Windows(AFWindows));
LowQualityWind = WindNotAnalyzed - AFWind;
percentAF = AFWind/TotWind *100;
percentLowQuality = LowQualityWind/TotWind *100;
percentNotAnalyzed = WindNotAnalyzed/TotWind*100;
filename = strcat(HRVparams.writedata, filesep,'RemovedWindowsStatsSummary_', HRVparams.time,'.csv');
try
T = readtable(filename);
catch
T = []; % If file does not exist yet
end
variables_names = [{'patID' 'TotWind' 'NotAnalyzed' 'PercentNotAnalyzed' 'PercentAFWind' 'PercentLowQualityWind'}]; % Add colum with patient ID
patid_array = string({sub_id});
variables_vals = [patid_array TotWind WindNotAnalyzed percentNotAnalyzed percentAF percentLowQuality]; % Add colum with patient ID
% Write AF results to a table
T = [T ; array2table(variables_vals,'VariableNames',variables_names)];
% Use writetable to geberate csv file with the results
writetable(T,filename);