2017-02-06 19 views
0

survival 패키지와 함께 생존 분석을 수행하는 방법에 대한 정보 (인터넷 및 교과서)가 상당히 있습니다. 그러나 왼쪽에 데이터가 생략되었을 때이 작업을 수행하는 방법에 대한 정보는 찾지 못했습니다.간격 검열 데이터 : 콕스 비례 위험 및 R에서 서바이벌 차이


문제의 배경 :

내가 게시 생존 데이터와 설정 자체 제작 데이터를 가지고있다. 일반적으로 이벤트 시간과 마지막 후속 조치 (오른쪽 검열) 날짜가 주어집니다. 그러나 그 사건이 360 일 이전에 일어났다 고 한 연구 만 있습니다. 그래서 나는 왼쪽으로 검열 된이 데이터를 남겼습니다.

내가하고 싶은 것 :

내가 왼쪽 절단, 이벤트 및 우측 절단 세트 전체 데이터를 분석하고자합니다. 나는성에 의하여 카플란 - 마이어 곡선을 플롯하려면 다음

  1. 는 콕스 회귀 분석을 수행 로그 순위 테스트를

내가 필요한 것 :

내가 type = interval2으로 Surv 객체를 만들 수 있습니다. 그러나 이것은 생존 패키지의 survdiff이나 coxph을 계산하지 못합니다.

intcox 패키지가 CRAN에서 제거되었으며 icenReg 또는 interval 패키지에서 검색 한 것을 찾지 못했습니다.


누구든지 내 문제를 해결하는 방법이나 이에 대한 실제 정보를 어디에서 찾을 수 있습니까? 나는 이미이 일을 쓰고있다.

감사합니다.

+0

스택 오버플로의 통계 페이지에서 더 나은 대답을 얻을 수 있다고 생각합니다. 여기에서보십시오 : http://stats.stackexchange.com/ – sconfluentus

+0

Inter2 형 검열은 Therneau가 일반적으로 "전 시간"예측을 제공하는 것을 피하는 "시간 의존 공변량"과 동일한 문제를 만듭니다. 공변량 시간의 기간에 걸쳐 일정한 것으로 가정 될 수 없으므로 두 특정 경우에 대해 현명한 "diff"가있을 수 없습니다. 항상 값이있는 두 경우에 대해 잘 정의 된 공변량 집합을 지정하려면 해당 공변량을 나타내는 R 코드를 게시하십시오.(실용적인 정보는 Rhelp Archives에 Therneau의 게시물이 될 것입니다.) –

답변

0

ic_sp 기능을 사용하여 오른쪽 및 왼쪽 검열 용 Cox-PH 모델을 icenReg에 맞출 수 있습니다. 이 즉,

fit <- ic_sp(Surv(L, R, type = 'interval2') ~ treatment, data = myData) 

또는 더 간결

fit <- ic_sp(cbind(L, R) ~ treatment, data = myData) 

로그 순위 테스트와 작은 icenReg에서 사용할 수 없습니다를 표준 Surv 응답 변수를 사용하여 들어갈 수 있지만, interval 패키지에서 찾을 수 있습니다 .