demo/chest.R

In gRaven: Bayes Nets: 'RHugin' Emulation with 'gRain'

chest<-hugin.domain()

yn <- c("yes","no")
add.node(chest,"asia",states=yn)
add.node(chest,"tub",states=yn)
add.node(chest,"smoke",states=yn)
add.node(chest,"lung",states=yn)
add.node(chest,"bronc",states=yn)
add.node(chest,"either",states=yn)
add.node(chest,"xray",states=yn)
add.node(chest,"dysp",states=yn)

add.edge(chest,"tub","asia")
add.edge(chest,"lung","smoke")
add.edge(chest,"bronc","smoke")
add.edge(chest,"either",c("lung","tub"))
add.edge(chest,"xray","either")
add.edge(chest,"dysp",c("bronc","either"))

set.table(chest,"asia",c(0.01,0.99))
set.table(chest,"tub",c(0.05,0.95,0.01,0.99))
set.table(chest,"smoke",c(0.5,0.5))
set.table(chest,"lung",c(0.1,0.9,0.01,0.99))
set.table(chest,"bronc",c(0.6,0.4,0.3,0.7))
set.table(chest,"either",c(1,0,1,0,1,0,0,1))
set.table(chest,"xray",c(0.98,0.02,0.05,0.95))
set.table(chest,"dysp",c(0.9,0.1,0.7,0.3,0.8,0.2,0.1,0.9))

chest

compile(chest)

summary(chest,jt=TRUE)

get.normalization.constant(chest)

set.finding(chest,"smoke","no")
set.finding(chest,"dysp","yes")

get.normalization.constant(chest)

get.belief(chest,'lung')
get.belief(chest,'bronc')
map.configurations(chest,get.nodes(chest),0.001)

chest2<-clone.domain(chest)
retract(chest2)
summary(chest2)

simulate(chest2,10)

list.domains()