교수님, 안녕하세요.
다름이 아니라 sphagetti plot을 그리려고 하는데 잘되지가 않아서 문의드립니다.
환자의 경과에 따른 혈압변화를 나타내고자 sphagetti plot을 그리려고합니다.
가능할지요.....
Comment 3
-
cardiomoon
2019.07.23 08:29
-
cardio_H
2019.07.23 09:04
위의 파일은 이런 그래프를 그리고 싶어 올린 예시 그래프입니다.
-
shanmdphd
2019.08.14 11:21
이렇게 하면 어떨까요?
# setup ----
library(tidyverse)
library(lubridate)
df <- read.csv('AMI_registry_BP.csv', as.is=TRUE) %>%
as_tibble()
tidy_df <- function(df, column = 'Date', subtext) {
df %>%
select(id, starts_with(column)) %>%
gather(visit, column, -id) %>%
filter(grepl('[0-9]$', visit)) %>%
mutate(visit = sub(subtext, '', visit) %>% as.numeric)
}
# combine all data ----
df_raw <- left_join(tidy_df(df, 'Date', 'Date_') %>% rename(date = column),
tidy_df(df, 'SBP', 'SBP_') %>% rename(sbp = column),
by = c('id', 'visit')) %>%
left_join(tidy_df(df, 'DBP', 'DBP_') %>% rename(dbp = column),
by = c('id', 'visit')) %>%
mutate(date = lubridate::ymd(date)) %>%
filter(!is.na(date) & !is.na(sbp) & !is.na(dbp)) %>%
print()
# cleaning strange data ----
df_final <- df_raw %>%
filter(!id %in% (df_raw %>% filter(sbp > 300) %>% .$id)) %>%
filter(!id %in% (df_raw %>% filter(dbp > 200) %>% .$id))
# SBP plot ----
df_final %>%
ggplot(aes(date, sbp, group = id)) +
geom_line(alpha = 0.1) +
theme_bw()
ggsave('systolic-bp.png')
# DBP plot ----
df_final %>%
ggplot(aes(date, dbp, group = id)) +
geom_line(alpha = 0.1) +
theme_bw()
ggsave('diastolic-bp.png')
No. | Subject | Author | Date |
---|---|---|---|
674 | 안녕하세요 mytable 함수 관련하여 질문드립니다. [1] | 지옥소년 | 2019.07.26 |
» | 교수님, sphagetti plot으로 문의드립니다. [3] | cardio_H | 2019.07.22 |
672 | 교차분석 결과로 로지스틱 회귀분석 가능한지 질문드립니다 [1] | 김도 | 2019.07.22 |
671 | 그래프 출력 오류 [2] | 김날 | 2019.07.17 |
670 | mytable 함수 추가 질문 드립니다 [3] | 수수 | 2019.07.17 |
669 | 틀리지 않는법 | cardiomoon | 2019.07.13 |
668 | mytable에서 하는 통계법 질문 [1] | 수수 | 2019.07.11 |
667 | 여러 군의 평균 값 비교에 대해 질문드립니다. [2] | jera0131 | 2019.07.10 |
666 | mytable에서의 P value [2] | RR | 2019.07.08 |
665 | 정회원신청합니다. [1] | 박소이 | 2019.07.08 |
664 | mytable 및 ztable 관련 질문드립니다 [1] | jera0131 | 2019.07.05 |
663 | 정회원 서버접속이 안됩니다. [1] | 세희맘 | 2019.07.05 |
662 | 정회원 로그인이 안됩니다 [1] | jaljary | 2019.07.05 |
661 | 정회원 서버와 web r studio 접속이 안됩니다. [1] | jera0131 | 2019.07.05 |
660 | 정회원서버 접속 실패 [1] | cellory | 2019.07.05 |
659 | 정회원서버 [1] | 박소이 | 2019.07.02 |
658 | 정회원 신청드립니다. [3] | cellory | 2019.07.02 |
657 | survival analysis | diddbwls9 | 2019.07.01 |
656 | 정회원 서버 접속 실패 | Sally | 2019.06.28 |
655 | 후원금 입금 관련? [2] | duehdaud | 2019.06.25 |
Back-calculated age와 Back-calculated length가 뭐지요?