R 코드 | 여러개의 위험비(HR)와 신뢰구간(CI)을 하나의 forest plot으로 그리기

아래와 같이 정리한 위험비(HR, 95% CI) 데이터 forestplotdata.csv가 있다. 위 데이터를 R의 forestplot 패키지를 이용하여 원하는 모양으로 직접 설정한 포레스트 플롯 forest plot을 아래와 같이 그려 보았다.

작성한 R 코드는 아래와 같다. library(forestplot) # 파일에서 데이터를 읽어온다 data <- read.csv("forestplotdata.csv", header = T) # 두 개의 새로운 문자열 벡터 np와 hr을 생성 # np: 이벤트와 전체 케이스 수, 비율을 결합하여 하나의 문자열로 # hr: 위험비(Hazard Ratio)와 95% 신뢰구간(CI)을 결합하여 하나의 문자열로 np <- paste(data$event, " / ", data$n," (",data$rate,")",sep="") hr <- paste(data$HR, " (", data$lower, "-", data$upper, ")",sep="") # ...