## ── 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