0byt3m1n1
Path:
C:
/
Program Files
/
R
/
R-4.4.1
/
library
/
survival
/
tests
/
[
Home
]
File: fr_rat1.R
options(na.action=na.exclude) # preserve missings options(contrasts=c('contr.treatment', 'contr.poly')) #ensure constrast type library(survival) # Tests using the rats data # # (Female rats, from Mantel et al, Cancer Research 37, # 3863-3868, November 77) rfit <- coxph(Surv(time,status) ~ rx + frailty(litter), rats, method='breslow', subset= (sex=='f')) rfit rfit$iter rfit$df rfit$history[[1]] rfit1 <- coxph(Surv(time,status) ~ rx + frailty(litter, theta=1), rats, method='breslow', subset=(sex=="f")) rfit1 rfit2 <- coxph(Surv(time,status) ~ frailty(litter), rats, subset=(sex=='f')) rfit2