inst/extdata/ordered_gaussian.R

In robertness/bninfo: Queries and information theoritic operations on Bayesian networks.

library(clusterGeneration)
# A big tree-like graph --------------------------
library(clusterGeneration)
# A MVN Network Simulated with --------------------------
num.vars <- 40
num.obs <- 1000
## Mean 0 data from a random MVN covariance matrix
net <- simGaussianNet(num.vars, c(.5, 1.5), "ordered")
## Data simulated directly from the Bayesian network.  
boot <- rbn(hairy, n = num.obs) %>% 
  boot.strength(R = 500, m = nrow(hairy.data),
                algorithm = "hc", 
                algorithm.args = list(score = "bic-g"))
ordered.boot <- list(truth = net, boot = boot)
save(ordered.boot, file = "data/ordered_boot.Rdata")