# NOT RUN {
yn <- c("yes","no")
a <- cptable(~asia, values=c(1,99), levels=yn)
t.a <- cptable(~tub + asia, values=c(5,95,1,99), levels=yn)
plist <- compileCPT(list(a, t.a ))
bn <- grain(plist)
bnc <- compile(bn, propagate=FALSE)
bncp <- compile(bn, propagate=TRUE)
## New p(tub | asia)
z <- c(20, 80, 1, 99)
bn2 <- setCPT(bn, list(tub=z))
bnc2 <- setCPT(bnc, list(tub=z))
bncp2 <- setCPT(bncp, list(tub=z))
# }
Run the code above in your browser using DataLab