Re: [R] Conditional gap time frailty cox model for recurrent events

2016-09-06 Thread Therneau, Terry M., Ph.D.
You can ignore the message below. The maximizing routine buried within the frailty() command buried with coxph() has a maximizer that is not the brightest. It sometimes gets lost but then finds its way again. The message is from one of those. It likely took a not-so-good update step, and too

Re: [R] conditional gap time frailty cox model for recurrent events

2016-09-05 Thread Göran Broström
Dear Elisabetta, I have no direct answer to your question, but a suggestion: Use the 'coxme' function (in the package with the same name). In the help page for 'frailty' (survival) you will find: "The coxme package has superseded this method. It is faster, more stable, and more flexible." Ht

[R] conditional gap time frailty cox model for recurrent events

2016-09-05 Thread Elisabetta Petracci
Dear users, I am fitting a conditional gap time frailty cox model weighting observations by means of inverse probability time dependent weigths. Attached find the self-explaining dataset. I have used the following sintax: coxph(Surv(gaptstart,gaptstop,status)~treat+strata(nrecord01)+frailty(id,d