if the dataset follow the normal distribution,
run Pearson's correlation coefficient
otherwise,
run Spearman's rank correlation coefficient
R script;
>cor(dataset, method = "pearson")
>cor(var1, var2, method = "spearman")
>cor(var1, var2, method = "spearman")