From 4b0448d05efdd0b1522fddca3e2e3fcbed5d5f3b Mon Sep 17 00:00:00 2001 From: Waldir Leoncio Date: Fri, 27 Sep 2024 12:59:55 +0200 Subject: [PATCH] Fixed calculation of logml --- R/laskeLoggis.R | 6 ++---- 1 file changed, 2 insertions(+), 4 deletions(-) diff --git a/R/laskeLoggis.R b/R/laskeLoggis.R index bb6ced9..52eae31 100644 --- a/R/laskeLoggis.R +++ b/R/laskeLoggis.R @@ -2,8 +2,6 @@ laskeLoggis <- function(counts, sumcounts, adjprior) { npops <- size(counts, 3) replicated_adjprior <- array(adjprior, c(nrow(adjprior), ncol(adjprior), npops)) sum1 <- sum(sum(sum(lgamma(counts + replicated_adjprior)))) - sum3 <- sum(sum(lgamma(adjprior))) - sum(sum(lgamma(1 + sumcounts))) - logml2 <- sum1 - npops * sum3 - loggis <- logml2 - return(loggis) + sum2 <- npops * sum(sum(lgamma(adjprior))) + sum(sum(lgamma(1 + sumcounts))) + sum1 - sum2 }