diff options
Diffstat (limited to 'experiments/process.py')
| -rw-r--r-- | experiments/process.py | 10 |
1 files changed, 6 insertions, 4 deletions
diff --git a/experiments/process.py b/experiments/process.py index ffa3162..f339fd1 100644 --- a/experiments/process.py +++ b/experiments/process.py @@ -16,12 +16,14 @@ if __name__ == "__main__": 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(-3,1,num=16) + alphas = np.logspace(-6,1,num=5) # deltas = np.arange(.00003, .0003, .00002) # parameter of the structural component - deltas = np.logspace(-4,-1,num=12) - lmbda = 1. + 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 in product(alphas, deltas): + 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") |