coxme 오브젝트 (R)로 진행 절차

Question

R에 Cox 혼합 모델을 구축해야하는 심장 이식 데이터 작업 중입니다. 이 목적으로 설계되었으므로 coxme 패키지를 사용하기로 선택했지만 모델을 선택하는 데 많은 어려움이 있습니다. 전진 절차를 사용하고 싶지만 coxme 개체에는 적용 할 수 없습니다. 나는 대답을 step()의 코드에 접근하여 수정할 수 있도록 여기에 또 다른 질문을했다. 따라서 coxme 모델에 적용 할 수있다. 그러나 이것은 매우 어렵습니다 (적어도 저에게는) 기능이 매우 복잡합니다 (이것은 단지 exctractAIC()을 수정하는 것만은 아닙니다. 적응할 다른 많은 기능이 있습니다).

1) coxme 모델에 기대를 실행할 수있는 간단한 방법이 있나요 :

그래서 여기 내 질문입니까?

2) 다른 패키지를 사용해야합니까?

Answer 1

step()을 수정하는 데 오랜 시간이 걸린 후에 처음부터 간단하게 정확한 전진 함수를 작성하기로 결정했습니다. 그래서 누군가 나처럼 똑같은 질문을한다면, 정말로 그를 권합니다.

1) 처음부터 전체 기능을 코딩하는 것을 두려워하지 않아야합니다! (원래 함수의 80 %는 포워드 프로 시저에 유용하지 않으며 다시 코딩되지 않습니다.)

2) coxme 패키지를 사용합니다.

그러나 혼란을 막을 수있는 한 가지는 다음과 같습니다. coxme 모델에서 AIC/BIC라고하는 것은 사실 null 모델과 비교하여 테스트 된 모델의 개선입니다. 즉, 최대화해야하며 사용을 결정하면 최소화하지 않아야합니다 (원하는 경우 실제 AIC/BIC로 계산을 변경할 수 있음).