Nothing
FctPersMonNonMixExp <- function(data, lambda.hat, group.name) {
## Calculates the values of some function of the person months
## in the non-mixture model with exponential survival, i. e.
## S(t) = c^[1 - exp(- lambda * t)], lambda > 0, 0 < c < 1, t >= 0.
##
## Args:
## data: Data frame which consists of at least three columns for the group
## in the first (all of the same group),
## status (1 = event, 0 = censored) in the second
## and event time in the third column.
## lambda.hat: Estimator of the parameter lambda.
## group.name: Name of the group.
##
## Results:
## Returns the value of some function of the person months.
# function of person months, and their verification
o.stroke <- sum(1 - exp(- lambda.hat * data[, 3]))
if (o.stroke <= 0) {
stop(paste("Number of person months in", group.name, "must be bigger than 0.",
call. = FALSE))
}
return(o.stroke)
}
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.