def compute_bayes():
fa.save_log_evidence(b, edep=True)
fa.save_log_evidence(b_disk, edep=True)
fa.save_log_evidence(b_poiss, edep=True)
fa.save_log_evidence(b_just_NFW, edep=True)
fa.save_log_evidence(b_iso, edep=True)
if plus_IPS:
fa.save_log_evidence(b_plus_IPS, edep=True)
print "D&N Vs. D"
cbc = compute_bayes_class(
b.plots_dir_for_run + "lge.txt", b_disk.plots_dir_for_run + "lge.txt", b.plots_dir_for_run + "bayes.txt"
)
print "D&N Vs. Nothing"
cbc_poiss = compute_bayes_class(
b.plots_dir_for_run + "lge.txt", b_poiss.plots_dir_for_run + "lge.txt", b_poiss.plots_dir_for_run + "bayes.txt"
)
print "D&N Vs. I"
cbc_I = compute_bayes_class(
b.plots_dir_for_run + "lge.txt", b_iso.plots_dir_for_run + "lge.txt", b_iso.plots_dir_for_run + "bayes.txt"
)
print "D&N Vs. N"
cbc_just_nfw = compute_bayes_class(
b.plots_dir_for_run + "lge.txt",
b_just_NFW.plots_dir_for_run + "lge.txt",
b_just_NFW.plots_dir_for_run + "bayes.txt",
)
if plus_IPS:
print "D&N&I Vs. D&N"
cbc_plus_IPS = compute_bayes_class(
b_plus_IPS.plots_dir_for_run + "lge.txt",
b.plots_dir_for_run + "lge.txt",
b_plus_IPS.plots_dir_for_run + "bayes.txt",
)
def compute_bayes():
if just_disk:
fa.save_log_evidence(b_disk,edep=True)
if poiss:
fa.save_log_evidence(b_poiss,edep=True)
if plus_IPS:
fa.save_log_evidence(b_plus_IPS,edep=True)
if plus_IPS and just_disk:
cbc_plus_IPS = compute_bayes_class(b_plus_IPS.plots_dir_for_run + 'lge.txt',b_disk.plots_dir_for_run + 'lge.txt', b_plus_IPS.plots_dir_for_run + 'bayes.txt')
def compute_bayes():
fa.save_log_evidence(b,edep=True)
fa.save_log_evidence(b_disk,edep=True)
fa.save_log_evidence(b_poiss,edep=True)
fa.save_log_evidence(b_just_NFW,edep=True)
fa.save_log_evidence(b_iso,edep=True)
if plus_IPS:
fa.save_log_evidence(b_plus_IPS,edep=True)
print 'D&N Vs. D'
cbc = compute_bayes_class(b.plots_dir_for_run + 'lge.txt',b_disk.plots_dir_for_run + 'lge.txt', b.plots_dir_for_run + 'bayes.txt')
print 'D&N Vs. Nothing'
cbc_poiss = compute_bayes_class(b.plots_dir_for_run + 'lge.txt',b_poiss.plots_dir_for_run + 'lge.txt', b_poiss.plots_dir_for_run + 'bayes.txt')
print 'D&N Vs. I'
cbc_I = compute_bayes_class(b.plots_dir_for_run + 'lge.txt',b_iso.plots_dir_for_run + 'lge.txt', b_iso.plots_dir_for_run + 'bayes.txt')
print 'D&N Vs. N'
cbc_just_nfw = compute_bayes_class(b.plots_dir_for_run + 'lge.txt',b_just_NFW.plots_dir_for_run + 'lge.txt', b_just_NFW.plots_dir_for_run + 'bayes.txt')
if plus_IPS:
print 'D&N&I Vs. D&N'
cbc_plus_IPS = compute_bayes_class(b_plus_IPS.plots_dir_for_run + 'lge.txt',b.plots_dir_for_run + 'lge.txt', b_plus_IPS.plots_dir_for_run + 'bayes.txt')