summaryrefslogtreecommitdiffstats
path: root/experiments/process.py
diff options
context:
space:
mode:
authorThibaut Horel <thibaut.horel@gmail.com>2015-09-14 23:08:02 -0400
committerThibaut Horel <thibaut.horel@gmail.com>2015-09-14 23:08:02 -0400
commitab0b1f3cefedb35327a19ec1b6afd560bfdf802d (patch)
treeb777f3e2c0ac0e712d8c5faab5107b1d236e2c3a /experiments/process.py
parent960676226862d2d68c7a9c04c56d4f8157803025 (diff)
downloadcriminal_cascades-ab0b1f3cefedb35327a19ec1b6afd560bfdf802d.tar.gz
Import supplements and repo reorganization
Diffstat (limited to 'experiments/process.py')
-rw-r--r--experiments/process.py63
1 files changed, 0 insertions, 63 deletions
diff --git a/experiments/process.py b/experiments/process.py
deleted file mode 100644
index f339fd1..0000000
--- a/experiments/process.py
+++ /dev/null
@@ -1,63 +0,0 @@
-import sys
-from ml import ml
-from ml2 import ml2
-from ml3 import ml3
-import numpy as np
-from cPickle import load
-from itertools import product
-from math import exp
-import time
-
-
-if __name__ == "__main__":
- time.sleep(3)
- if len(sys.argv) < 2:
- sys.exit("usage: {0} <file>".format(sys.argv[0]))
- root_victims, victims, non_victims, age = load(open(sys.argv[1]))
-
- # alphas = np.arange(1, 200, 10) # parameter of the time component
- alphas = np.logspace(-6,1,num=5)
- # deltas = np.arange(.00003, .0003, .00002) # parameter of the structural component
- deltas = np.logspace(-6,-1,num=5)
- lmbdas = np.logspace(-8,-3,num=5)
- # lmbda = 1
- with open("out.log", "w") as fh:
- for alpha, delta, lmbda in product(alphas, deltas, lmbdas):
- # for alpha, delta in product(alphas, deltas):
- beta, roots, ll = ml(root_victims, victims, non_victims, alpha, delta, lmbda)
- print "\t".join(map(str, [alpha, delta, beta, roots, ll]))
- fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n")
- fh.flush()
-
- # alphas = np.arange(.15, .25, .01) # parameter of the time component
- # # alphas = np.logspace(-4,-.5,num=15)
- # deltas = np.arange(.08, .1, .01) # parameter of the structural component
- # # deltas = np.logspace(-2,-.1,num=15)
- # lmbda = 0.10#np.logspace(-12,-3,num=10)
- # with open("out.log", "w") as fh:
- # for alpha, delta in product(alphas, deltas):
- # lmbda, roots, ll = ml2(root_victims, victims, non_victims, alpha, delta, lmbda)
- # print "\t".join(map(str, [alpha, delta, lmbda, roots, ll]))
- # fh.write("\t".join(map(str, [alpha, delta, lmbda, roots, ll])) + "\n")
- # fh.flush()
-
- # alpha = 0.0373
- # delta = 0.00518
- # beta, roots, ll = ml(root_victims, victims, non_victims, alpha, delta, 1.)
- # print "\t".join(map(str, [alpha, delta, beta, roots, ll]))
-
-# from ml import ml
-# alphas = np.arange(.01, .05, .01) # parameter of the time component
-# # alphas = np.logspace(-4,1,num=5)
-# deltas = np.arange(.001, .02, .003) # parameter of the structural component
-# # deltas = np.logspace(-10,-1,num=8)
-# lmbda = 1.
-# with open("out.log", "w") as fh:
-# for alpha, delta in product(alphas, deltas):
-# beta, roots, ll = ml(root_victims, victims, non_victims, alpha, delta, lmbda)
-# print "\t".join(map(str, [alpha, delta, beta, roots, ll]))
-# fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n")
-# fh.flush()
-
-
-