Search
Question: ReactomePA / Reactome: Crash when plotting
0
gravatar for martin.busch
7 days 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

pls make sure you are using latest release version.

 

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

 

ADD REPLYlink written 6 days ago by Guangchuang Yu990

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) {
      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")
}

# 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]

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[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]

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")

 

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 6 days 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.

 

ADD REPLYlink written 6 days ago by Guangchuang Yu990

ADD REPLYlink written 6 days ago by Guangchuang Yu990

Thanks for your help, I will try to implement your script.

 

Best,

Martin

ADD REPLYlink written 6 days ago by martin.busch0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 398 users visited in the last hour