Hi All, In absence of any reply, I am posting a slightly modified question. What is "x" in hazard.ratio in the command below?
example(hazard.ratio) binscores<-cut(scores.train,c(-1000,-1,1,1000),c("low","intermediate","high")) * hazard.ratio(x=?, surv.time=train$ProgFreeSurv, surv.event=train$PfSStatus, strat=binscores) * Input to "x" is ("x a vector of risk predictions"), what does it mean? I have computed scores (binscore) above from a 100-gene signature I have derived from an expression data. I need to compute hazard rations of the groups (low and high) defined by me using "cut" command above. Thanks, Angel. ------ I am interested in calculating hazard ratio of the stratified groups defined by me. library(survcomp) binscores <- cut(scores.train,c(-1000,-1,1,1000),c("low","intermediate","high")) dd <- data.frame("surv.time"=OS, "surv.event"= status, "strat"=binscores) km.coxph.plot(formula.s=Surv(surv.time, surv.event) ~ strat) How do I calculate hazard ratio between the groups defined by me as above? Thanks very much, Angel [[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.