[R] Frailty with a Weibull PH Model

Singh, Jatinder SinghJatinder at PRAIntl.com
Tue Dec 13 19:13:51 CET 2005


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,

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


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 at praintl.com

More information about the R-help mailing list