...
 
Commits (2)
......@@ -25,11 +25,13 @@ def sim(mu1, mu2, pi1, pi2):
return p025s, p975s
with open("results/publication-figure-2-b.tsv", "w", encoding="utf-8") as fig2:
fig2.write("π₁\tL-RI\tL-RI-CI5\tL-RI-CI95\n")
for pi1 in np.arange(0.0, 0.51, 0.01):
result = sim(10, 18, pi1, 0)
fig2.write(str(pi1) + "\t" + str(np.percentile(result[0], 50)) + "\t" + str(np.percentile(result[0], 5)) + "\t" + str(np.percentile(result[0], 95)) + "\n")
with open("results/publication-figure-2-c.tsv", "w", encoding="utf-8") as fig2:
fig2.write("µ₁\tL-RI\tL-RI-CI5\tL-RI-CI95\n")
for mu1 in np.arange(8.0, 14.1, 0.1):
result = sim(mu1, 18, 0.20, 0)
fig2.write(str(mu1) + "\t" + str(np.percentile(result[0], 50)) + "\t" + str(np.percentile(result[0], 5)) + "\t" + str(np.percentile(result[0], 95)) + "\n")