Monitoring physical activity with consumers wearables is one of the possibilities to control a patient’s self-care and adherence to recommendations. However, clinically approved methods, software, and data analysis technologies to collect data and make it suitable for practical use for patient care are still lacking.
This study aimed to analyze the potential of patient physical activity monitoring using Fitbit physical activity trackers and find solutions for possible implementation in the health care routine.
Thirty patients with impaired fasting glycemia were randomly selected and participated for 6 months. Physical activity variability was evaluated and parameters were calculated using data from Fitbit Inspire devices.
Changes in parameters were found and correlation between clinical data (HbA1c, lipids) and physical activity variability were assessed. Better correlation with variability than with body composition changes shows the potential to include nonlinear variability parameters analysing physical activity using mobile devices. Less expressed variability shows better relationship with control of prediabetic and lipid parameters.
Evaluation of physical activity variability is essential for patient health, and these methods used to calculate it is an effective way to analyze big data from wearable devices in future trials.
Lack of physical activity and obesity belong to the biggest health problems we face. This leads to increased mortality related to cardiovascular diseases and diabetes. One of the conditions caused by lack of physical activity and malnutrition is prediabetes. Prediabetes is a transitional hyperglycaemic condition when glucose levels are elevated but not high enough to meet the diagnostic criteria for diabetes [1, 2].
However, it is a reversible condition, and prompt diagnosis and preventive action can stop or delay development or progression of its complications . One of the actions to prevent these chronic diseases is increasing and promoting of physical activity.
The World Health Organization (WHO) 2010 guidelines state that an adult without significant mobility impairments should assign at least 150 minutes per week to moderate-intensity aerobic exercise or at least 75 minutes to high-intensity physical exercise or an equivalent combination of physical activity. To improve health and prevent chronic diseases, it is recommended at least to double time indicated in the basic guidelines for physical activity [4, 5, 6].
Healthy lifestyle professionals and disease prevention program developers are actively looking for ways to adapt and introduce more active lifestyle habits to the general public. Walking is considered a form of moderate physical activity and is a great way to maintain the required amount of physical activity [7, 8, 9, 10]. It not only facilitates easy access from point A to point B but also contributes to better health. The American College of Sports Medicine has systematized WHO recommendations and compared the usual recommended physical activity with 7,000–9,000 steps per day [11, 12, 13]. Recent years went through enormous mobile devices breakthroughs, and now we can count our steps using smartphones, smartwatches, and activity trackers  Smart bracelets measure the daily number of steps with sufficient accuracy and are suitable for long-term physical activity assessment . Analysis of the data collected by these smart devices has shown that, on average, people walk 4,000 to 6,000 steps per day, which is only two-thirds of the recommended physical activity .
Now we have a simple way to test the physical activity change in free-living conditions using consumer wearables, such as smartwatches and physical activity trackers. Patients have the opportunity to perform self-monitoring of their physical activity. Some studies show us how much time people spend sedentary or physically active. Studies also show that physical activity during the week varies and can differ significantly [17, 18, 19, 20].
The aim of this study was to evaluate how physical activity variability, in the long run, impacts people’s health in free-living conditions and to analyse what impact physical activity variability in 6 months has on their glycated hemoglobin, cholesterol, and body composition. Our study included people with prediabetes who have a higher risk of developing diabetes.
During the literature review, we haven’t found the same or similar research looking for physical activity variability impact to people’s health in free-living conditions, so we have got some intriguing results for future research.
The aims of this study are: 1) to analyze the potential of patient physical activity monitoring using Fitbit devices; 2) to find variables associated with lipid profile and HbA1c changes in 6 months; and 3) to find solutions for possible implementation in the health care routine and health information system.
Vilnius University Hospital Santaros Klinikos (VUHSK) Family Medicine Centre patients visiting their family doctor were included in our study. The selected study population had increased fasting glucose levels (5.6 to 6.9 mmol/l), did not have skeletal or muscular diseases, and conditions that compromise patients’ ability to move. The participants had to have a smartphone with the latest operating system, such as Apple iOS 13 or later and Android OS 8.0 or later, an email address and be able to use it.
30 subjects were selected. Subjects were informed about the study, its possibilities and benefits for their health. Informed consent form has been signed.
Blood tests were performed during the first consultation: total cholesterol, triglycerides, low-density lipoproteins, high-density lipoproteins, glycated hemoglobin. We also performed body mass composition analysis with the medical bioimpedance device X-Contact 356. Also, we provided physical activity recommendations and consulted on a healthy diet.
Subjects were given Fitbit Inspire physical activity trackers (PAT) and they were asked to connect the bracelet to their smartphones through the Fitbit application. If the patient failed to perform these actions, the PAT was connected by a family doctor.
After 6 months, during a consultation at the family medicine center, the subject was asked to transfer data from his online Fitbit account about his physical activity during the study period to the principal researcher’s main computer.
During this visit, we repeated blood tests: total cholesterol, triglycerides, low-density lipoproteins, high-density lipoproteins, glycated hemoglobin, and bodyweight analysis. The collected information was analysed and discussed with patients in order to suggest lifestyle changes.
The study was approved by the Vilnius Regional Biomedical Research Ethics Committee (approval no. 2019/6-1143-634).
2.2Physical activity monitoring using Fitbit
Participants received Fitbit Inspire PAT (free of charge) during the consultation with the family doctor and agreed to allow researchers to collect data from their accounts. Subjects installed the Fitbit app on their smartphones and created their own user accounts with a personal email address. Then using a “Bluetooth” connection on the smartphones they connected their smart bracelets and synchronized with their personal Fitbit account. If the participant had difficulty with creating a personal account or connecting a bracelet, we helped him to do this and we taught him how to use the device . Subjects had to periodically synchronise data with their phone throughout the study period by activating the Fitbit application on their smartphone and regularly charge their bracelet every 6–7 days.
After 6 months, during the routine visit at the family doctor’s office, we asked patients to log in to their accounts at https://accounts.fitbit.com/login. Then we used the Fitbit data export platform to download the data covering 30 days interval from start to finish of the study (Figs 1 and 2).
Data obtained from the smart bracelets were: physical activity recording data, steps per day, distance, minutes sedentary, minutes weakly active, active, minutes fairly active, minutes very active.
2.3Poincaré plot and physical activity variability parameters
In this study, the Poincaré plot was used to analyze physical activity variability. The Poincaré plot is a two-dimensional scatterplot of paired values of time series data points, where the x-axis represents the current value at data point and the y-axis represent value at data point with a time delay of (Fig. 3).
This method is widely used to study physiological signals , especially in the analysis of heart rate variability [23, 24, 25]. but there are also promising applications in other areas such as measuring hand grip strength , assessing blood pressure variability , monitoring temporal blood glucose fluctuations [28, 29, 30], studying reaction time during repeated testing , etc.
To analyze physical activity variability, the Poincaré plot method with a time delay of one day () was applied and additional measures for selected physical activity time series accessed from patients’ Fitbit data were calculated. An average value (AVG), standard deviation (SD) and some standard Poincaré plot parameters representing short-term variability (SD1), long-term variability (SD2), the ratio of SD1 and SD2 (SD12) and area of fitting ellipse (AFE) were selected to analyze variability, as defined in other studies [22, 23, 24, 25, 26, 27, 28, 29, 30, 31]:
The Shapiro-Wilk test was used for normality testing. Data was presented as mean and standard deviation (mean SD) or median and quartiles (median [Q1–Q3]) when appropriate. Paired test or Wilcoxon’s matched pairs signed rank test was used for baseline and follow-up data comparison. Pearson’s or Spearman’s correlation was performed to assess association between changes of selected variables and physical activity variability parameters. Analyses were conducted using jamovi 1.6  and R 4.1.0 , and R packages RHRV , corrplot . The list of abbreviations used in this paper is presented in Table 1.
|HgbA1c %||Glycated hemoglobin, percentage||%|
|HgbA1c mmol/l||Glycated hemoglobin, millimoles per liter||mmol/l|
|Bio-electrical impedance analysis|
|BMI||Body Mass Index|
|LBM||Lean body mass||kg|
|MBF||Mass of body fat||kg|
|SLM||Soft lean mass||kg|
|SMM||Skeletal muscle mass||kg|
|TBW||Total body water||kg|
|PBF||Percent body fat||%|
|VFL||Visceral fat level||units|
|VFA||Visceral fat area||cm2|
|WHR||Waist to hip ratio||Ratio|
|SLM Lt. ARM||Left arm soft lean mass||kg|
|SLM Rt.ARM||Right arm soft lean mass||kg|
|SLM Lt.LEG||Left leg soft lean mass||kg|
|SLM Rt.LEG||Right leg soft lean mass||kg|
|SLM Trunk||Trunk soft lean mass||kg|
|MBF Lt. ARM||Left arm mass of body fat||kg|
|MBF Rt.ARM||Right arm mass of body fat||kg|
|MBF Lt.LEG||Left leg mass of body fat||kg|
|MBF Rt.LEG||Right leg mass of body fat||kg|
|MBF Trunk||Trunk mass of body fat||kg|
|E.C.W/T.B.W||Extra-cellular water to total body water ratio||Ratio|
|Fitbit activity tracker data|
|Data||Day of physical activity registration||DD-MM-YYYY|
|S||Steps per day|
|D||Distance per day||km|
|MLA||Minutes lightly active||minutes|
|MFA||Minutes fairly active||minutes|
|MVA||Minutes very active||minutes|
Mean age of study subjects was 53.8 9.1 (median [Q1–Q3]: 54 [49.5–61]), there were 9 males and 21 females. After intervention, significant differences in patient weight, body composition and self-reported physical activity were observed (Table 2).
|Baseline ( 30)||Follow-up ( 30)||value||Effect size|
|HgbA1c [%]#||5.61 0.352||5.65 0.395||0.||367||0.||1674|
|HgbA1c [mmol/l]#||37.7 3.74||38.2 4.36||0.||331||0.||1804|
|TC [mmol/l]#||5.70 1.14||5.90 1.25||0.||338||0.||1779|
|TG [mmol/l]$||1.81 [1.19–2.50]||1.16 [1.16–2.14]||0.||271||0.||2344|
|HDL [mmol/l]$||1.29 [1.12–1.52]||1.11 [1.11–1.61]||0.||593||0.||1140|
|LDL [mmol/l]#||3.44 1.04||3.67 1.23||0.||247||0.||2156|
|Weight [kg]#||87.6 16.9||86.1 16.4||0.||022||0.||4406|
|BMI [kg/m2]$||32.0 [26.4–34.6]||27.0 [27.0–35.0]||0.||054||0.||4115|
|LBM [kg]#||54.2 12.1||56.0 11.1||0.||181||0.||2505|
|MBF [kg]#||32.0 9.72||30.2 9.05||0.||001||0.||7617|
|SLM [kg]#||50.7 10.1||51.2 10.3||0.||080||0.||3312|
|SMM [kg]#||30.4 6.07||30.7 6.21||0.||076||0.||3357|
|Mineral [kg]#||4.93 0.937||4.86 0.926||0.||055||0.||3651|
|Protein [kg]$||10.0 [8.67–12.0]||8.90 [8.90–12.4]||0.||026||0.||4759|
|TBW [kg]#||40.0 7.87||40.4 8.03||0.||141||0.||2764|
|PBF [%]#||36.1 6.86||34.7 6.58||0.||001||0.||8971|
|VFL [Units]$||16.0 [13.0–18.0]||12.3 [12.3–17.0]||0.||001||0.||87619|
|VFA [cm2]#||169 79.1||148 69.1||0.||001||0.||9162|
|WHR#||0.948 0.0912||0.926 0.0864||0.||001||0.||7798|
|AC [cm]#||98.1 13.0||95.9 12.2||0.||001||0.||7618|
|SLM Lt. ARM [kg]#||3.50 0.718||3.51 0.728||0.||751||0.||0584|
|SLM Rt.ARM [kg]#||3.52 0.736||3.53 0.745||0.||795||0.||0478|
|SLM Lt.LEG [kg]#||9.07 2.17||9.41 1.97||0.||095||0.||3155|
|SLM Rt.LEG [kg]#||9.32 1.90||9.51 1.98||0.||008||0.||5158|
|SLM Trunk [kg]#||25.1 4.96||25.2 5.02||0.||258||0.||2105|
|MBF Lt. ARM [kg]#||1.93 0.612||1.86 0.586||0.||008||0.||5227|
|MBF Rt. ARM [kg]#||1.91 0.625||1.82 0.603||0.||005||0.||5518|
|MBF Lt. LEG [kg]#||5.84 1.75||5.50 1.62||0.||001||0.||7810|
|MBF Rt. LEG [kg]#||5.83 1.76||5.47 1.62||0.||001||0.||8146|
|MBF Trunk [kg]#||16.5 4.98||15.5 4.65||0.||001||0.||7666|
|E.C.W/T.B.W#||0.401 0.008||0.399 0.00839||0.||052||0.||3700|
Values are presented as mean SD or median [Q1–Q3]. Two-sided value 0.05. Paired -test, Cohen’s d (effect size). Wilcoxon’s matched pairs signed rank test, Rank biserial correlation (effect size).
3.2Physical activity monitoring using Fitbit
In order to analyze physical activity variability, the Poincaré plots (Fig. 4) were represented and corresponding parameters were calculated for the exported Fitbit data (Table 3). Different patterns of physical activity variability were observed in different patients and different Fitbit data variables. SD1/SD2 ratio varied between 0.697 and 0.870 indicating higher long-term than short-term variability.
|Fitbit variable||Parameter||Mean SD||Median [Q1–Q3]||value|
|Steps per day||SD1||3457||1120||3375||[2601–4157]||0.||631|
|Minutes lightly active [min/day]||SD1||70.0||29.1||59.6||[50.1–75.0]||0.||001|
|Minutes fairly active [min/day]||SD1||18.6||7.53||17.1||[14.8–20.6]||0.||086|
|Minutes very active [min/day]||SD1||14.5||8.27||12.0||[9.24–17.7]||0.||002|
3.3Poincaré plot and physical activity variability (PHAV) parameters
To measure association between physical activity variability (PHAV) parameters accessed from the Poincaré plot, Pearson’s and Spearman’s correlation coefficients were calculated (Fig. 5).
Statistically significant negative correlations between changes (delta follow-up baseline) of lipid values and variability of Steps count, distance and minutes of very active physical activity were found, showing a decrease of lipid levels with higher long-term variability and area of an ellipse, but not short-term variability. Additionally, decreased body composition measurements were associated with minutes being fairly or very active. However, the SD1/SD2 ratio showed an opposite relationship with body composition measures indicating that more balanced long-term and short-term variability and possibly higher short-term variability could be associated with positive changes. Also, a similar association between HbA1c mmol/l and Steps and Distance SD1/SD2 ratio was observed.
3.4Vision of remote patient physical monitoring system
In previous paragraphs, we have described in detail how data were collected and analysed. However, many data operations (extraction, integration and analysis) were done manually. We have a vision to improve the proposed remote patient physical monitoring system in order the system could be used for all patients and we could remotely monitor patient health on a day-to-day basis. In this section, we will describe a vision of an integrated remote patient physical monitoring system compatible with existing SANTA-HIS and Lithuanian national eHealth framework.
VUHSK has developed hospital information system (SANTA-HIS), which integrates electronic health record, laboratory, picture archiving and communication system (PACS), registry of informed consent, staff and resource management, document management, quality management system, adverse events management system and many other systems that are necessary for effective health care services. SANTA-HIS system is integrated with Lithuanian national health information system which is called Electronic Health Services and Cooperation Infrastructure Information System (ESPBI IS). ESPBI IS was launched for patients and health care practitioners from 2015 and already has huge amounts of clinical data.
Figure 6 shows the vision of a remote patient physical monitoring system. The elements of the systems, which are shown in green, are already developed and used on a daily basis by health care practitioners.
Remote patient physical monitoring system and external medical devices integration system (shown in red) should be implemented and integrated into existing SANTA-HIS in order to automate data operations.
We strongly believe that we can benefit from the systems integrations because we will have all needed data in one place, and we can make detailed analyses and have a clear picture of the patient’s health. Systems integration will let health care practitioners to review and use patient data collected in other clinics and laboratories. Objective patient physical activity data and other parameters are essential in clinical work because it helps us know how the patient feels under real-life conditions and adheres doctor’s recommendations. We can find out why a patient can’t reach those goals and help solve the problem faster and more effectively.
Nonlinear physical activity variability parameters showed potential to include variability measurement in monitoring patients with prediabetic status. Using these parameters, we can visualize physical activity variability (Poincaré plot). Implementation architecture suggested in our paper can be used for implementation of these parameters into clinical practice and improve telemonitoring availability for physicians. In future research, the target value for variability in daily steps and other physical activity parameters should be evaluated not only in the larger group of prediabetic patients but also in healthy subjects to understand the impact on clinical parameters and health improvement.
Conflict of interest
None to report.
Association AD. 2. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes – 2020. Diabetes Care [Internet]. 2020 Jan 1 [cited 2021 Aug 7]; 43(Supplement 1): S14–31. Available from: https://care.diabetesjournals.org/content/43/Supplement_1/S14.
Association AD. Diagnosis and Classification of Diabetes Mellitus. Diabetes Care [Internet]. 2014 Jan 1 [cited 2021 Aug 7]; 37(Supplement 1): S81–90. Available from: https://care.diabetesjournals.org/content/37/supplement_1/s81.
Tuso P. Prediabetes and Lifestyle Modification: Time to Prevent a Preventable Disease. The Permanente Journal [Internet]. 2014 Jun 1 [cited 2021 Aug 4]; 18(3): 88. Available from: /pmc/articles/PMC4116271/.
Professor Wendy Brown Professor AJ, Activity P, Adrian Bauman Sesquicentenary Professor of Public Health Director PE, Fiona P, Burton NW. Development of Evidence-based Physical Activity Recommendations for Adults (18–64 years). 2013.
Physical activity [Internet] [cited 2021 Aug 7]. Available from: https://www.who.int/news-room/fact-sheets/detail/physical-activity.
HHS. Physical Activity Guidelines for Americans 2 nd edition.
Williams PT, Thompson PD. Walking Versus Running for Hypertension, Cholesterol, and Diabetes Mellitus Risk Reduction. Arteriosclerosis, Thrombosis, and Vascular Biology [Internet]. 2013 May [cited 2021 Aug 7]; 33(5): 1085–91. Available from: doi: 10.1161/ATVBAHA.112.300878.
Murtagh EM, Murphy MH, Boone-Heinonen J. Walking – the first steps in cardiovascular disease prevention. Current opinion in cardiology [Internet]. 2010 Sep [cited 2021 Aug 7]; 25(5): 490. Available from: /pmc/articles/PMC3098122/.
IM L, DM B. The importance of walking to public health. Medicine and science in sports and exercise [Internet]. 2008 [cited 2021 Aug 7]; 40(7 Suppl). Available from: https://pubmed.ncbi.nlm.nih.gov/18562968/.
Accumulating brisk walking for fitness, cardiovascular risk, and psychological health [Internet] [cited 2021 Aug 7]. Available from: https://wlv.openrepository.com/handle/2436/11110.
Kraus WE, Janz KF, Powell KE, Campbell WW, Jakicic JM, Troiano RP, et al. Daily Step Counts for Measuring Physical Activity Exposure and Its Relation to Health. Medicine and Science in Sports and Exercise [Internet]. 2019 Jun 1 [cited 2021 Aug 7]; 51(6): 1206–12. Available from: https://journals.lww.com/acsm-msse/Fulltext/2019/06000/Daily_Step_Counts_for_Measuring_Physical_Activity.15.aspx.
HHS. 2018 Physical Activity Guidelines Advisory Committee Scientific Report. 2018.
Daily Steps and Health Walking Your Way to Better Health [Internet] [cited 2021 Aug 7]. Available from: https://www.acsm.org/blog-detail/acsm-certified-blog/2019/06/14/walking-10000-steps-a-day-physical-activity-guidelines.
Althoff T, Sosič R, Hicks JL, King AC, Delp SL, Leskovec J. Large-scale physical activity data reveal worldwide activity inequality. Nature 2017 547:7663 [Internet]. 2017 Jul 10 [cited 2021 Aug 7]; 547(7663): 33–9. Available from: https://www.nature.com/articles/nature23018.
Degroote L, de Bourdeaudhuij I, Verloigne M, Poppe L, Crombez G. The accuracy of smart devices for measuring physical activity in daily life: Validation study. JMIR mHealth and uHealth. 2018 Dec 1; 6(12).
Aoyagi Y, Shephard RJ. Steps Per Day. Sports Medicine 2009 39:6 [Internet]. 2012 Nov 25 [cited 2021 Aug 7]; 39(6): 423–38. Available from: https://link.springer.com/article/10.2165/00007256-200939060-00001.
WE K, KF J, KE P, WW C, JM J, RP T, et al. Daily Step Counts for Measuring Physical Activity Exposure and Its Relation to Health. Medicine and science in sports and exercise [Internet]. 2019 Jun 1 [cited 2021 Aug 7]; 51(6): 1206–12. Available from: https://pubmed.ncbi.nlm.nih.gov/31095077/.
Rowlands AV, Gomersall SR, Tudor-Locke C, Bassett DR, Kang M, Fraysse F, et al. Introducing novel approaches for examining the variability of individuals’ physical activity. Journal of Sports Sciences. 2015; 33(5): 457–66.
Gill DP, Blunt W, Boa Sorte Silva NC, Stiller-Moldovan C, Zou GY, Petrella RJ. The HealtheSteps lifestyle prescription program to improve physical activity and modifiable risk factors for chronic disease: a pragmatic randomized controlled trial. BMC Public Health 2019 19:1 [Internet]. 2019 Jun 28 [cited 2021 Aug 8]; 19(1): 1–15. Available from: https://link.springer.com/articles/10.1186/s12889-019-7141-2.
Smith L, Hamer M, Ucci M, Marmot A, Gardner B, Sawyer A, et al. Weekday and weekend patterns of objectively measured sitting, standing, and stepping in a sample of office-based workers: the active buildings study. BMC Public Health 2015 15:1 [Internet]. 2015 Jan 17 [cited 2021 Aug 8]; 15(1): 1–9. Available from: https://link.springer.com/articles/10.1186/s12889-014-1338-1.
How do I set up my Fitbit device? [Internet]. Help.fitbit.com. 2021 [cited 6 August 2021]. Available from: https://help.fitbit.com/articles/en_US/Help_article/1873.htm#How.
Golińska AK. Poincaré Plots in Analysis of Selected Biomedical Signals. Studies in Logic, Grammar and Rhetoric. 201; 35(1): 11–27.
Brennan M, Palaniswami M, Kamen P. Do existing measures of Poincaré plot geometry reflect nonlinear features of heart rate variability? IEEE Transactions on Biomedical Engineering. 2001; 48(11): 1342–7.
Schneider M, Kraemmer MM, Weber B, Schwerdtfeger AR. Life events are associated with elevated heart rate and reduced heart complexity to acute psychological stress. Biological Psychology. 2021; 163: 108116.
Ahn JM. Heart Rate Variablity Assessment using the Poincaré Plot Parameters during Sleep. International Journal of Engineering Research and Technology. 2020; 13(8): 1921.
Iconaru EI, Ciucurel C. Hand grip strength variability during serial testing as an entropic biomarker of aging: A Poincaré plot analysis. BMC Geriatrics. 2020; 20(1).
Packiasabapathy S, Prasad V, Rangasamy V, Popok D, Xu X, Novack V, et al. Cardiac surgical outcome prediction by blood pressure variability indices Poincaré plot and coefficient of variation: An observational study [Internet]. Research Square. 2020. Available from: doi: 10.21203/rs.2.15887/v4.
Rechenberg K, Griggs S, Jeon S, Redeker N, Yaggi HK, Grey M. Sleep and glycemia in youth with type 1 diabetes [Internet]. Journal of Pediatric Health Care. Mosby; 2020 [cited 2021Aug9]. Available from: https://www.sciencedirect.com/science/article/pii/S0891524519305528.
Crenier L. Poincaré plot quantification for assessing glucose variability from continuous glucose monitoring systems and a new risk marker for hypoglycemia: Application to type 1 diabetes patients switching to continuous subcutaneous insulin infusion. Diabetes Technology & Therapeutics. 2014; 16(4): 247–54.
Kohnert K-D, Heinke P, Vogt L, Augstein P, Salzsieder E. Applications of variability analysis techniques for continuous glucose monitoring derived time series in diabetic patients. Frontiers in Physiology. 2018; 9.
Iconaru EI, Ciucurel MM, Georgescu L, Tudor M, Ciucurel C. The applicability of the Poincaré plot in the analysis of variability of reaction time during serial testing. International Journal of Environmental Research and Public Health. 2021; 18(7): 3706.
The jamovi project (2021). jamovi (Version 1.6) [Computer Software]. Retrieved from https://www.jamovi.org.
Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/.
Leandro Rodriguez-Linares, Xose Vila, Maria Jose Lado, Arturo Mendez, Abraham Otero and Constantino Antonio Garcia (2020). RHRV: Heart Rate Variability Analysis of ECG Data. R package version 4.2.6. https://CRAN.R-project.org/package=RHRV.
Taiyun Wei and Viliam Simko (2021). R package ‘corrplot’: Visualization of a Correlation Matrix (Version 0.90). Retrieved from https://github.com/taiyun/corrplot.