The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: overlapPermTest will not allow randomizing from MASKED hg19
gravatar for peds09
7 months ago by
peds090 wrote:

Using overlapPermTest to compare two GRanges objects - a and b. a is a ChIP seq broad peaks GEO dataset downloaded directly from Annotation Hub in Bioconductor. b is a GRanges list of widths of peaks of my own ChIP seq data. I tried the following code and it worked

random_list = overlapPermTest(a, b, ntimes=1000, genome="BSgenome.Hsapiens.UCSC.hg19", count.once=TRUE)

giving expected output as listed in regioneR tutorial.

But when I try to do the randomization on MASKED hg19 genome - which I downloaded here:



And I run the following:

random_list = overlapPermTest(a, b, ntimes=1000, genome="BSgenome.Hsapiens.UCSC.hg19.masked", count.once=TRUE)

It stalled x 3.  No error message, just stopped computing and crashed my R session.

Any suggestions appreciated,

Thank you in advance,


# Session info:

R version 3.4.4 (2018-03-15)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04 LTS

Matrix products: default
BLAS: /usr/lib/x86_64-linux-gnu/blas/
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/

 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            

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

other attached packages:
 [1] yeastCC_1.18.0                   
 [2] shiny_1.1.0                      
 [3] SummarizedExperiment_1.8.1       
 [4] DelayedArray_0.4.1               
 [5] matrixStats_0.53.1               
 [6] Biobase_2.38.0                   
 [7] BSgenome.Hsapiens.UCSC.hg19_1.4.0
 [8] AnnotationHub_2.10.1             
 [9] regioneR_1.10.0                  
[10] BSgenome_1.46.0                  
[11] Biostrings_2.46.0                
[12] XVector_0.18.0                   
[13] memoise_1.1.0                    
[14] rtracklayer_1.38.3               
[15] GenomicRanges_1.30.3             
[16] GenomeInfoDb_1.14.0              
[17] IRanges_2.12.0                   
[18] S4Vectors_0.16.0                 
[19] BiocGenerics_0.24.0              
[20] BiocInstaller_1.28.0             

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.17                  later_0.7.3                  
 [3] compiler_3.4.4                bitops_1.0-6                 
 [5] tools_3.4.4                   zlibbioc_1.24.0              
 [7] digest_0.6.15                 bit_1.1-14                   
 [9] jsonlite_1.5                  RSQLite_2.1.1                
[11] lattice_0.20-35               rlang_0.2.1                  
[13] pkgconfig_2.0.1               Matrix_1.2-12                
[15] DBI_1.0.0                     curl_3.2                     
[17] yaml_2.1.19                   GenomeInfoDbData_1.0.0       
[19] httr_1.3.1                    bit64_0.9-7                  
[21] grid_3.4.4                    R6_2.2.2                     
[23] AnnotationDbi_1.40.0          XML_3.98-1.11                
[25] BiocParallel_1.12.0           magrittr_1.5                 
[27] blob_1.1.1                    promises_1.0.1               
[29] htmltools_0.3.6               Rsamtools_1.30.0             
[31] GenomicAlignments_ 1.14.2      xtable_1.8-2                 
[33] mime_0.5                      interactiveDisplayBase_1.16.0
[35] httpuv_1.4.4.2                RCurl_1.95-4.10  

regioner • 115 views
ADD COMMENTlink written 7 months ago by peds090
Please log in to add an answer.


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