coxph는 연도별로 계층화됩니다.

Question

나는 이것이 아주 쉽게해야한다고 생각하지만, 나는 그 주위에서 머리를 터뜨릴 수 없습니다.

나는 다음과 같은 코드가 있습니다

library(survival) 
cox <- coxph(Surv(SURV, DEAD)~YEAR, data) 
summary(cox) 

을하지만 아래로 개별 년에 결과 분할을하고 싶습니다.

COXREG surv /STATUS=dead(1) /CONTRAST (year)=Indicator(1) 
/METHOD=ENTER year /PRINT=CI(95) 
/CRITERIA=PIN(.05) POUT(.10) ITERATE(20). 
EXECUTE. 

및 STATA에서 같은 일 : 여기

는 SPSS 구문 및 솔루션의 모습 무엇을 여기

xi: stcox i.year 

의 출력입니다 ,515,

str(data) 

Answer 1

당신은 어떻게이 "데이터를"준 재현 예를 구성하는 우리에게 STR (데이터) 또는 표시되지 않았다 원하는 대부분의 문서 페이지 하단의 사례가있다. 나는 "YEAR"이 숫자 벡터로 판명 난다고 생각합니다. 요소 변수라면 요격과 n-1 계수를 보았을 것입니다. Interecpt 계수는 "연도"와 같고 다른 계수는 연도 (n) 값과 일치합니다. SPSS 엔진에 "연도"가 "표시기"이지만 R 엔진에 동일한 호의를 제공하지 않았다고 말했습니다.

이 시도 : '데이터'

data$year.ind <- factor(data$year) # equivalent of SPSS INDICATOR 
            # or SAS /CLASS 
cox.mdl <- coxph(Surv(SURV, DEAD)~YEAR, data) 
as.matrix(coef(coc.mdl) 
summary(cox.mdl) 

Answer 2

R은 종종 더 많은 자유를 허용하는 연산 결과 표시를 분할한다. coxph (?predict.coxph)의 예측 기능이 필요하다고 가정합니다.

는

predict(cox, type="terms")