Hi,
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,
Hi Peter,
Apologies - I ran the script using S-Plus 6.2. I am happy to run it in R, but
my feeling is that there may be something wrong with the code itself. I have
included the dataset and script, in case you can help.
Jindi Singh.
Jatinder Singh
Senior Manager, Analysis and Reporting
PRA
Hi,
I am trying to convert the following NLMIXED code to NLME, but am
running into problems concerning 'Singularity in backsolve'. As I am new
to R/S-Plus, I thought I may be missing something in the NLME code.
NLMIXED
***
proc nlmixed data=kidney.kidney;
parms delta=0.03 gamma=1.1