Nothing
confband <- function(object, newdata, level = .95, ...)
UseMethod("confband")
confband.mlt <- function(object, newdata, level = 0.95,
type = c("trafo", "distribution", "survivor", "cumhazard"),
K = 20, cheat = K, ...) {
stopifnot(!missing(newdata))
stopifnot(is.data.frame(newdata))
type <- match.arg(type)
if (nrow(newdata) > 1) {
ret <- lapply(1:nrow(newdata), function(i)
confband(object = object, newdata = newdata[i,,drop = FALSE],
level = level, type = type, K = K, cheat = K, ...))
return(ret)
}
stopifnot(nrow(newdata) == 1)
y <- variable.names(object, "response")
q <- mkgrid(object, n = K)[[y]]
nd <- newdata[rep(1, length(q)),,drop = FALSE]
nd[[y]] <- q
X <- model.matrix(object$model$model, data = nd)
ci <- confint(multcomp::glht(multcomp::parm(coef(object), vcov(object)),
linfct = X), ...)$confint
### use quantile obtained for K contrasts for larger number
### of contrasts
if (cheat > K)
return(confband(object = object, newdata = newdata, level = level,
type = type, K = cheat, calpha = attr(ci, "calpha")))
if (type == "distribution") ci <- object$model$todistr$p(ci)
if (type == "survivor") ci <- 1 - object$model$todistr$p(ci)
if (type == "cumhazard") ci <- -log(1 - object$model$todistr$p(ci))
ci <- cbind(q = q, ci)
return(ci)
}
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.