Nothing
### confint.influenceTest.R ---
##----------------------------------------------------------------------
## Author: Brice Ozenne
## Created: jun 1 2018 (12:15)
## Version:
## Last-Updated: Dec 21 2021 (12:29)
## By: Brice Ozenne
## Update #: 37
##----------------------------------------------------------------------
##
### Commentary:
##
### Change Log:
##----------------------------------------------------------------------
##
### Code:
## * confint.influenceTest (documentation)
##' @title Confidence Intervals and Confidence Bands for the Difference Between Two Estimates
##' @description Confidence intervals and confidence Bands for the difference between two estimates.
##' @name confint.influenceTest
##'
##' @param object A \code{influenceTest} object, i.e. output of the \code{influenceTest} function.
##' @param parm not used. For compatibility with the generic method.
##' @param level [numeric, 0-1] Level of confidence.
##' @param transform [character] the transformation used to improve coverage of the confidence intervals.
##' Can be \code{"none"} or \code{"atanh"}.
##' @param n.sim [integer, >0] the number of simulations used to compute the quantiles for the confidence bands.
##' @param seed [integer, >0] seed number set before performing simulations for the confidence bands.
##' If not given or NA no seed is set.
##' @param ... not used.
##'
##' @details Except for the cumulative hazard,
##' the confidence bands and confidence intervals are automatically restricted to the interval [-1;1].
##'
##' @author Brice Ozenne
##'
## * confint.influenceTest (code)
##' @rdname confint.influenceTest
##' @method confint influenceTest
##' @export
confint.influenceTest <- function(object,
parm = NULL,
level = 0.95,
n.sim = 1e4,
transform = "none",
seed = NA,
...){
## ** check arguments
dots <- list(...)
if(length(dots)>0){
txt <- names(dots)
txt.s <- if(length(txt)>1){"s"}else{""}
stop("unknown argument",txt.s,": \"",paste0(txt,collapse="\" \""),"\" \n")
}
if(object$type == "cumhazard"){
object$transform <- match.arg(transform, c("none"))
}else{
object$transform <- match.arg(transform, c("none","atanh"))
}
## ** compute se, CI/CB
outCIBP <- transformCIBP(estimate = object$delta,
se = object$delta.se,
iid = object$delta.iid,
null = 0,
conf.level = level,
n.sim = n.sim,
seed = seed,
type = object$transform,
min.value = switch(object$transform,
"none" = -1,
"atanh" = NULL),
max.value = switch(object$transform,
"none" = 1,
"atanh" = NULL),
ci = TRUE,
band = object$band,
method.band = "maxT-simulation",
alternative = "two.sided",
p.value = TRUE)
names(outCIBP) <- paste0("delta.", names(outCIBP))
object[names(outCIBP)] <- outCIBP
## export
object$conf.level <- level
return(object)
}
######################################################################
### confint.influenceTest.R ends here
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.