diff options
Diffstat (limited to 'experiments/process.py')
| -rw-r--r-- | experiments/process.py | 26 |
1 files changed, 21 insertions, 5 deletions
diff --git a/experiments/process.py b/experiments/process.py index 33c83b6..ffa3162 100644 --- a/experiments/process.py +++ b/experiments/process.py @@ -15,10 +15,10 @@ if __name__ == "__main__": sys.exit("usage: {0} <file>".format(sys.argv[0])) root_victims, victims, non_victims, age = load(open(sys.argv[1])) - # alphas = np.arange(.02, .1, .01) # parameter of the time component - alphas = np.logspace(-2.5,-1,num=6) - # deltas = np.arange(.02, .1, .01) # parameter of the structural component - deltas = np.logspace(-3,-2,num=6) + # alphas = np.arange(1, 200, 10) # parameter of the time component + alphas = np.logspace(-3,1,num=16) + # deltas = np.arange(.00003, .0003, .00002) # parameter of the structural component + deltas = np.logspace(-4,-1,num=12) lmbda = 1. with open("out.log", "w") as fh: for alpha, delta in product(alphas, deltas): @@ -42,4 +42,20 @@ if __name__ == "__main__": # 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]))
\ No newline at end of file + # 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() + + + |