library(bnlearn) source('my.bn.inference.r') net = model2network("[x1][x2|x1][x3|x1][x4|x2:x3]") cpt1.x1 = matrix(c(0.7, 0.3), ncol = 2, dimnames = list(NULL, c('true', 'false'))) cpt1.x2 = c(0.1, 0.9, 0.3, 0.7) dim(cpt1.x2) = c(2, 2) dimnames(cpt1.x2) = list("x2" = c("true", "false"), "x1" = c("true", "false")) cpt1.x3 = c(0.5, 0.5, 0.2, 0.8) dim(cpt1.x3) = c(2, 2) dimnames(cpt1.x3) = list("x3" = c("true", "false"), "x1" = c("true", "false")) cpt1.x4 = c(0.9, 0.1, 0.5, 0.5, 0.4, 0.6, 0.1, 0.9) dim(cpt1.x4) = c(2, 2, 2) dimnames(cpt1.x4) = list("x4" = c("true", "false"), "x2" = c("true", "false"), "x3" = c("true", "false"))