install.packages("BiocManager") BiocManager::install(c("graph", "RBGL", "Rgraphviz"),force=TRUE) install.packages("gRim") library(gRim) train <- read.csv(file="train.csv",sep=",") train$Survived <- factor(train$Survived) train$Pclass <- ordered(train$Pclass) train$Embarked <- factor(train$Embarked,levels=c("C","Q","S")) train$Age <- cut(train$Age,breaks=c(0,10,20,30,45,60,90)) train$Fare <- cut(train$Fare,breaks=c(0,20,50,100,200,600)) train$Parch <- cut(train$Parch,breaks=c(-0.5,0.5,6)) train$SibSp <- cut(train$SibSp,breaks=c(-0.5,0.5,8)) n <- dim(train)[1] model.loglin <- dmod(~.^.,data=train[,c(2,3,5,6,7,8,10,12)]) iplot(model.loglin) model.loglin <- backward(model.loglin,crit="aic",k=log(n)) iplot(model.loglin) model.loglin <- dmod(~.^1,data=train[,c(2,3,5,6,7,8,10,12)]) iplot(model.loglin) model.loglin <- forward(model.loglin,crit="aic",k=log(n)) iplot(model.loglin)