R/frontend-averaging.R

In bnlearn: Bayesian Network Structure Learning, Parameter Learning and Inference

# build the averaged network structure using arc strengths and a
# significance threshold.
averaged.network = function(strength, threshold) {

  # check the main argument.
  check.bn.strength(strength, valid = c("bootstrap", "bayes-factor"))
  # check the strength threshold.
  threshold = check.threshold(threshold, strength)

  avg = averaged.network.backend(strength = strength, threshold = threshold)

  # add the metadata for the print() method.
  avg$learning$algo = "averaged"
  avg$learning$args = list(threshold = threshold)

  return(avg)

}#AVERAGED.NETWORK

# compute the inclusion threshold from a set of arc strengths.
inclusion.threshold = function(strength) {

  # check the (only) argument.
  check.bn.strength(strength, valid = c("bootstrap", "bayes-factor"))

  threshold(strength = strength, method = "l1")

}#INCLUSION.THRESHOLD