Nothing
confint.DirichletRegModel <- function(
object
, parm
, level = .95
, ...
, type = c("all", "beta", "gamma")
, exp = FALSE
){
e <- exp
type <- match.arg(type)
if(any(level <= 0) | any(level >= 1)) stop("level must be in (0, 1)")
level <- sort(level)
repar <- if(object$parametrization == "alternative") TRUE else FALSE
if(!repar) if(type == "gamma") warning("type = gamma not meaningful for the common parametrization - ignored")
co <- object$coefficients
se <- object$se
names(co) <- object$coefnames
names(se) <- object$coefnames
res <- list(
level = level
, type = type
, coefficients = coef(object)
, se = se
, e = e
, repar = repar
)
rci <- lapply(level, function(L){
cbind(qnorm( (1 - L)/2, co, se),
qnorm(L + (1 - L)/2, co, se))
})
res$ci <- lapply(seq_along(rci), function(i) list())
if(repar){
Xc <- rev(rev(cumsum(c(1, object$n.vars)))[-1])
Zc <- c(1, ncol(object$Z)) + rev(Xc)[1] - 1
for(ll in seq_along(rci)){
inti <- 0
for(i in seq_len(object$dims)){
res$ci[[ll]][i] <- if(i == object$base){
list(NULL)
} else {
inti <- inti + 1
list(rci[[ll]][Xc[inti]:(Xc[inti+1]-1), , drop = FALSE])
}
}
res$ci[[ll]][[length(res$ci[[ll]]) + 1]] <- rci[[ll]][Zc[1]:Zc[2] , , drop = FALSE]
names(res$ci[[ll]]) <- c(object$varnames, "gamma")
}
} else {
Xc <- cumsum(c(1, object$n.vars))
for(ll in seq_along(rci)){
inti <- 0L
for(i in 1:object$dims){
inti <- inti + 1L
res$ci[[ll]][[i]] <- rci[[ll]][Xc[inti]:(Xc[inti+1]-1) , , drop = FALSE]
}
}
}
if(e){
for(ll in 1:length(rci)){
for(these in which(!unlist(lapply(res$ci[[1]], is.null)))){
res$ci[[ll]][[these]] <- exp(res$ci[[ll]][[these]])
}
}
}
class(res) <- "DirichletRegConfint"
return(res)
}
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.