Nothing
#' Kaplan-Meier product-limit estimate of survival
#'
#' @description This function provides survival estimates using the product-limit
#' Kaplan-Meier estimator.
#'
#' @usage KM(time, status, t)
#'
#' @param time Survival time of the process.
#' @param status Censoring indicator of the survival time of the process; 0 if
#' the survival time is censored and 1 otherwise.
#' @param t The time for obtaining survival estimates.
#' @return Vector with Kaplan-Meier estimate of survival.
#' @references E. Kaplan and P. Meier. Nonparametric estimation from incomplete
#' observations. Journal of the American Statistical Association, 53:457-481,
#' 1958.
#' @examples
#' require(survival)
#' data("bladder4state")
#'
#' obj<- multidf(gap1=bladder4state$y1, event1=bladder4state$d1,
#' gap2=bladder4state$y2, status=bladder4state$d2,
#' size=bladder4state$size)
#'
#' obj2<-obj[[1]]
#' KM(time = obj2$time, status = obj2$status, t = 20)
#' fit <- survfit(Surv(obj2$time, obj2$status) ~ 1, data = obj2)
#' summary(fit, time = 20)$surv
#'
#' @author Gustavo Soutinho and Luis Meira-Machado
KM <- function(
time,
status,
t
) {
if (missing(time)) stop("Argument 'time' is missing with no default");
if (missing(status)) stop("Argument 'status' is missing with no default");
if (missing(t)) stop("Argument 't' is missing with no default");
return(
.C(
"KaplanMeierValueSort", as.double(time), as.integer(status),
as.integer( length(time) ), as.double(t), p=as.double(1),
PACKAGE="survivalREC"
)$p
);
} # KM
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.