From 034c71328460d460deef74bf591c8d6cd28d5f93 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Wed, 8 Jul 2015 00:41:29 -0400 Subject: updated temporal model for different time scales, ran some tests --- experiments/process.py | 13 +++++++------ 1 file changed, 7 insertions(+), 6 deletions(-) (limited to 'experiments/process.py') diff --git a/experiments/process.py b/experiments/process.py index 6c71934..33c83b6 100644 --- a/experiments/process.py +++ b/experiments/process.py @@ -16,12 +16,13 @@ if __name__ == "__main__": 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(-1,-.5,num=10) + alphas = np.logspace(-2.5,-1,num=6) # deltas = np.arange(.02, .1, .01) # parameter of the structural component - deltas = np.logspace(-2,-1,num=10) + deltas = np.logspace(-3,-2,num=6) + 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, age, alpha, delta) + 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() @@ -38,7 +39,7 @@ if __name__ == "__main__": # fh.write("\t".join(map(str, [alpha, delta, lmbda, roots, ll])) + "\n") # fh.flush() - # alpha = .016 - # delta = 0.077 - # beta, roots, ll = ml2(root_victims, victims, non_victims, alpha, delta, 1.) + # 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 -- cgit v1.2.3-70-g09d2