Plot Wevelet Variation
Power spectrum of HRV data with
Plotting Wevelet analysis
Figure 1 OIGAWA MARATHON 2023 Wevelet analysis, 2023-10-29
Power spectrum of HRV data with
Kizen Sasaki
2023 November 30
2023 December 08
Figure 1 OIGAWA MARATHON 2023 Wevelet analysis, 2023-10-29
---
title: "Plot Wevelet Variation"
subtitle: "Power spectrum of HRV data with "
author:
- Kizen Sasaki
date: 2023-11-30
date-modified: today
date-format: 'YYYY MMMM DD'
language:
title-block-published: First version
titel-block-modified: This version
format:
html:
html-math-method: katex
code-tools: true
code-fold: true
code-link: true
pdf:
geometry:
- top=30mm
- left=30mm
docx: default
code-line-numbers: true
---
```{r}
#| echo: false
#| message: false
library(RHRV)
csv_xls <- "~/Documents/brevet/csv_xls/time_hr_1.csv"
data <- read.csv(csv_xls, header = TRUE)
hr <- data[,2]
rr <- 1/hr
rr = rr*60*1000
```
```{r}
#| echo: false
#| message: false
hrv.data <- CreateHRVData()
hrv.data <- SetVerbose(hrv.data, TRUE)
hrv.data$Beat$RR = rr
hrv.data$Beat$niHR = hr
hrv.data$Beat$Time = cumsum(rr)/1000
hrv.data <- BuildNIHR(hrv.data)
hrv.data <- FilterNIHR(hrv.data)
hrv.data <- InterpolateNIHR(hrv.data, freqhr = 4)
hrv.data <- CreateFreqAnalysis(hrv.data)
```
```{r}
#| message: false
#| echo: false
hrv.data <- CalculatePowerBand(hrv.data, indexFreqAnalysis = 1,
type = "wavelet", wavelet = "la8", bandtolerance = 0.01,
relative = FALSE)
```
### Plotting Wevelet analysis
@fig-wevelet OIGAWA MARATHON 2023 Wevelet analysis, 2023-10-29
```{r}
#| echo: false
#| message: false
#| fig-width: 8
#| fig-height: 11.31
#| label: fig-wevelet
#| fig.cap: "Plot obtained with the PlotPowerBand for the wavelet analysis."
#| fig-subcap:
#| - ymax = 700, ymaxratio = 50
#| - ymax = 1000, ymaxratio = 50
#| layout-ncol: 2
PlotPowerBand(hrv.data, indexFreqAnalysis = 1,
ymax = 700, ymaxratio = 50)
PlotPowerBand(hrv.data, indexFreqAnalysis = 1,
ymax = 1000, ymaxratio = 50)
```