[R] no "rcorrp.cens" in hmisc package

Frank E Harrell Jr f.harrell at vanderbilt.edu
Mon Jan 4 13:56:50 CET 2010


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



More information about the R-help mailing list