mersmers comparison

Load the two analyses

savefile <- paste0("combined_", mers, ".rda")
old_data <- new.env()
load(file=file.path("small_previous", savefile), envir=old_data)
old_data <- old_data[["ret"]]
new_data <- new.env()
load(file=file.path("small_new", savefile), envir=new_data)
new_data <- new_data[["ret"]]

merged <- merge(x=old_data[["data"]][[1]][, c("limma_logfc", "deseq_logfc", "edger_logfc")],
                y=new_data[["data"]][[1]][, c("limma_logfc", "deseq_logfc", "edger_logfc")],
                by="row.names")
cor.test(x=merged[, "limma_logfc.x"], y=merged[, "limma_logfc.y"])

## 
##  Pearson's product-moment correlation
## 
## data:  merged[, "limma_logfc.x"] and merged[, "limma_logfc.y"]
## t = 11, df = 2700, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1642 0.2369
## sample estimates:
##    cor 
## 0.2008

cor.test(x=merged[, "deseq_logfc.x"], y=merged[, "deseq_logfc.y"])

## 
##  Pearson's product-moment correlation
## 
## data:  merged[, "deseq_logfc.x"] and merged[, "deseq_logfc.y"]
## t = 11, df = 2700, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1731 0.2455
## sample estimates:
##    cor 
## 0.2095

cor.test(x=merged[, "edger_logfc.x"], y=merged[, "edger_logfc.y"])

## 
##  Pearson's product-moment correlation
## 
## data:  merged[, "edger_logfc.x"] and merged[, "edger_logfc.y"]
## t = 11, df = 2700, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1777 0.2499
## sample estimates:
##    cor 
## 0.2141

tt <- plot_linear_scatter(df=merged[, c(2,5)])

## Used Bon Ferroni corrected t test(s) between columns.

tt$scatter