From dfc9659180f73a3f5efe15f3365ce181818e94d8 Mon Sep 17 00:00:00 2001
From: Waldir Leoncio <waldir.leoncio@gmail.com>
Date: Tue, 14 Jul 2020 15:33:00 +0200
Subject: [PATCH] translated laskeLoggis()

---
 R/greedyMix.R   | 12 +-----------
 R/laskeLoggis.R |  9 +++++++++
 2 files changed, 10 insertions(+), 11 deletions(-)
 create mode 100644 R/laskeLoggis.R

diff --git a/R/greedyMix.R b/R/greedyMix.R
index 6624a4d..4aa71fc 100644
--- a/R/greedyMix.R
+++ b/R/greedyMix.R
@@ -978,14 +978,4 @@ greedyMix <- function(
 # 	else
 # 		popDiff = diff([0 pops npops+1]);
 # 		emptyPop = min(find(popDiff > 1));
-# 	end
-
-# 	%------------------------------------------------------
-
-# 	function loggis = laskeLoggis(counts, sumcounts, adjprior)
-# 	npops = size(counts,3);
-
-# 	logml2 = sum(sum(sum(gammaln(counts+repmat(adjprior,[1 1 npops]))))) ...
-# 		- npops*sum(sum(gammaln(adjprior))) - ...
-# 		sum(sum(gammaln(1+sumcounts)));
-# 	loggis = logml2;
\ No newline at end of file
+# 	end
\ No newline at end of file
diff --git a/R/laskeLoggis.R b/R/laskeLoggis.R
new file mode 100644
index 0000000..7bcd26f
--- /dev/null
+++ b/R/laskeLoggis.R
@@ -0,0 +1,9 @@
+laskeLoggis <- function(counts, sumcounts, adjprior) {
+	npops <- size(counts, 3)
+
+	sum1 <- sum(sum(sum(gammaln(counts + repmat(adjprior, c(1, 1, npops))))))
+	sum3 <- sum(sum(lgamma(adjprior))) - sum(sum(lgamma(1 + sumcounts)))
+	logml2 <- sum1 - npops * sum3
+	loggis <- logml2
+	return(loggis)
+}
\ No newline at end of file