https://github.com/cran/ensembleBMA
Tip revision: 65bc4d8d72e5060987d37deb8279042bbcd033ff authored by Chris Fraley on 08 May 2008, 00:00:00 UTC
version 3.0-2
version 3.0-2
Tip revision: 65bc4d8
quantBMAnormal.R
`quantBMAnormal` <-
function(alpha, WEIGHTS, MEAN, SD)
{
lower <- 0
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
}