from dfg4grn import *
settings = ReadConfig(settings, "./config/default_values/dfg4grn.cfg")
settings = ReadConfig(settings, settings["dfg4grn"]["config"])
grid = Generate_Grid("dfg4grn", None, settings, ["eta_z", "lambda_w", "tau"], 25).test_list
jobman = JobManager(settings)
for p in grid:
settings["dfg4grn"]["eta_z"] = p[0]
settings["dfg4grn"]["lambda_w"] = p[1]
settings["dfg4grn"]["tau"] = p[2]
dfg = DFG4GRN(
timeseries_storage, trans_factors, settings, "dfg4grn-" + str(p[0]) + "." + str(p[1]) + "." + str(p[2])
)
jobman.queueJob(dfg)
print jobman.queue
jobman.runQueue()
jobman.waitToClear()
accs = []
for job in jobman.finished:
print job.alg.gene_list
print job.alg.gather_output(settings)
jobnet = Network()
jobnet.read_netmatrix(job.alg.network, job.alg.gene_list, True)
accs.append(jobnet.calcAcc(goldnet))
print accs
