R/probDistance.r

Defines functions probDistance.COPResult probDistance

probDistance <- function(result, numSimulations = BLCOPOptions("numSimulations")) 
{
    monteCarloSample <- rmnorm(numSimulations, result@priorMean, result@priorCovar)
    mean(abs(dmnorm(monteCarloSample, result@priorMean, result@priorCovar,log=TRUE) - dmnorm(monteCarloSample, result@posteriorMean, 
        result@posteriorCovar, log = TRUE)))
}

setGeneric("probDistance")

probDistance.COPResult <- function(result, numSimulations = BLCOPOptions("numSimulations") )
{
    show("Not implemented yet...")
}

setMethod("probDistance", signature(result = "COPResult"), probDistance.COPResult)

Try the BLCOP package in your browser

Any scripts or data that you put into this service are public.

BLCOP documentation built on Jan. 26, 2021, 1:05 a.m.