diff options
| author | Thibaut Horel <thibaut.horel@gmail.com> | 2015-09-14 23:08:02 -0400 |
|---|---|---|
| committer | Thibaut Horel <thibaut.horel@gmail.com> | 2015-09-14 23:08:02 -0400 |
| commit | ab0b1f3cefedb35327a19ec1b6afd560bfdf802d (patch) | |
| tree | b777f3e2c0ac0e712d8c5faab5107b1d236e2c3a /experiments/process.py | |
| parent | 960676226862d2d68c7a9c04c56d4f8157803025 (diff) | |
| download | criminal_cascades-ab0b1f3cefedb35327a19ec1b6afd560bfdf802d.tar.gz | |
Import supplements and repo reorganization
Diffstat (limited to 'experiments/process.py')
| -rw-r--r-- | experiments/process.py | 63 |
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() - - - |
