From b26412ed5a3e08e9dc32fc73c27c42be54d82aa8 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Sat, 20 Jun 2015 19:09:29 -0400 Subject: altered weight_success to show dist and dt between each parent and child --- experiments/process.py | 28 ++++++++++++++-------------- 1 file changed, 14 insertions(+), 14 deletions(-) (limited to 'experiments/process.py') diff --git a/experiments/process.py b/experiments/process.py index 330a49e..49fe7a0 100644 --- a/experiments/process.py +++ b/experiments/process.py @@ -12,18 +12,18 @@ if __name__ == "__main__": sys.exit("usage: {0} ".format(sys.argv[0])) root_victims, victims, non_victims, age = load(open(sys.argv[1])) - # alphas = np.arange(5e-4, 3e-3, 1e-4) # parameter of the time component - alphas = np.logspace(-5,-1,num=20) - # deltas = np.arange(0.03, 0.1, 0.005) # parameter of the structural component - deltas = np.logspace(-5,-1,num=20) - with open("out.log", "w") as fh: - for alpha, delta in product(alphas, deltas): - beta, roots, ll = ml2(root_victims, victims, non_victims, age, alpha, delta) - 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(5e-4, 3e-3, 1e-4) # parameter of the time component + # alphas = np.logspace(-5,-1,num=20) + # # deltas = np.arange(0.03, 0.1, 0.005) # parameter of the structural component + # deltas = np.logspace(-5,-1,num=20) + # with open("out.log", "w") as fh: + # for alpha, delta in product(alphas, deltas): + # beta, roots, ll = ml2(root_victims, victims, non_victims, age, alpha, delta) + # print "\t".join(map(str, [alpha, delta, beta, roots, ll])) + # fh.write("\t".join(map(str, [alpha, delta, beta, roots, ll])) + "\n") + # fh.flush() - # alpha = 1/10. - # delta = .5 - # beta, roots, ll = ml2(root_victims, victims, non_victims, age, alpha, delta) - # print "\t".join(map(str, [1./alpha, delta, beta, roots, ll])) \ No newline at end of file + alpha = 0.004 + delta = 0.06 + beta, roots, ll = ml2(root_victims, victims, non_victims, age, alpha, delta) + print "\t".join(map(str, [1./alpha, delta, beta, roots, ll])) \ No newline at end of file -- cgit v1.2.3-70-g09d2