data(kidney)
# type 'help(kidney)' for a description of the data set
kidney$sex <- kidney$sex - 1
mod <- parfm(Surv(time,status) ~ sex + age, cluster = "id",
data = kidney, dist = "exponential", frailty = "gamma")
tau(mod)
Run the code above in your browser using DataLab