# load the data and build the correct network from the model string.
data(insurance)
res = empty.graph(names(insurance))
modelstring(res) = paste("[Age][Mileage][SocioEcon|Age][GoodStudent|Age:SocioEcon]",
"[RiskAversion|Age:SocioEcon][OtherCar|SocioEcon]",
"[VehicleYear|SocioEcon:RiskAversion][MakeModel|SocioEcon:RiskAversion]",
"[SeniorTrain|Age:RiskAversion][HomeBase|SocioEcon:RiskAversion]",
"[AntiTheft|SocioEcon:RiskAversion][RuggedAuto|VehicleYear:MakeModel]",
"[Antilock|VehicleYear:MakeModel][DrivingSkill|Age:SeniorTrain]",
"[CarValue|VehicleYear:MakeModel:Mileage][Airbag|VehicleYear:MakeModel]",
"[DrivQuality|RiskAversion:DrivingSkill][Theft|CarValue:HomeBase:AntiTheft]",
"[Cushioning|RuggedAuto:Airbag][DrivHist|RiskAversion:DrivingSkill]",
"[Accident|DrivQuality:Mileage:Antilock][ThisCarDam|RuggedAuto:Accident]",
"[OtherCarCost|RuggedAuto:Accident][MedCost|Age:Accident:Cushioning]",
"[ILiCost|Accident][ThisCarCost|ThisCarDam:Theft:CarValue]",
"[PropCost|ThisCarCost:OtherCarCost]", sep = "")
# there are too many nodes for plot(), use graphviz.plot().
graphviz.plot(res)
Run the code above in your browser using DataLab