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)
+		}
+	}
+}