R help - Jan 2010 - no "rcorrp.cens" in hmisc package

Dear,

I wanna to compare AUC generated by two distribution models using the same
sample.

I tried improveProb function's example code below.
set.seed(1)
library(survival)
x1 <- rnorm(400)
x2 <- x1 + rnorm(400)
d.time <- rexp(400) + (x1 - min(x1))
cens <- runif(400,.5,2)
death <- d.time <= cens
d.time <- pmin(d.time, cens)
rcorrp.cens(x1, x2, Surv(d.time, death))

#rcorrp.cens(x1, x2, y) ## no censoring
set.seed(1)
x1 <- runif(1000)
x2 <- runif(1000)
y <- sample(0:1, 1000, TRUE)
rcorrp.cens(x1, x2, y)
improveProb(x1, x2, y)

 However, a error message jumped out, saying "rcorrp.cens" not found.

Please kindly help and thank you.

Elaine

	[[alternative HTML version deleted]]

You omitted library(Hmisc) or require(Hmisc)
Frank

elaine kuo wrote:
>  Dear,
> 
> I wanna to compare AUC generated by two distribution models using the same
> sample.
> 
> I tried improveProb function's example code below.
> set.seed(1)
> library(survival)
> x1 <- rnorm(400)
> x2 <- x1 + rnorm(400)
> d.time <- rexp(400) + (x1 - min(x1))
> cens <- runif(400,.5,2)
> death <- d.time <= cens
> d.time <- pmin(d.time, cens)
> rcorrp.cens(x1, x2, Surv(d.time, death))
> 
> #rcorrp.cens(x1, x2, y) ## no censoring
> set.seed(1)
> x1 <- runif(1000)
> x2 <- runif(1000)
> y <- sample(0:1, 1000, TRUE)
> rcorrp.cens(x1, x2, y)
> improveProb(x1, x2, y)
> 
>  However, a error message jumped out, saying "rcorrp.cens" not
found.
> 
> Please kindly help and thank you.
> 
> Elaine
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.
> 

-- 
Frank E Harrell Jr   Professor and Chairman        School of Medicine
                      Department of Biostatistics   Vanderbilt University