From 8abd467ca64f849ae6c19e7659ed8465ec47f782 Mon Sep 17 00:00:00 2001
From: Waldir Leoncio <waldir.leoncio@gmail.com>
Date: Mon, 19 Oct 2020 14:08:25 +0200
Subject: [PATCH] Added computePopulationLogml.R

---
 R/computePopulationLogml.R | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)
 create mode 100644 R/computePopulationLogml.R

diff --git a/R/computePopulationLogml.R b/R/computePopulationLogml.R
new file mode 100644
index 0000000..bf9782d
--- /dev/null
+++ b/R/computePopulationLogml.R
@@ -0,0 +1,24 @@
+computePopulationLogml <- function(pops, adjprior, priorTerm) {
+	# Palauttaa length(pops)*1 taulukon, jossa on laskettu korikohtaiset
+
+	x <- size(COUNTS, 1)
+	y <- size(COUNTS, 2)
+	z <- length(pops)
+
+	popLogml <- squeeze(
+		sum(
+			sum(
+				reshape(
+					lgamma(
+						repmat(adjprior, c(1, 1, length(pops))) +
+							COUNTS[, , pops]
+					),
+					c(x, y, z)
+				),
+				1
+			),
+			2
+		)
+	) - sum(lgamma(1 + SUMCOUNTS[pops, ]), 2) - priorTerm
+	return(popLogml)
+}
\ No newline at end of file