Esempio n. 1
0
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",
        )
Esempio n. 2
0
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')
Esempio n. 3
0
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')