R의 이산 시간 위험 모델 (loglog)

survival의 패키지 R가 나타납니다 연속 시간 생존 모델에 초점을 맞 춥니 다. 보완 로그-로그 모델 인 비례 위험 모델의 이산 시간 버전을 추정하는 데 관심이 있습니다. 간단한 오른쪽 검열을 통해 상당히 간단한 생존 모델이 있습니다.

이 모델을 추정하는 한 가지 방법은 "죽지 않은"각 기간 동안 각 관측치에 대해 별도의 행이있는 데이터 세트를 작성하는 것입니다. 그런 다음 링크가 있는 glm모델을 cloglog사용할 수 있습니다.

이 방법은 메모리가 매우 비효율적입니다. 실제로 내 컴퓨터의 메모리에 비해 너무 큰 데이터 세트를 생성 할 수 있습니다.

두 번째 방법은 MLE을 직접 코딩하는 것입니다. 그것은 충분히 간단 할 것이지만,이 생존 모델을 통조림으로 한 패키지가 있기를 바라고 있습니다. 공동 작업이 쉬워지고 패키지를 사용하는 코딩 오류를 피하는 것이 더 쉬울 것입니다.

누구든지 그러한 패키지를 알고 있습니까?

각 관측치에 대해 여러 개의 행을 갖는 것은 중복되어 보일 수 있지만 그렇지 않을 수도 있습니다. 모형에 시변 공변량이있는 경우 각 관측치-월에는 반드시 고유 한 행이 필요합니다. 시변 공변량의 하나의 특정 예는 경과 시간이다. 이 변수는 모델에 거의 확실하게 포함되어야하므로 각 관측 기간에 대해 별도의 행을 갖는 것이 좋습니다. 따라서 제안 된 첫 번째 방법은 아마도 가장 좋은 방법 일 것입니다.

이는와 이블 분포가있는 연속 시간 비례 위험 모델과 다릅니다. 여기서 경과 시간이 유일한 시변 공변량 인 경우 생존 모델을 각 관측치에 대해 한 줄로 단순화 할 수 있습니다 ( 예 : 여기 참조 ). Cox 비례 위험 모델에 대해서도 유사한 결과가 나타납니다.