Re: [R] about the summary(cph.object)

2009-08-02 Thread zhu yao
Thanks for your help! 2009/8/2 David Winsemius > Here's what I would do. Let's assume that you are presenting the results of > the example on the cph help page. I agree with you that the results should > be presented on the hazard ratio scale. The Design package provides > appropriate plotting

Re: [R] about the summary(cph.object)

2009-08-02 Thread David Winsemius
Here's what I would do. Let's assume that you are presenting the results of the example on the cph help page. I agree with you that the results should be presented on the hazard ratio scale. The Design package provides appropriate plotting tools for creation of publication quality graphics

Re: [R] about the summary(cph.object)

2009-08-01 Thread Frank E Harrell Jr
zhu yao wrote: Thx for your reply. In this example, age was transformed with rcs. So the output was different between f and summary(f). If I need to publicate the results, how do I explation the hazard ratio of age? David explained this. Nonlinearity in age does not complicate the explanation

Re: [R] about the summary(cph.object)

2009-08-01 Thread zhu yao
Thx for your reply. In this example, age was transformed with rcs. So the output was different between f and summary(f). If I need to publicate the results, how do I explation the hazard ratio of age? 2009/8/1 David Winsemius > > On Jul 31, 2009, at 11:24 PM, zhu yao wrote: > > Could someone ex

Re: [R] about the summary(cph.object)

2009-08-01 Thread David Winsemius
On Jul 31, 2009, at 11:24 PM, zhu yao wrote: 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 <- 1

[R] about the summary(cph.object)

2009-07-31 Thread zhu yao
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=='F