[R] Frailty with a Weibull PH Model

From: Singh, Jatinder <SinghJatinder_at_PRAIntl.com>
Date: Wed 14 Dec 2005 - 05:13:51 EST


I am trying to fit a Weibull PH model with a gaussian frailty term (kidney data). Is the following approach correct?

  1. Fit the model using survReg with the frailty term added in.

Kidn1<-survreg(formula = Surv(rtime, event)~ age+sex+gn+an+pkn+frailty.gaussian(patient), data=kidney, dist='weibull')

2. Based on Venables and Ripley (350-353), divide the negative of the coefficients by the scale parameter to obtain the Weibull PH model coefficients.


Jatinder Singh
Senior Manager, Analysis and Reporting
PRA International
300-730 View Street
Victoria, B.C. V8W 3Y7
Tel: 250-483-4416
Fax: 250 483 4588
e-mail: singhjatinder@praintl.com

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Dec 14 07:27:47 2005

This archive was generated by hypermail 2.1.8 : Wed 14 Dec 2005 - 09:34:06 EST