evidence=rep(NA,length(net.1))
names(evidence) <- names(net.1)
evidence[4]='true'
res <- my.bn.inference(net.1,query,evidence)
##res not shown
evidence[4]='false'
my.bn.inference(net.1,query,evidence)
## $logp.evi
## [1] -0.3816998
## $logp.query
## [1] -0.8187104
## $p
## [1] 0.6459646
## $evidence
## [1] "x4=false"
## $query
## [1] "x1=true,x2=false"