Question: ReactomePA / Reactome: Crash when plotting
0
16 months ago by
martin.busch0 wrote:

Hi everybody,

ReactomePA used to work fine in my handy, however, when trying to plot results it recently started to crash Rstudio completely. I cannot see an appearent reason, e.g. I can dotplot a result, but using emapplot of the same result one line of code later it completely crashes. Ocassionally I can perform two plots (dotplot and emapplot) of one dataset, but when using the next one it starts to crash again. Is there any known issue? as.data.frame(result) always works nicely and displays resulting terms... it seems a bit random. Any ideas on how to track down the bug? Any help appreciated - thank you in advance!

Best,

Martin

written 16 months ago by martin.busch0

pls make sure you are using latest release version.

If the issue still exists when using latest release version, pls provide reproducible example.

Thank your for getting in touch. The version I am using is 1.24.0. An example would be the following:

rm(list=ls())
library(ReactomePA)

dataSets <- c("http://martin-busch.com/work/set1.csv",
"http://martin-busch.com/work/set2.csv",
"http://martin-busch.com/work/set3.csv",
"http://martin-busch.com/work/set4.csv"
)

for (dataSet in dataSets) {
result <- enrichPathway(data$x, organism = "human", pvalueCutoff = 0.05, pAdjustMethod = "BH", qvalueCutoff = 0.2) dotplot(result) emapplot(result, color = "pvalue") } # Manual dataSet <- dataSets[1] data <- read.table(dataSet, header=T) result <- enrichPathway(data$x, organism = "human", pvalueCutoff = 0.05, pAdjustMethod = "BH", qvalueCutoff = 0.2)

dotplot(result)
emapplot(result, color = "pvalue")

dataSet <- dataSets[2]

result <- enrichPathway(data$x, organism = "human", pvalueCutoff = 0.05, pAdjustMethod = "BH", qvalueCutoff = 0.2) dotplot(result) emapplot(result, color = "pvalue") dataSet <- dataSets[3] data <- read.table(dataSet, header=T) result <- enrichPathway(data$x, organism = "human", pvalueCutoff = 0.05, pAdjustMethod = "BH", qvalueCutoff = 0.2)

dotplot(result)
emapplot(result, color = "pvalue")

dataSet <- dataSets[4]

result <- enrichPathway(data$x, organism = "human", pvalueCutoff = 0.05, pAdjustMethod = "BH", qvalueCutoff = 0.2) dotplot(result) emapplot(result, color = "pvalue") The loop runs without any problem, but when I go "manually" trough the datasets, the first two plots are no problem, but the next one causes an immediate crash of Rstudio (1.1.447, Win10/64bit). ADD REPLYlink written 16 months ago by martin.busch0  > res <- lapply(dataSets, function(dataSet) { + data <- read.table(dataSet, header=T) + enrichPathway(data$x, organism = "human", pvalueCutoff = 0.05, pAdjustMethod = "BH", qvalueCutoff = 0.2)
+ })

> emp = lapply(res, emapplot)
> cowplot::plot_grid(plotlist=emp, ncol=2)


can't reproduce any issue.

Best,

Martin