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 = 3.6, df = 280, p-value = 3e-04
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.09902 0.32404
## sample estimates:
##    cor 
## 0.2144

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 = 4.5, df = 280, p-value = 1e-05
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1481 0.3679
## sample estimates:
##    cor 
## 0.2614

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 = 4.7, df = 280, p-value = 5e-06
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.1577 0.3763
## sample estimates:
##    cor 
## 0.2705

tt <- plot_linear_scatter(df=merged[, c(3, 6)])

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

tt$scatter