An example of bnlearn
library(bnlearn) norm <- rnorm(4000) Data <- matrix(norm, nrow=1000, ncol=4, byrow=T) colnames(Data) <- c("Height", "BMI", "SBP", "FBS") Data <- [Link](Data) Data2 <- Data Data2$Height <- 170 + Data$Height*10 Data2$SBP <- 120 + Data$SBP*10 Data2$BMI <- 22 + Data$Height*2 + Data$BMI*2 + Data$SBP*2 Data2$FBS <- 90 + (Data2$BMI-22)*5 + Data$FBS*10 Data2 <- [Link](Data2) res.Data2 = gs(Data2) plot(res.Data2) res.arc2 <- [Link](res.Data2, "FBS","BMI") [Link](res.arc2, Data2)