-
Notifications
You must be signed in to change notification settings - Fork 3
/
shanmdphd.Rmd
56 lines (41 loc) Β· 1.24 KB
/
shanmdphd.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
---
title: "BMIμ λ°λ₯Έ Primary coronary intervention μ entry κ²½λ‘μ κ²°μ "
author: "νμ±ν"
header-includes:
- \usepackage{kotex}
mainfont: NanumGothic
output:
pdf_document:
latex_engine: xelatex
---
## μλ‘
μ΄ λΆμμμλ 857λͺ
μ acute coronary syndrome νμμμ primary coronary intervention μ entry κ²½λ‘μ κ²°μ μ΄ BMIμ λ°λΌ μν₯μ λ°λμ§ μμλ³΄κ³ μ νλ€.
## λ°©λ²
`moonBook::acs` μλ£λ₯Ό μ¬μ©νμλ€.
## κ²°κ³Ό
ν° μ°¨μ΄κ° μμλ€.
```{r}
library(tidyverse)
library(readxl)
raw_acs <- read_excel('acs.xlsx') %>%
mutate(BMI = as.numeric(BMI))
raw_acs %>%
group_by(entry) %>%
summarise(median_BMI = median(BMI, na.rm = TRUE),
sd_BMI = sd(BMI, na.rm = TRUE))
```
```{r}
qplot(x = entry, y = BMI, data = raw_acs, geom = 'boxplot')
```
μ±λ³κ³Ό μ§λ¨μ κ³ λ €ν΄λ λμ λλ μ°¨μ΄λ μμλ€.
```{r}
raw_acs %>%
group_by(Dx, sex, entry) %>%
summarise(median_BMI = median(BMI, na.rm = TRUE),
sd_BMI = sd(BMI, na.rm = TRUE))
```
```{r}
qplot(x = entry, y = BMI, data = raw_acs, geom = 'boxplot', facets = sex ~ Dx)
```
## κ²°λ‘
BMIλ₯Ό λΉλ‘―νμ¬ μ§λ¨, μ±λ³ λ±μ μν₯μ PCIμ μ§μ
κ²½λ‘μ μν₯μ λ―ΈμΉμ§ μμλ€.