query=rep(NA,length(net.1)) names(query) <- names(net.1) query['lung'] <- 'true' evi=rep(NA,length(net.1)) names(evi) <- names(net.1) evi['dysp'] <- 'true' res <- my.cn.inference(net.1,net.2,query,evi) cat('Query p(',res$query,'|',res$evidence, ') -- interval result: [',res$min.p,',',res$max.p,']\n') evi['bronc'] <- 'true' res <- my.cn.inference(net.1,net.2,query,evi) cat('Query p(',res$query,'|',res$evidence, ') -- interval result: [',res$min.p,',',res$max.p,']\n') evi['bronc'] <- 'false' evi['xrays'] <- 'true' res <- my.cn.inference(net.1,net.2,query,evi) cat('Query p(',res$query,'|',res$evidence, ') -- interval result: [',res$min.p,',',res$max.p,']\n') evi['smoke'] <- 'true' res <- my.cn.inference(net.1,net.2,query,evi) cat('Query p(',res$query,'|',res$evidence, ') -- interval result: [',res$min.p,',',res$max.p,']\n')