--- a
+++ b/bin/correlation.r
@@ -0,0 +1,28 @@
+#!/usr/bin/env Rscript
+
+library("ggplot2")
+library("ggpubr")
+
+prc_data_norm <- read.csv(file="../data/LungCancerDataset_AllRecords_NORM.csv",head=TRUE,sep=",",stringsAsFactors=FALSE)
+
+ggscatter(prc_data_norm, x = "Stage", y = "Metastasis", 
+          add = "reg.line", conf.int = TRUE, 
+          cor.coef = TRUE, cor.method = "pearson",
+          xlab = "Metastasis", ylab = "Stage")
+
+cat("Pearson correlation between Stage and Metastasis: ")
+cor(as.numeric(prc_data_norm$Stage), as.numeric(prc_data_norm$Metastasis), method = c("pearson"))
+
+cat("Kendall correlation between Stage and Metastasis: ")
+cor(as.numeric(prc_data_norm$Stage), as.numeric(prc_data_norm$Metastasis), method = c("kendall"))
+
+cat("Spearman correlation between Stage and Metastasis: ")
+cor(as.numeric(prc_data_norm$Stage), as.numeric(prc_data_norm$Metastasis), method = c("spearman"))
+
+
+
+
+
+res <- cor.test(prc_data_norm$Stage, prc_data_norm$Metastasis, 
+                    method = "pearson")
+res
\ No newline at end of file