diff --git a/src/analysis/core.jl b/src/analysis/core.jl
index e5f6db7f..9aef83b4 100644
--- a/src/analysis/core.jl
+++ b/src/analysis/core.jl
@@ -82,13 +82,19 @@ end
 
 
 """
-    trainGigaSOM(som::Som, dInfo::LoadedDataInfo;
-                 kernelFun::Function = gaussianKernel,
-                 metric = Euclidean(),
-                 somDistFun = distMatrix(Chebyshev()),
-                 knnTreeFun = BruteTree,
-                 rStart = 0.0, rFinal=0.1, radiusFun=expRadius(-5.0),
-                 epochs = 20)
+    trainGigaSOM(
+        som::Som,
+        dInfo::LoadedDataInfo;
+        kernelFun::Function = gaussianKernel,
+        metric = Euclidean(),
+        somDistFun = distMatrix(Chebyshev()),
+        knnTreeFun = BruteTree,
+        rStart = 0.0,
+        rFinal = 0.1,
+        radiusFun = expRadius(-5.0),
+        epochs = 20,
+        eachEpoch = (e, r, codes) -> nothing,
+    )
 
 # Arguments:
 - `som`: object of type Som with an initialised som
@@ -102,6 +108,9 @@ end
 - `rFinal`: target radius at the last epoch, defaults to 0.1
 - `radiusFun`: Function that generates radius decay, e.g. `linearRadius` or `expRadius(10.0)`
 - `epochs`: number of SOM training iterations (default 10)
+- `eachEpoch`: a function to call back after each epoch, accepting arguments
+  `(epochNumber, radius, codes)`. For simplicity, this gets additionally called
+  once before the first epoch, with `epochNumber` set to zero.
 """
 function trainGigaSOM(
     som::Som,
@@ -114,6 +123,7 @@ function trainGigaSOM(
     rFinal = 0.1,
     radiusFun = expRadius(-5.0),
     epochs = 20,
+    eachEpoch = (e, r, codes) -> nothing,
 )
 
     # set the default radius
@@ -127,6 +137,8 @@ function trainGigaSOM(
 
     codes = som.codes
 
+    eachEpoch(0, rStart, codes)
+
     for j = 1:epochs
         @debug "Epoch $j..."
 
@@ -145,6 +157,8 @@ function trainGigaSOM(
 
         wEpoch = kernelFun(dm, r)
         codes = (wEpoch * numerator) ./ (wEpoch * denominator)
+
+        eachEpoch(epoch, r, codes)
     end
 
     som.codes = copy(codes)
@@ -154,46 +168,19 @@ end
 
 """
     trainGigaSOM(som::Som, train;
-                 kernelFun::Function = gaussianKernel,
-                 metric = Euclidean(),
-                 somDistFun = distMatrix(Chebyshev()),
-                 knnTreeFun = BruteTree,
-                 rStart = 0.0, rFinal=0.1, radiusFun=expRadius(-5.0),
-                 epochs = 20)
+                 kwargs...)
 
 Overload of `trainGigaSOM` for simple DataFrames and matrices. This slices the
-data using `DistributedArrays`, sends them the workers, and runs normal
-`trainGigaSOM`. Data is `undistribute`d after the computation.
+data, distributes them to the workers, and runs normal `trainGigaSOM`. Data is
+`undistribute`d after the computation.
 """
-function trainGigaSOM(
-    som::Som,
-    train;
-    kernelFun::Function = gaussianKernel,
-    metric = Euclidean(),
-    somDistFun = distMatrix(Chebyshev()),
-    knnTreeFun = BruteTree,
-    rStart = 0.0,
-    rFinal = 0.1,
-    radiusFun = expRadius(-5.0),
-    epochs = 20,
-)
+function trainGigaSOM(som::Som, train; kwargs...)
 
     train = Matrix{Float64}(train)
 
     #this slices the data into parts and and sends them to workers
     dInfo = distribute_array(:GigaSOMtrainDataVar, train, workers())
-    som_res = trainGigaSOM(
-        som,
-        dInfo,
-        kernelFun = kernelFun,
-        metric = metric,
-        somDistFun = somDistFun,
-        knnTreeFun = knnTreeFun,
-        rStart = rStart,
-        rFinal = rFinal,
-        radiusFun = radiusFun,
-        epochs = epochs,
-    )
+    som_res = trainGigaSOM(som, dInfo; kwargs...)
     undistribute(dInfo)
     return som_res
 end