library(bnlearn)

options(echo=TRUE)
args <- commandArgs(trailingOnly = TRUE)
print(args)

data <- read.table(args[1], sep = ",", header = TRUE)
col_names <- names(data)
data[, col_names] <- lapply(data[,col_names] , factor)

start = random.graph(nodes = names(data), num = 100)
netlist = lapply(start, function(net) {
  hc(data, score = "bic", iss = 1, start = net) })

save.image(file = "/home/otu/hdogan/bnp_R/bnp_R_image.gzip", compress=TRUE)
q()