Question: RankProd crashing - stability issues?
0
gravatar for bhgyu
22 months ago by
bhgyu30
bhgyu30 wrote:

Hi guys,

Anyone know why RankProd crashes R Studio when there are more than say 100 sample cases? 

It consistently crashes R Studio for experiments with a larger number of subjects. 

No error code or anything. R just encounters a fatal error. 

Cheers

> sessionInfo()
R version 3.3.2 (2016-10-31)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X El Capitan 10.11.5

locale:
[1] en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] RankProd_3.0.0 Rmpfr_0.6-1    gmp_0.5-13.1  

loaded via a namespace (and not attached):
[1] tools_3.3.2
> 

 

Crash report: 

Process:               RStudio [64386]
Path:                  /Applications/RStudio.app/Contents/MacOS/RStudio
Identifier:            org.rstudio.RStudio
Version:               1.0.136 (1.0.136)
Code Type:             X86-64 (Native)
Parent Process:        ??? [1]
Responsible:           RStudio [64386]
User ID:               501

Date/Time:             2017-06-16 16:13:17.807 +0100
OS Version:            Mac OS X 10.11.5 (15F34)
Report Version:        11
Anonymous UUID:        C80C5A19-D1DA-B4D5-393C-1C60B848093A

Sleep/Wake UUID:       56D8F5A1-BDD4-48C9-BD19-AF2EF89B7265

Time Awake Since Boot: 500000 seconds
Time Since Wake:       3500 seconds

System Integrity Protection: enabled

Crashed Thread:        0  Dispatch queue: com.apple.main-thread

Exception Type:        EXC_BAD_ACCESS (SIGSEGV)
Exception Codes:       KERN_INVALID_ADDRESS at 0x0000000011288c38
Exception Note:        EXC_CORPSE_NOTIFY

VM Regions Near 0x11288c38:
--> 
    __TEXT                 0000000100000000-00000001002d7000 [ 2908K] r-x/rwx SM=COW  /Applications/RStudio.app/Contents/MacOS/RStudio

Application Specific Information:
objc_msgSend() selector name: _hasActiveAppearanceIgnoringKeyFocus

rankprod crash • 526 views
ADD COMMENTlink modified 22 months ago • written 22 months ago by bhgyu30

These types of question usually require a simple reproducible example, and information about the specific software you use. Please update your question to (a) include a reproducible example, e.g., using a data set from RankProd or simple simulated data and (b) include the output of the sessionInfo() command. For the later, I have

> sessionInfo()
R version 3.4.0 Patched (2017-05-24 r72729)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.2 LTS

Matrix products: default
BLAS: /home/mtmorgan/bin/R-3-4-branch/lib/libRblas.so
LAPACK: /home/mtmorgan/bin/R-3-4-branch/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] RankProd_3.2.0 Rmpfr_0.6-1    gmp_0.5-13.1  

loaded via a namespace (and not attached):
[1] compiler_3.4.0

Edit your question with this example and information, and ADD COMMENT to indicate that you have done so.

ADD REPLYlink written 22 months ago by Martin Morgan ♦♦ 23k

Ok thank you.

ADD REPLYlink modified 22 months ago • written 22 months ago by bhgyu30

You need to show us how you produce the problem -- a few lines of code that I can copy and paste into my own session and see what happens. You need to be using the current version of Bioconductor (R version 3.4.0, BiocInstaller version 1.26.0). It looks like you have an old version of R and RankProd. Install the current R (version 3.4.0) and be sure to follow standard Bioconductor installation instructions

source("https://bioconductor.org/biocLite.R")
biocLite("RankProd")

Since you already have a confused installation, be sure to try and validate your updated installation with

BiocInstaller::biocValid()

 

ADD REPLYlink written 22 months ago by Martin Morgan ♦♦ 23k

Hi Martin, 

Thank you. I did the above yesterday. I updated R to 3.4 and RankProd to 3.2

Bear in mind though, this all does work with a data set that has only 20 samples. Add 180 and it keeps crashing. 

I also should add that I did gave this a go on a more powerful computer with the same problem (2.4Ghz, Core i7).

I run the following:

norm <- read.delim("norm.txt", header=T)

GSE.cl <- c(rep(1, each=90), rep(0, each=90))

origin <- rep(1,6)

symbol <- mapIds(huex10sttranscriptcluster.db,keys = as.vector(rownames(norm)),keytype="PROBEID",column = "ENTREZ")
gnames.27 <- symbol
RP <- RP.advance(norm, GSE.cl, origin, gene.names=gnames.27)

 

PS: Validation is True.

ADD REPLYlink modified 22 months ago • written 22 months ago by bhgyu30

When I try and reproduce your code, I run into the following problem, for obvious reasons

> norm <- read.delim("norm.txt", header=T)
Error in file(file, "rt") : cannot open the connection
In addition: Warning message:
In file(file, "rt") :
  cannot open file 'norm.txt': No such file or directory

Maybe you can use rnorm() or similar to make your example reproducible by others, or use an existing data set (e.g., library(ALL)) that shows the problem.

 

ADD REPLYlink written 22 months ago by Martin Morgan ♦♦ 23k

Sorry,

Let's try this: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE85426

The data set has to be large to test this so I cannot think of anything but an external one. Usually examples in R are rather small. 

ADD REPLYlink modified 22 months ago • written 22 months ago by bhgyu30

Hi,

I'm having the same issue. Did you find a solution?

Thanks

ADD REPLYlink written 16 months ago by tcalvo30
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 16.09
Traffic: 90 users visited in the last hour