R/contabs.netwBMA.R

contabs.netwBMA <-
function (network, reference, known=NULL, thresholds = NULL)
{

   prelim <- contabs.prelim( network, reference, known)

   if (is.null(thresholds)) {
     prob <- prelim$network[,3]
#    thresholds <- sort(unique(prob))
     thresholds <- sort(as.numeric(names(table(prob))))
   }

   if (any(thresholds > 1)) stop("thresholds: use probabilities rather than percentages")

   contabs( network = prelim$network, reference = prelim$reference, 
            size = prelim$size, thresholds = thresholds)
}

Try the networkBMA package in your browser

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

networkBMA documentation built on Jan. 28, 2021, 2:02 a.m.