https://github.com/cran/ensembleBMA
Raw File
Tip revision: a685fedf9ca178c112c82c8b14539132a04b88b6 authored by Chris Fraley on 15 December 2009, 00:00:00 UTC
version 4.5
Tip revision: a685fed
quantBMAnormal.R
quantBMAnormal <-
function(alpha, WEIGHTS, MEAN, SD)
{
  lower <- min(MEAN-6*SD)
  upper <- max(MEAN+6*SD)

  if (cdfBMAnormal(lower, WEIGHTS, MEAN, SD, 0) > alpha) return(NA)
  if (cdfBMAnormal(upper, WEIGHTS, MEAN, SD, 0) < alpha) return(NA)

  z <- uniroot(cdfBMAnormal, lower = lower, upper = upper,
               WEIGHTS=WEIGHTS, MEAN=MEAN, SD=SD, offset = alpha)

  z$root
}

back to top