if (FALSE) {
data(hdfail)
# Select only Western Digital hard drives
dat <- subset(hdfail, grepl("WDC", model))
fit.hd <- fitfrail(Surv(time, status) ~ temp + rer + rsc
+ psc + cluster(model),
dat, frailty="gamma", fitmethod="score")
fit.hd
}
Run the code above in your browser using DataLab