https://github.com/cran/quantreg
Raw File
Tip revision: 61e5d2e598fb5e53345fa9414a615a4d7b524808 authored by Roger Koenker on 18 October 2020, 10:50:03 UTC
version 5.74
Tip revision: 61e5d2e
KMvCRQ.R
# Example Comparison of Kaplan-Meier vs crq fitting
# The red crq estimate should overplot the black KM Survival Curve.
if (requireNamespace("survival", quietly = TRUE)){
    n <- 100
    y <- rchisq(n,3)
    c <- rchisq(n,5)
    Y <- pmin(y,c)
    d <- (y < c)
    Surv <- survival::Surv
    plot(survival::survfit(Surv(Y,d)~1))
    f <- crq(Surv(Y,d)~1, method = "Portnoy", grid = "pivot")
    x <- f$sol[2,]
    p <- 1-f$sol[1,]
    p <- c(p,p[length(p)])
    par(col = "red")
    fs <- plot(stepfun(x, p),do.points = FALSE, add = TRUE)
}
back to top