# calculate the mean of the beta density given coefficients and testvalues
#' mean_beta
#'
#' @param betaobj betamle object
#' @param testvals vector
#' @param brange vector
#'
#' @return vector
#' @export
#'
mean_beta <- function(betaobj, testvals, brange=c(0,1)) {
nu <- betaobj$coef_nu
omega <- betaobj$coef_omega
xb_nu <- nu %*% testvals
xb_omega <- omega %*% testvals
rval <- exp(xb_nu)/(exp(xb_nu) + exp(xb_omega))
rval*(brange[2]-brange[1])+brange[1]
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.