안녕하세요, 항상 너무 많은 도움을 받고 있습니다.
생존분석 중 곡선을 그렸는데,
구한 데이터에서 time=0의 survival rate은 1(100%)로 나옵니다만, 이상하게도 생존곡선을 그려보면 시작점이 1(100%)이 아니어서 이를 문의드립니다.
[Survival rate by time] - 해당 결과값은 위의 그림 중 2번째 그래프에 해당되는 것입니다.
제가 time의 단위를 쪼개서(원래 year 이었는데 day로 쪼갬) 그래프 앞단을 확장해서 보니
다음과 같이 시작점 day 0의 값이 없고 day1부터 곡선이 시작됨을 확인했습니다. 혹시 이를 고칠 수 있는 방법이 있을까요?
Comment 1
-
cardiomoon
2022.08.20 10:07
말씀하신대로 문제가 있습니다. 그래서 제가 만든 autoReg 패키지의 adjustedPlot을 쓰시면 해결됩니다. 제가 쓴 "의료인을 위한 R생존분석"책을 참조하세요.
다음 코드를 실행해보셔요
library(survival)
fit <- coxph(Surv(time, status) ~ age + sex, lung)
library(autoReg)
adjustedPlot(fit,xnames=c("sex"))
adjustedPlot(fit,xnames=c("age","sex"),facet="sex")