###Install the package (bnlearn) install.packages("bnlearn") library(bnlearn) ###Load the data into your workspace d=read.csv("simulatedbundleddata.csv") d=d[,-1] for(i in 1:ncol(d)){ d[,i]=as.factor(d[,i]) } ###Paste the following functions in gssequence = function(Matrix){ print("For each of the following determine it's timing") sort=as.matrix(colnames(Matrix)) sort=cbind(sort,sort) for( i in 1:nrow(sort)){ x=sprintf ("Enter time for %s,",sort[i,1]) sort[i,2]=readline(x) } blacklist=matrix(nrow=1,ncol=2) colnames(blacklist)=c("from","to") w=as.matrix(as.numeric(sort[,2])) m=max(w) for ( i in 1:nrow(sort)) { if (w[i,1] > 1) { for ( kk in 1:nrow(sort)){ if ( w[kk,1] < w[i,1]){ blacklist=rbind(blacklist,c(sort[i,1],sort[kk,1])) } }}} blacklist=na.omit(blacklist) t=gs(d,blacklist=blacklist) plot(t) return(t) } hitonsequence = function(Matrix){ print("For each of the following determine it's timing") sort=as.matrix(colnames(Matrix)) sort=cbind(sort,sort) for( i in 1:nrow(sort)){ x=sprintf ("Enter time for %s,",sort[i,1]) sort[i,2]=readline(x) } blacklist=matrix(nrow=1,ncol=2) colnames(blacklist)=c("from","to") w=as.matrix(as.numeric(sort[,2])) m=max(w) for ( i in 1:nrow(sort)) { if (w[i,1] > 1) { for ( kk in 1:nrow(sort)){ if ( w[kk,1] < w[i,1]){ blacklist=rbind(blacklist,c(sort[i,1],sort[kk,1])) } }}} blacklist=na.omit(blacklist) t=si.hiton.pc(d,blacklist=blacklist) plot(t) return(t) } iambsequence = function(Matrix){ print("For each of the following determine it's timing") sort=as.matrix(colnames(Matrix)) sort=cbind(sort,sort) for( i in 1:nrow(sort)){ x=sprintf ("Enter time for %s,",sort[i,1]) sort[i,2]=readline(x) } blacklist=matrix(nrow=1,ncol=2) colnames(blacklist)=c("from","to") w=as.matrix(as.numeric(sort[,2])) m=max(w) for ( i in 1:nrow(sort)) { if (w[i,1] > 1) { for ( kk in 1:nrow(sort)){ if ( w[kk,1] < w[i,1]){ blacklist=rbind(blacklist,c(sort[i,1],sort[kk,1])) } }}} blacklist=na.omit(blacklist) t=iamb(d,blacklist=blacklist) plot(t) return(t) } ###Now run the functions on the data ##For grow-shrink use l=gssequence(d) ##For incremental association use l=iambsequence(d) ##For semi-interleaved HITON use l=hitonsequence(d)