diff --git a/software/DataAcquisition/drawDiagramSingle.R b/software/DataAcquisition/drawDiagramSingle.R
index 863bb4d..7895bc9 100644
--- a/software/DataAcquisition/drawDiagramSingle.R
+++ b/software/DataAcquisition/drawDiagramSingle.R
@@ -8,9 +8,8 @@ Sys.setenv(LANG = "en")
 drawDiagramSingle <- function(samples) {
 
     # date column should be type of date, not string, so parse it
-    png(filename = paste0("results-", format(Sys.time(), "%d_%b_%Y_%H_%M"), ".png"), width = 800, height = 900, bg = "white")
-
-    # create graphs
+    png(filename = paste0("results-", format(Sys.time(), "%d_%b_%Y_%H_%M"), ".png"), 
+        width = 800, height = 900, bg = "white")
 
     samples$Time = as.POSIXct(samples$Time, format="%Y-%m-%d %H:%M:%S.%OS")
 
@@ -21,8 +20,7 @@ drawDiagramSingle <- function(samples) {
                     grid = TRUE, 
                     col.line = "black",
                     ylab = "counts per minute",
-            ylim=c(min(samples$Counter), max(samples$Counter)),
-
+                    ylim=c(min(samples$Counter), max(samples$Counter)),
                     xlim=c(min(samples$Time), max(samples$Time)),
                     main = "Change over time")
 
diff --git a/software/DataAcquisition/main.R b/software/DataAcquisition/main.R
index a1cebbf..38675ca 100644
--- a/software/DataAcquisition/main.R
+++ b/software/DataAcquisition/main.R
@@ -9,3 +9,7 @@ samples <- read.delim(inputFileName, , sep=",")
 
 drawDiagramSingle(samples)
 
+# TODO move to a separate module
+library(Hmisc)
+rcorr(samples$Counter,samples$DMM, type="pearson")
+