Could someone explain the summary(cph.object)?
The example is in the help file of cph.
n <- 1000
set.seed(731)
age <- 50 + 12*rnorm(n)
label(age) <- "Age"
sex <- factor(sample(c('Male','Female'), n,
rep=TRUE, prob=c(.6, .4)))
cens <- 15*runif(n)
h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
dt <- -log(runif(n))/h
label(dt) <- 'Follow-up Time'
e <- ifelse(dt <= cens,1,0)
dt <- pmin(dt, cens)
units(dt) <- "Year"
dd <- datadist(age, sex)
options(datadist='dd')
Srv <- Surv(dt,e)
f <- cph(Srv ~ rcs(age,4) + sex, x=TRUE, y=TRUE)
summary(f)
Effects Response : Srv
Factor Low High Diff. Effect S.E. Lower 0.95 Upper 0.95
age 40.872 57.385 16.513 1.21 0.21 0.80 1.62
Hazard Ratio 40.872 57.385 16.513 3.35 NA 2.22 5.06
sex - Female:Male 2.000 1.000 NA 0.64 0.15 0.35 0.94
Hazard Ratio 2.000 1.000 NA 1.91 NA 1.42 2.55
Wat's the meaning of Effect, S.E. Lower, Upper?
Thanks
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
