[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