## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.2 ──
## ✔ ggplot2 3.4.0      ✔ purrr   1.0.0 
## ✔ tibble  3.1.8      ✔ dplyr   1.0.10
## ✔ tidyr   1.2.1      ✔ stringr 1.5.0 
## ✔ readr   2.1.3      ✔ forcats 0.5.2 
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## Loading required package: timechange
## 
## 
## Attaching package: 'lubridate'
## 
## 
## The following objects are masked from 'package:base':
## 
##     date, intersect, setdiff, union
## Warning: `qplot()` was deprecated in ggplot2 3.4.0.
## Warning: Removed 18 rows containing missing values (`geom_point()`).

## Warning: Removed 18 rows containing missing values (`geom_point()`).

## Warning: Removed 1 rows containing missing values (`geom_point()`).

## [1] 329  14

## [1] 118.8489
## [1] 73.31833
## [1] 119
## [1] 74
## [1] 74.29878
## [1] 72
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 18 rows containing non-finite values (`stat_bin()`).
## Warning: Removed 18 rows containing non-finite values (`stat_bin()`).

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 18 rows containing non-finite values (`stat_bin()`).
## Removed 18 rows containing non-finite values (`stat_bin()`).

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 1 rows containing non-finite values (`stat_bin()`).
## Warning: Removed 1 rows containing non-finite values (`stat_bin()`).

## 
##  Shapiro-Wilk normality test
## 
## data:  gdf$press_high
## W = 0.98783, p-value = 0.01028

相関

qplot(gdf$press_high,gdf$press_low)
## Warning: Removed 18 rows containing missing values (`geom_point()`).

cor(gdf$press_high,gdf$press_low,use = "complete.obs")
## [1] 0.7518274
qplot(gdf$press_high,gdf$BR_bpm)
## Warning: Removed 18 rows containing missing values (`geom_point()`).

cor(gdf$press_high,gdf$BR_bpm,use = "complete.obs")
## [1] -0.06123426
qplot(gdf$press_low,gdf$BR_bpm)
## Warning: Removed 18 rows containing missing values (`geom_point()`).

cor(gdf$press_low,gdf$BR_bpm,use = "complete.obs")
## [1] 0.04705072

血圧=u_date(d)+hour[t]+bpm+state1 bpm=u_date(d)+hour[t]+hstate+pressure+state1