From ba96cb38b2d34dae5ef93902a70c582a369ec1a0 Mon Sep 17 00:00:00 2001
From: Gregory Ashton <gregory.ashton@ligo.org>
Date: Tue, 30 Jun 2020 10:38:16 +1000
Subject: [PATCH] Fix bug in evidence combination
---
bilby/core/result.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/bilby/core/result.py b/bilby/core/result.py
index e756aaf75..75df60a37 100644
--- a/bilby/core/result.py
+++ b/bilby/core/result.py
@@ -1650,7 +1650,7 @@ class ResultList(list):
result_weights = np.exp(log_evidences - np.max(log_evidences))
log_errs = [res.log_evidence_err for res in self if np.isfinite(res.log_evidence_err)]
if len(log_errs) > 0:
- result.log_evidence_err = logsumexp(2 * np.array(log_errs), b=1. / len(self))
+ result.log_evidence_err = 0.5 * logsumexp(2 * np.array(log_errs), b=1. / len(self))
else:
result.log_evidence_err = np.nan
posteriors = list()
--
GitLab