다중 상태 생존 분석하여 R 패키지 "생존"

Question

나는 다음과 같이 R survival 패키지를 사용하여 상태 1에서 상태 2로 전환을 위해 콕스 비례 위험 모델을 추정하기 위해 노력하고

:

내가없이 모델을 추정하고

Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + 
      wSTA, data=Multistate, subset = from == 1, 
      frailty(x=id, distribution= "gamma")) 

인수가 frailty인데, 정상적으로 작동합니다. 내가 잘못 가고 어디

Error in if (any(ord > 1)) stop("Penalty terms cannot be in an interaction") : 
missing value where TRUE/FALSE needed 

사람이 설명해주십시오 수 : 나는 그것을 포함하고 때, 나는 다음과 같은 오류는 무엇입니까?

감사합니다, 저도 같은 문제를했지만, 결국 연약함 기간이 첨가제 것을 깨달았다

Answer 1

J! 즉, 쉼표 뒤에 오는 것이 아니라 더하기 다음에 오는 것이 좋습니다.

Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + wSTA + frailty(x=id, distribution= "gamma"), data=Multistate, subset = from == 1)