별도의 생존 곡선 인덱싱

Question

ggplot의 두 그룹 각각에 대해 Kaplan-Meier 생존 추정치를 플로팅하고 싶습니다.

이렇게하려면 각 그룹에 대해 별도의 생존 곡선을 얻어야합니다. survival 패키지의 survfit 함수는 멋지게 나뉘지만 각 함수에 대해 별도의 플롯을 인덱싱하는 방법을 알지 못합니다. I가 맞는 경우 예

sCurve$n.event 

들어

rearrest<-read.table("http://stats.idre.ucla.edu/stat/examples/alda/rearrest.csv", sep=",", header=T) 

이

(sCurve <- summary(arr1 <- survfit(Surv(months, abs(censor-1))~1, data = rearrest))) 

그것은이 내의 인덱스 요소 쉽게 그룹화 곡선이다 : 여기

샘플 데이터이고 이 시간 이외의 동일한 항목은 personal 변수의 값에 따라 그룹화됩니다. 나는 두 개의 멋진 생존 곡선 객체를 준비합니다.

(sCurveA <- summary(arr1 <- survfit(Surv(months, abs(censor-1))~personal, data = rearrest))) 

하나의 객체는 personal=0 다른 personal=1를 표시합니다. $, [], [[]]으로 색인을 생성하려고 시도했지만 번호 유형 색인과 named-는 모두 사용할 수 없습니다.

아무도 도와 줄 수 있습니까?

Answer 1

sCurveA $ strata는 그룹화 변수를 벡터로 제공합니다. 키 조각을 꺼내 ggplot의 data.frame에 던져 넣을 수 있습니다.

df = data.frame(Time = sCurveA$time, 
       Survival = sCurveA$surv, 
       Strata = sCurveA$strata) 

ggplot(df, aes(Time, Survival, col = Strata)) + 
    geom_line()