From 4cb6f238fca09a817d88b2a8414faa1608c5dfd2 Mon Sep 17 00:00:00 2001 From: Rick Sexton Date: Mon, 22 Jan 2024 19:36:47 -0700 Subject: [PATCH] corrected lnp --- basicrta/functions.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/basicrta/functions.py b/basicrta/functions.py index 9ade366..b8e3be3 100644 --- a/basicrta/functions.py +++ b/basicrta/functions.py @@ -166,7 +166,7 @@ def run(self): lnp[j] = np.log(tmp.take(s)).sum()+\ np.log(mcweights[j][uniqs]).sum()-\ (mcrates[j][uniqs]*rhypers[uniqs, 1]).sum()+\ - np.log(mcweights[j]**(whypers-1)).sum() + np.log(mcweights[j][uniqs]**(whypers[uniqs]-1)).sum() Ns[j][:] = np.array([len(inds[i]) for i in range(ncomp)]) Ts = np.array([x[inds[i]].sum() for i in range(ncomp)])