Nothing
#
# Ingmar Visser, 8-11-2019
#
# Description
#
# Computes confidence intervals for (dep)mix model parameters.
#
# Details
#
# Confidence intervals are computed through the variance-covariance matrix
# which in turn is computed using the hessian and the linear constraints
# of the model. See ?vcov, ?hessian and ?standarderror for more details on
# these underlying functions. The confidence intervals are computed using
# the normal theory assumptions. The desired significance level can be
# supplied through the 'level' argument.
#
# Value
#
#
setMethod("confint", "mix",
function(object, level=0.95, fixed=NULL, equal=NULL,
conrows=NULL, conrows.upper=NULL, conrows.lower=NULL,
tolerance=1e-6,
method="finiteDifferences", ...) {
vc <- vcov(object,fixed=fixed,equal=equal,
conrows=conrows,conrows.upper=conrows.upper,conrows.lower=conrows.lower,
tolerance=tolerance,method=method, ...)
ses <- sqrt(diag(vc$vcov))
pars <- getpars(object)
elements <- vc$elements
parsinc <- pars[which(elements=="inc")]
upper <- parsinc+qnorm(0.5+level/2)*ses
lower <- parsinc-qnorm(0.5+level/2)*ses
ret <- data.frame(par=pars, constr=elements, lower=NA, upper=NA)
ret$lower[which(elements=="inc")] <- lower
ret$upper[which(elements=="inc")] <- upper
colnames(ret)[3:4] <- c(formatperc(0.5-level/2, 3),formatperc(0.5+level/2, 3))
return(ret)
}
)
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.