Nothing
# -------------------------------------------------------------------------------------------
#' Computes confidence interval of Inverse Gamma distributed variable with given mean and
#' standard deviation.
#'
#' @param mu A \code{k x 1} vector of means.
#' @param sd A \code{k x 1} vector of standard deviations.
#' @param qlower A \code{k x 1} vector of lower quantiles.
#' @param qupper A \code{k x 1} vector of upper quantiles.
#'
#' @return A \code{2 x k} matrix containing the lower and upper bounds of the intervals.
#'
#' @importFrom stats qgamma
#' @keywords internal
.intervalIGamma <- function(mu, sd, qlower = 0.025, qupper = 1 - qlower) {
# convert to shape and scale parameters alpha and beta
alpha <- 2 + mu^2 / sd^2
beta <- mu * (alpha - 1)
# compute confidence interval
int <- matrix(c(
1 / qgamma(1 - qlower, shape = alpha, rate = beta),
1 / qgamma(1 - qupper, shape = alpha, rate = beta)
),
2, length(mu),
byrow = TRUE
)
if (int[1] > mu | int[2] < mu) {
warning("Variance constraint initialization: The supplied mean does not belong to the chosen interval of the Inverse Gamma distribution. \n")
}
return(int)
}
# ---------------------------------------------------------------------------------------------------
#' Capitalizes the first letter of a string.
#'
#' @param x A character.
#' @keywords internal
firstLetterUp <- function(x) {
substr(x, 1, 1) <- toupper(substr(x, 1, 1))
return(x)
}
# ---------------------------------------------------------------------------------------------------
#' Capitalizes the first letter of a string.
#'
#' @param a A multiple time series object.
#' @param a A matrix.
#' @keywords internal
matmult3d <- function(a,b) {
n <- NROW(a)
m <- NCOL(a)
k <- dim(b)[3]
if (k > 1) {
y <- do.call(rbind, lapply(1:n, function(x) a[x, ] %*% t(b[, , x])))
y <- ts(y, start = start(a), frequency = frequency(a))
} else {
y <- a %*% t(b[,,1])
}
y <- ts(y, start = start(a), frequency = frequency(a))
return(y)
}
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.