Reputation: 1
I am currently trying to implement the Prentice, Williams and Peterson extension of the Cox model in R on survival data but I want a parametric model and can not get it done. Is there maybe a package I am unaware of or has anyone else created a parametric PWP model before? I would greatly appreciate any help with this.
Thanks J
Upvotes: 0
Views: 283
Reputation: 1026
I don't see why you would use a parametric model.
This means specifying parametric forms for baseline hazards (and in PWP you have one for every event).
Regardless, you can use maximum likelihood for this, but you will probably have to program it yourself.
If you denote the hazard for the j-th event by lambda[0] (j), then the contribution of an individual i would be:
Just choose your parametric form for the lambdas, take the log, and plug into optim()
.
Upvotes: 1