genbankr throwing C stack usage error when reading in files >500Kb
0
0
Entering edit mode
Lauren • 0
@9541a822
Last seen 3.0 years ago
United States

In my hands, the readGenBank function in the genbankr package works very well with smaller GenBank files (<500 KB), but gives me a C stack usage error with any larger files.

In the example below chrI is 442Kb and chrIII is 638Kb. Am I doing something wrong and/or do the larger files need some pre-treatment?

I've tried other files as well - larger files always throw this error. Small files work just fine. I am currently running genbankr version 1.20.0 on R version 4.1.2.

Thanks!


    > library(genbankr)

    > chrI <- readGenBank("/Users/USER/Desktop/chrI.gb")
    Translation product seems to be missing for 1 of 101 CDS annotations. Setting to ''
    No exons read from genbank file. Assuming sections of CDS are full exons
    No transcript features (mRNA) found, using spans of CDSs

    > chrII <- readGenBank("/Users/USER/Desktop/chrII.gb")
    Error: C stack usage  7970208 is too close to the limit


    > sessionInfo()
    R version 4.1.2 (2021-11-01)
    Platform: x86_64-apple-darwin17.0 (64-bit)
    Running under: macOS Big Sur 11.5.2

    Matrix products: default
    LAPACK: /Library/Frameworks/R.framework/Versions/4.1/Resources/lib/libRlapack.dylib

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

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

    other attached packages:
     [1] BSgenome_1.60.0             rtracklayer_1.52.1          GenomicAlignments_1.28.0   
     [4] Rsamtools_2.8.0             Biostrings_2.60.2           XVector_0.32.0             
     [7] SummarizedExperiment_1.22.0 MatrixGenerics_1.4.3        matrixStats_0.61.0         
    [10] GenomeInfoDbData_1.2.6      GenomicFeatures_1.44.2      AnnotationDbi_1.54.1       
    [13] Biobase_2.52.0              GenomicRanges_1.44.0        GenomeInfoDb_1.28.4        
    [16] IRanges_2.26.0              S4Vectors_0.30.2            BiocGenerics_0.38.0        
    [19] genbankr_1.20.0
genbankr • 551 views
ADD COMMENT

Login before adding your answer.

Traffic: 739 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6