From dc9b4f7f9cb44ab171eb3e6d54737eb88e3167ef Mon Sep 17 00:00:00 2001
From: Waldir Leoncio <waldir.leoncio@gmail.com>
Date: Thu, 19 Nov 2020 08:44:38 +0100
Subject: [PATCH] Translated initializeGammaln

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

diff --git a/R/initializeGammaln.R b/R/initializeGammaln.R
new file mode 100644
index 0000000..ac1864f
--- /dev/null
+++ b/R/initializeGammaln.R
@@ -0,0 +1,9 @@
+initializeGammaln <- function(ninds, rowsFromInd, maxAlleles) {
+	#Alustaa GAMMALN muuttujan s.e. GAMMALN(i, j)=gammaln((i - 1) + 1/j)
+	GAMMA_LN <- zeros((1 + ninds) * rowsFromInd, maxAlleles)
+	for (i in 1:(ninds + 1) * rowsFromInd) {
+		for (j in 1:maxAlleles) {
+			GAMMA_LN[i, j] <- log_gamma((i - 1) + 1/j)
+		}
+	}
+}