On the improvement of heart rate prediction using the combination of singular spectrum analysis and copula-based analysis approach

Asieh Namazi

doi:10.7717/peerj.14601

On the improvement of heart rate prediction using the combination of singular spectrum analysis and copula-based analysis approach

Asieh Namazi

Department of Physical Education and Sport Science, Iran University of Science and Technology Tehran, Tehran, Iran

DOI: 10.7717/peerj.14601

Published: 2022-12-19
Accepted: 2022-11-29
Received: 2022-10-19

Academic Editor: Cosme Buzzachera

Subject Areas: Cardiology, Public Health, Sports Medicine
Keywords: Heart rate, Machine learning, Wearable sensors, Artificial intelligence, Prediction

Copyright: © 2022 Namazi
Licence: This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.

Cite this article: Namazi A. 2022. On the improvement of heart rate prediction using the combination of singular spectrum analysis and copula-based analysis approach. PeerJ 10:e14601 https://doi.org/10.7717/peerj.14601

Abstract

In recent years, many people have been working from home due to the exceptional circumstances concerning the coronavirus disease 2019 (COVID-19) pandemic. It has also negatively influenced general health and quality of life. Therefore, physical activity has been gaining much attention in preventing the spread of Severe Acute Respiratory Syndrome Coronavirus. For planning an effective physical activity for different clients, physical activity intensity and load degree needs to be appropriately adjusted depending on the individual’s physical/health conditions. Heart rate (HR) is one of the most critical health indicators for monitoring exercise intensity and load degree because it is closely related to the heart rate. Heart rate prediction estimates the heart rate at the next moment based on now and other influencing factors. Therefore, an accurate short-term HR prediction technique can deliver efficient early warning for human health and decrease the happening of harmful events. The work described in this article aims to introduce a novel hybrid approach to model and predict the heart rate dynamics for different exercises. The results indicate that the combination of singular spectrum analysis (SSA) and the Clayton Copula model can accurately predict HR for the short term.

Introduction

In recent years, many people have been working from home due to the exceptional circumstances concerning the coronavirus disease 2019 (COVID-19) pandemic. It has negatively influenced general health and quality of life. Physical inactivity imposes a significant economic burden on the health care system, and the COVID-19 pandemic has harmed physical activity (Puccinelli et al., 2021). Those who reduced their physical activity level had the highest mood disorders. Nowadays, physical activity has been gaining much attention in preventing the spread of Severe Acute Respiratory Syndrome Coronavirus. Over the last several decades, data obtained from the linear relationship between heart rate and oxygen consumption are taken to estimate maximal oxygen uptake ( $V_{O} 2 m a x$ ) (Jlstrand & Ryhming, 1954). For planning an effective physical activity for different clients, physical activity intensity and load degree needs to be appropriately adjusted depending on the individual physical/health conditions (Böning, Gönen & Maassen, 1984; Hoogendoorn et al., 1999).

Heart rate (HR) is one of the most critical health indicators for monitoring the physical load degree because it is closely related to the heart rate (Dong, 2016). Therefore, studying heart rate and its influencing factors is of great importance. The human body does not adapt to activities immediately, and the response heart rate is delayed and increases after a certain time of physical activity (Coats et al., 1992). Heart rate prediction estimates the heart rate at the next moment based on now and other influencing factors. Therefore, an accurate short-term HR prediction technique can deliver efficient early warning for human health and decrease the happening of harmful events or premature death (Malik & Camm, 1990; Janssen, 2012,). However, traditional heart rate prediction methods cannot meet the requirements of high-precision dynamic prediction. Thus, reliable and effective heart rate prediction is significant for preventing and controlling some cardiovascular diseases (Lei, Cai & Hua, 2021). In the last decade, various heart rate prediction models have been examined. Some models are based on a Hammerstein model’s differential equations or variants. Cheng et al. (2008) present a nonlinear condition space model for forecasting a person’s heart rate behavior based on running speed on a treadmill. Paradiso et al. (2012) employ the same model to control the bicycle heart rate. Mohammad et al. (2011) introduced the fuzzy Takagi-Sugeno model dealing with 12 parameters and is commonly used to optimization of physical activity for older, untrained people. Other model-based systems exist for running and cycling on different training devices (Su et al., 2007; Koenig et al., 2009; Leitner et al., 2014).

The work described in this article aims to introduce a novel hybrid approach to model and predict the heart rate dynamics for different exercises. Modiri et al. (2018) introduced a hybrid approach, a combination of deterministic and stochastic methods (the combination of singular spectrum analysis (SSA) and Copula-based analysis) to predict space geodetic data for the first time. The main advantage of using the SSA+Copula approach is fully describing all relations underlying different variables. SSA is not limited to using predetermined functions such as a sine wave as a base. Instead, SSA uses data-driven basis functions to extract basic components of the data series and applies a classification method to examine the relationship between the derived elements (Broomhead & King, 1986; Vautard, Yiou & Ghil, 1992). The Copula-based analysis approach exploits linear and nonlinear dependencies between variables. It is a powerful and efficient tool for dealing with multi-dimensional data and modeling the dependency structure between parameters (Joe, 1997). This study applied the combination of SSA and Copula for the first time as a novel deterministic-stochastic tool for HR prediction using only HR past data. In this method, the deterministic part is estimated by SSA, whereas Copula is used for modeling the stochastic part. This study can support the development of automated algorithms for real-time health monitoring of athletes and preventive identification of individuals at increased risk for sudden sports-related events, as suggested by Sbrollini et al. (2019).

The remainder of this article is organized as follows. “Methodology” discusses how to describe the methodology and prediction algorithms. Next, “Data Set And Data Analysis” illustrates the data set and analysis where the input data is fully described. Then, “Results” discusses the results, and finally, the conclusion is written in “Conclusions”.

Methodology

In this study, the integration of Copula-based analysis and SSA is developed for precisely predicting heart rate.

Singular spectrum analysis

SSA is a non-parametric technique used for many tasks, such as trend detection and extraction, de-noising, forecasting, and change-point detection. SSA works with arbitrary statistical processes, whether linear or nonlinear, stationary or non-stationary, Gaussian or non-Gaussian (Ghil et al., 2002; Hassani & Thomakos, 2010; Modiri et al., 2018).

The SSA trend extraction method can be described as having four stages (Hoseini et al., 2020; Hoseini, 2022):

Stage 1: Embedding. A sliding window with length L forms a trajectory matrix X over the time series $f_{i}$ described as follows:

(1) $\begin{matrix} X & = (x_{i j})_{i, j = 1}^{L, K} = & [\begin{matrix} f_{1} & f_{2} & f_{3} & \dots & f_{K} \\ f_{2} & f_{3} & f_{4} & \dots & f_{K + 1} \\ f_{3} & f_{4} & f_{5} & \dots & f_{K + 2} \\ ⋮ & ⋮ & ⋮ & ⋱ & ⋮ \\ f_{L} & f_{L + 1} & f_{L + 2} & \dots & f_{N} \end{matrix}], {\begin{matrix} 1 < L < K \\ K = N - L + 1 \end{matrix} \end{matrix}$ where, X is a symmetric matrix with identical elements on the anti-diagonals (Korobeynikov, 2009). Also, the K represents the indices of the L-dimensional time-lagged vectors.

Stage 2: Decomposition. The singular value decomposition (SVD) is applied and decomposed into the trajectory matrix X (Wall, Rechtsteiner & Rocha, 2003; Modiri, 2022), as $X = U Σ V^{T}$ . The matrices U and V are singular and orthonormal vectors. Also, $Σ$ denotes a diagonal matrix of the singular values. Then, the SVD of the X matrix can be obtained from eigenvalues and eigenvectors of the lag-covariance matrix defined as $S = X X^{T}$ of size $L \times L$ . S is named empirical orthogonal functions (EOFs) and contains eigenvalues denoted by $λ_{1} \geq λ_{2} \geq \dots \geq λ_{L} \geq 0$ . Then, the corresponding orthonormal eigenvectors of S can be indicated by $U_{1}, U_{2}, \dots U_{L}$ and $V_{1}, V_{2}, \dots V_{L}$ for left and right, respectively. Thus, the SVD obtained a representation $X = X_{1} + X_{2} + \dots + X_{d}$ , where the trajectory matrix can be written as: $X_{i} = \sqrt{λ_{i}} U_{i} V_{i}^{T}$ ( $i = 1, \dots, d)$ . The $\sqrt{λ_{i}}$ is the corresponding singular value of $X_{i}$ and d is the number of non-zero eigenvalues of $λ_{i}$ . Each of the matrices $X_{i}$ describes some components of the time series, and summing them should deliver $X$ .

Stage 3: Grouping. The indices are split into two groups, trend and residual, as follows:

(2) ${\begin{array}{l} X_{t r e n d} = X_{1} + X_{2} + \dots + X_{I} = ({\hat{x}}_{i, j})_{i, j = 1}^{L, K} \\ X_{r e s i d u a l} = X_{I + 1} + X_{I + 2} + \dots + X_{d} \end{array}$ where $I \leq d$ .

Stage 4: Reconstruction. To reconstruct the trend, averaging the anti-diagonal entries of the matrix $X_{t r e n d}$ in Eq. (2) should be done. Let L < K, then the trend of time series $G = (g_{1}, g_{2}, \dots, g_{N})$ is:

(3) $g_{i} = {\begin{matrix} \frac{1}{i} Σ_{m = l}^{i} {\hat{X}}_{m, i - m + 1} 1 \leq i \leq L \\ \frac{1}{L} Σ_{m = l}^{L} {\hat{X}}_{m, i - m + 1} L \leq i \leq K \\ \frac{1}{N - i + 1} Σ_{m = i - K + 1}^{N - K + 1} {\hat{X}}_{m, i - m + 1} K \leq i \leq N \end{matrix}$ where ${\hat{x}}_{i, j}$ is an estimation of the element $f_{i + j - 1}$ of the original time series.

Furthermore, the reconstructed trajectory matrix might not keep the Hankel condition, and thus ${\hat{X}}_{i, j} \neq {\hat{X}}_{j, i}$ (Li, Liu & Razavilar, 1997).

Copula-based analysis

The Copula-based analysis exploit the dependency structure between random variables independent of their marginal distributions. It is a flexible tool that greatly improves the capture of the whole correlation pattern (Embrechts, McNeil & Straumann, 2002). The Copula-based analysis model is used for various types of studies, e.g., economics, hydrology, meteorology, and geodesy (Patton, 2006, 2009; Bárdossy & Li, 2008; Bárdossy & Pegram, 2009; Laux et al., 2011; Modiri, 2015; Modiri et al., 2018, 2020).

A Copula model is revealed as a function C from $[0, 1]^{2}$ to [0, 1]. The general properties of Copula can be itemized as follows (Genest & Rivest, 1993; Nelsen, 2007):

$C (u, 0) = C (0, v)$ ,
$C (u, 1) = u$ and $C (1, v) = v$ ,
The Copula is unique when the marginals are continuous functions.

Besides, it should be mentioned that the variables transformed by any monotonic increasing functions will not affect its Copula. The empirical Copula is an estimator for the unknown theoretical Copula and is defined on the rank space as (Genest & Favre, 2007):

(4) $C_{e} (u, v) = \frac{1}{n} \sum_{i = 1}^{n} 1 (\frac{r_{i}}{n + 1} \leq u, \frac{s_{i}}{n + 1} \leq v)$ where, $(r_{1}), \dots, (r_{n})$ mean the pairs of ranks of the variable $(x_{1}), (x_{2}), \dots, (x_{n})$ . $(s_{1}), \dots, (s_{n})$ denote the pairs of ranks of the variable $(y_{1}), (y_{2}), \dots, (y_{n})$ . $n$ is the length of the data vector. A total of 1 is the indicator function. The most commonly used Copula model is the Archimedean Copula model, which demonstrates a-/symmetrical upper/lower tail dependency and can be estimated directly in a simple form.

(5) $C (u, v) = ϕ^{- 1} {ϕ (u) + ϕ (v), θ}$ where $θ$ is the Copula parameter and the function $ϕ$ is the Copula generator with the following characteristics (Nelsen, 2007).

Clayton Copula. The Clayton Copula is an asymmetric Archimedean Copula, which can be used for the variables with a lower tail dependence structure. In order to define the Clayton Copula model, first, its generator should be estimated as follow (Clayton, 1978):

(6) $ϕ^{C l} (x) = \frac{1}{θ} (t^{- θ} - 1)$ where, $θ$ is the Clayton Copula parameter given over $- 1 \leq θ$ . Therefore, the Clayton Copula can be simplified to Nelsen (2007):

(7) $C_{θ} (u, v) = max [(u^{- θ} + v^{- θ} - 1), 0]^{- \frac{1}{θ}}$

Frank Copula. The Frank Copula is mainly used to model the symmetrical dependency between the variables. So, the generator of the Frank Copula model is indicated by Nelsen (2007):

(8) $ϕ^{F r} (t) = - \ln {\frac{e^{- θ t} - 1}{e^{- θ} - 1}}$

The Frank Copula parameter is defined over $- \infty < θ < \infty$ and its formula is given by Nelsen (2007):

(9) $C_{θ} (u, v) = \frac{1}{θ} \ln (1 + \frac{(e^{- θ u} - 1) (e^{- θ v})}{e^{- θ} - 1})$

In the case of independence Copula, the Frank Copula parameter is equal to zero, and the higher parameter shows a linear dependency structure between variables.

Gumbel Copula. The Gumbel Copula is also an asymmetric Archimedean Copula, which is used to capture the upper tail dependence. The Gumbel Copula generator is shown as follows:

(10) $ϕ^{G u} (t) = (- \ln t)^{θ}$

The Gumbel Copula parameter is defined over $1 \leq θ$ and its formula is given by Nelsen (2007): (11) $C_{θ} (u, v) = e^{- {((- \ln {(u)}^{θ}) + (- \ln {(v)}^{θ}))}^{\frac{1}{θ}}}$

In the case of independence structure between the variables, the Gumbel Copula parameter is equal to 1 and if the $θ \to \infty$ is the comonotonicity Copula. All described Copula models are obtained by estimating the Copula parameters and, in this study, are derived from the maximum likelihood (ML) estimation described by Joe (1997).

Prediction algorithm

The Copula + SSA is developed as a novel hybrid method to predict different parameters. It uses training data to model the data’s underlying structure.

As shown in Fig. 1, the prediction algorithm consists of two parts; SSA models the deterministic part, and the anomaly (the residual of subtracted data by the reconstructed time series) part is modeled by Coupla-based analysis. Then, the SSA extrapolates the deterministic component, and the residual part is predicted by Copula-based analysis. In the end, the SSA extrapolated trend and predicted anomaly would be combined (Modiri, 2021).

Figure 1: Copula + SSA algorithm for modeling and predicting the data by using the information within the past data.

Download full-size image

DOI: 10.7717/peerj.14601/fig-1

Data set and data analysis

In this article, the HR data are from the online available data set: ‘Sport Database: Cardiorespiratory data acquired through wearable sensors while practicing sports’ (Sbrollini et al., 2019). All subjects gave their informed consent prior to data collection and acquisitions, which were undertaken in compliance with the ethical principles of Helsinki Declaration and approved by the institutional expert committee. The characteristics of the HR signal is shown in Table 1.

Table 1:

Characteristics of the HR signals.

Signal	Sampling frequency	Amplitude range
Heart rate	1 Hz	25–240 bpm

DOI: 10.7717/peerj.14601/table-1

In this study, the data from ten subjects who did running test is used. A particular survey protocol was characterized for running. It contains a 6.1 Km route around the city of Ancona, Italy (see Fig. 2 of the reference Sbrollini et al. (2019)). The protocol contains four phases with different slopes:

Initial flat phase (1.3 km long with a 0)
Uphill phase (1.2 km long with a $6.8 %$ slope).
Downhill phase (1 km long with a $7.2 %$ slope).
Final flat phase (2.6 km long with a $0 %$ slope)

All subjects were supposed to be healthy (i.e., no previous history of cardiorespiratory diseases and not taking any drug) at the acquisition time. Also, the demographic data of the Sport Database is given by Table 2. All subjects gave their informed consent prior to data collection and acquisitions, which were undertaken in compliance with the ethical principles of the Helsinki Declaration and approved by the institutional expert committee.

Table 2:

Demographic data of sport database.

Number of subjects	Number of CRD	Gender M/F	Age (years)	Weight (kg)	Height (cm)	Smoking NO/YES	Alcohol consumption NO/SOMETIMES	Weekly training rate
10	10	9/1	22 ± 3	70 ± 6	179 ± 7	5/5	1/9	3 ± 1

DOI: 10.7717/peerj.14601/table-2

As an example, Fig. 2 shows the HR of subject 1, where this subject has an initial flat phase: from 00:00:00 to 00:07:45, then has an uphill phase: from 00:07:45 to 00:18:32, after that, the downhill phase: from 00:18:32 to 00:25:47, in the end, the subject starts the final flat phase: from 00:25:47 to 00:43:54.

Figure 2: The original HR data (black) and the reconstructed HR data (blue).

Download full-size image

DOI: 10.7717/peerj.14601/fig-2

The measured HR data could be divided into two parts. The first part deals with the trend and deterministic term of HR. So, SSA is employed to model the HR signal’s main behavior (trend + deterministic term). Then, the difference between the HR and SSA estimated data is modeled using the Copula-based analysis approach. After that, the main behavior of HR is extrapolated using the SSA a priori model. Correspondingly, the anomaly part is predicted using the Copula-based model. Finally, the anomaly solution is added to the SSA-forecasted data.

Therefore, the analysis of the data is divided into two main steps.

SSA trend estimation and extrapolation

As it shown in Fig. 2, the first 1,500 epochs are taken as training data with sliding step 1. Then, the first 60 singular value is used to reconstruct the HR. Then, the reconstructed HP is extrapolated up to the next 30 s.

Copula anomaly modeling and prediction

The subtracted original HR and SSA reconstructed HR is formed using a window length of 30 and epoch delay of 1 epoch (see Fig. 3):

Figure 3: The difference between original and reconstructed time series for HR.

Download full-size image

DOI: 10.7717/peerj.14601/fig-3

The marginal distribution of each column of the matrix is computed.
It is transformed into the rank space.
The empirical Copula between column i and i + 1 is computed.
Eventually, the Clayton, Frank, and Gumbel Copula are fitted to empirical Copula by applying the appropriate goodness of fit test. After that, 3,000 random data are sampled from the conditional Copula CDF and transformed back to the data space using inverse marginal. Thus, the mean value of 3,000 random data is computed.

Finally, the HR predicted data is the sum of the results of predicted periodic terms using SSA and predicted anomaly using the Copula-based model (see Fig. 4).

Figure 4: A total of 30 s HR prediction using 1,500 past epochs.

Download full-size image

DOI: 10.7717/peerj.14601/fig-4

Results

We utilized 1,500 epochs of HR data for the 30-epochs-ahead prediction. To verify the reliability of this method, the results were compared with the original data. As a result, almost 502 sets of predicted HR have been done. So, Fig. 5 shows the mean value of MAE for different prediction algorithms (SSA+Clayton (green), SSA+Frank (blue), SSA+Clayton (red)).

Figure 5: Mean value of MAE of HR prediction for 508 data set with the unit is [bpm].

Download full-size image

DOI: 10.7717/peerj.14601/fig-5

In Fig. 5, the results of SSA+Copula are displayed in green, blue, and red for Clayton, Frank, and Gumbel Copula, respectively. The result of SSA + Clayton Copula shows slightly better performance.

The results demonstrated in this study, SSA + Copula technique performed better compared to the techniques shown in the recent studies by Alharbi et al. (2021), Monkaresi, Calvo & Yan (2013), Oyeleye et al. (2022). The mean of MAE is lower than 2 bpm up to the first seven epochs and slightly increases up to 3.3 bpm for the subsequent 30 epochs. Among the proposed techniques, SSA + Clayton Copula method performs better than SSA + Frank and SSA + Gumbel Copula.

Figure 6 presents the absolute error of 30-epochs-ahead prediction for almost 805 predicted data sets. The proposed method did not always perform better, and the results were not as good as we expected. This may have been caused by changes in the amplitudes of the deterministic terms where the SSA could not capture all features. In order to predict more precisely, we would have to increase the interval of training time.

Figure 6: Absolute errors of the predicted HR using SSA+Clayton Copula compared with original HRdata. The unit is [bpm].

Download full-size image

DOI: 10.7717/peerj.14601/fig-6

Conclusions

This study proposed a new prediction technique to predict HR, one of the most important health indicators for planning physical activity. Then, novel prediction techniques have been used to detect early warnings for human health. The best-developed model for each case with the lowest MAE is used to evaluate the prediction methods in real-time. The results illustrate that the proposed method could efficiently and precisely predict HR. As demonstrated, the Copula-based analysis successfully increases the accuracy of HR prediction by modeling the stochastic part of the HR. However, as described in Modiri et al. (2018), the primary error contributions come from the SSA extrapolation part. So, further investigations about the SSA training time will be required to clarify this issue. Also, the Copula-based analysis can potentially add more physical parameters for HR prediction. Therefore, future investigations are needed to get benefit from other parameters.

[1] Alharbi A, Alosaimi W, Sahal R, Saleh H. 2021. Real-time system prediction for heart rate using deep learning and stream processing platforms. Complexity 2021:1-9

[2] Bárdossy A, Li J. 2008. Geostatistical interpolation using copulas. Water Resources Research 44(7):W07412

[3] Bárdossy A, Pegram G. 2009. Copula based multisite model for daily precipitation simulation. Hydrology and Earth System Sciences 13(12):2299-2314

[4] Broomhead DS, King GP. 1986. Extracting qualitative dynamics from experimental data. Physica D: Nonlinear Phenomena 20(2–3):217-236

[5] Böning D, Gönen Y, Maassen N. 1984. Relationship between work load, pedal frequency, and physical fitness. International Journal of Sports Medicine 5(2):92-97

[6] Cheng TM, Savkin AV, Celler BG, Su SW, Lu Wang. 2008. Nonlinear modeling and control of human heart rate response during exercise with various work load intensities. IEEE Transactions on Biomedical Engineering 55(11):2499-2508

[7] Clayton DG. 1978. A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence. Biometrika 65(1):141-151

[8] Coats A, Adamopoulos S, Radaelli A, McCance A, Meyer TE, Bernardi L, Solda PL, Davey P, Ormerod O, Forfar C. 1992. Controlled trial of physical training in chronic heart failure. Exercise performance, hemodynamics, ventilation, and autonomic function. Circulation 85(6):2119-2131

[9] Dong J-G. 2016. The role of heart rate variability in sports physiology. Experimental and Therapeutic Medicine 11(5):1531-1536

[10] Embrechts P, McNeil A, Straumann D. 2002. Correlation and dependence in risk management: properties and pitfalls. Risk Management: Value at Risk and Beyond 1:176-223

[11] Genest C, Favre A-C. 2007. Everything you always wanted to know about copula modeling but were afraid to ask. Journal of Hydrologic Engineering 12(4):347-368

[12] Genest C, Rivest L-P. 1993. Statistical inference procedures for bivariate archimedean copulas. Journal of the American Statistical Association 88(423):1034-1043

[13] Ghil M, Allen MR, Dettinger MD, Ide K, Kondrashov D, Mann ME, Robertson AW, Saunders A, Tian Y, Varadi F, Yiou P. 2002. Advanced spectral methods for climatic time series. Reviews of Geophysics 40(1):1-3

[14] Hassani H, Thomakos D. 2010. A review on singular spectrum analysis for economic and financial time series. Statistics and Its Interface 3(3):377-397

[15] Hoogendoorn WE, Van Poppel MN, Bongers PM, Koes BW, Bouter LM. 1999. Physical load during work and leisure time as risk factors for back pain. Scandinavian Journal of Work, Environment & Health 25(5):387-403

[16] Hoseini M. 2022. On the remote sensing of the atmosphere and ocean using direct and reflected GNSS signals.

[17] Hoseini M, Alshawaf F, Nahavandchi H, Dick G, Wickert J. 2020. Towards a zero-difference approach for homogenizing GNSS tropospheric products. GPS Solutions 24(1):1-12

[18] Janssen I. 2012. Health care costs of physical inactivity in canadian adults. Applied Physiology, Nutrition, and Metabolism 37(4):803-806

[19] Jlstrand P, Ryhming I. 1954. A nomogram for calculation of aerobic capacity from pulse rate during submaximal work. Journal of Applied Physiology 7(2):218-221

[20] Joe H. 1997. Multivariate models and multivariate dependence concepts. Boca Raton: CRC press.

[21] Koenig AC, Somaini L, Pulfer M, Holenstein T, Omlin X, Wieser M, Riener R. 2009. Model-based heart rate prediction during lokomat walking.

[22] Korobeynikov A. 2009. Computation-and space-efficient implementation of SSA. ArXiv preprint.

[23] Laux P, Vogl S, Qiu W, Knoche HR, Kunstmann H. 2011. Copula-based statistical refinement of precipitation in RCM simulations over complex terrain. Hydrology and Earth System Sciences 15(7):2401-2419

[24] Lei T, Cai Z, Hua L. 2021. Training prediction and athlete heart rate measurement based on multi-channel PPG signal and SVM algorithm. Journal of Intelligent & Fuzzy Systems 40(4):7497-7508

[25] Leitner T, Kirchsteiger H, Trogmann H, del Re L. 2014. Model based control of human heart rate on a bicycle ergometer.

[26] Li Y, Liu KR, Razavilar J. 1997. A parameter estimation scheme for damped sinusoidal signals based on low-rank hankel approximation. IEEE Transactions on Signal Processing 45(2):481-486

[27] Malik M, Camm AJ. 1990. Heart rate variability. Clinical Cardiology 13(8):570-576

[28] Modiri S. 2015. Copula-based analysis of correlation structures in case of GRACE coefficients. PhD thesis, Universitätsbibliothek der Universität Stuttgart

[29] Modiri S. 2021. On the improvement of Earth orientation parameters estimation: using modern space geodetic techniques. Germany: Technische Universitaet Berlin.

[30] Modiri E. 2022. Climate change and concurrency of extreme events. In: Climate Change. Berlin: Springer. 151-170

[31] Modiri S, Belda S, Heinkelmann R, Hoseini M, Ferrándiz JM, Schuh H. 2018. Polar motion prediction using the combination of SSA and Copula-based analysis. Earth, Planets and Space 70(1):1-18

[32] Modiri S, Belda S, Hoseini M, Heinkelmann R, Ferrándiz JM, Schuh H. 2020. A new hybrid method to improve the ultra-short-term prediction of LOD. Journal of Geodesy 94(2):1-14

[33] Mohammad S, Guerra TM, Grobois JM, Hecquet B. 2011. Heart rate control during cycling exercise using Takagi-Sugeno models. IFAC Proceedings Volumes 44(1):12783-12788

[34] Monkaresi H, Calvo RA, Yan H. 2013. A machine learning approach to improve contactless heart rate monitoring using a webcam. IEEE Journal of Biomedical and Health Informatics 18(4):1153-1160

[35] Nelsen RB. 2007. An introduction to Copulas. Berlin: Springer Science & Business Media.

[36] Oyeleye M, Chen T, Titarenko S, Antoniou G. 2022. A predictive analysis of heart rates using machine learning techniques. International Journal of Environmental Research and Public Health 19(4):2417

[37] Paradiso M, Pietrosanti S, Scalzi S, Tomei P, Verrelli CM. 2012. Experimental heart rate regulation in cycle-ergometer exercises. IEEE Transactions on Biomedical Engineering 60(1):135-139

[38] Patton AJ. 2006. Modelling asymmetric exchange rate dependence. International Economic Review 47(2):527-556

[39] Patton AJ. 2009. Copula-based models for financial time series. In: Handbook of Financial Time Series. Berlin: Springer. 767-785

[40] Puccinelli PJ, da Costa TS, Seffrin A, de Lira CAB, Vancini RL, Nikolaidis PT, Knechtle B, Rosemann T, Hill L, Andrade MS. 2021. Reduced level of physical activity during COVID-19 pandemic is associated with depression and anxiety levels: an internet-based survey. BMC Public Health 21(1):1-11