## Graphical model from Dawid and Evett (1997) ## http://www.math.leidenuniv.nl/~gillrd/teaching/graphical/DawidEvett.net ## http://www.math.leidenuniv.nl/~gillrd/teaching/graphical/DawidEvettSN.net ## "SN" stands for "short names", better for some plots library(Rgraphviz) ## From bioconductor repository library(gRain) ## From CRAN repository network <- loadHuginNet("DawidEvett.net") ## Created by me with GeNIe networkSN <- loadHuginNet("DawidEvettSN.net") plot(network) network names(network) plot(networkSN) networkSN names(networkSN) networkC <- compile(network) networkC names(networkC) networkCP <- propagate(networkC) # plot(networkCP) networkCP names(networkCP) names(network$cptlist) network$cptlist[10] networkCP querygrain(networkCP,"CSuspectOneOfOffenders") networkCPev <-setFinding(networkCP,"Y1TypeOfFibreTuft","Special") networkCPev <-setFinding(networkCPev,"Y2BloodtypeSprayOnJumper","Special") networkCPev <-setFinding(networkCPev,"X1SuspectsBloodtype","Other") networkCPev <-setFinding(networkCPev,"X2GuardsBloodtype","Special") networkCPev <-setFinding(networkCPev,"X3TypeOfSuspectsJumper","Special") networkCPev <-setFinding(networkCPev,"NumberOffenders","Num_4") networkCPev <-setFinding(networkCPev,"WArrestingPoliceOfficer","Suspect_was_offender") networkCPev <-setFinding(networkCPev,"Guard1HowManyOffenders","Says_4") networkCPev <-setFinding(networkCPev,"Guard2WasPunchedOnTheNose","Yes") querygrain(networkCPev,"CSuspectOneOfOffenders")