Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts

  • ajsn6c
    started a topic MEDIPS by Lukas Chavez in R

    MEDIPS by Lukas Chavez in R

    I'm running a few quality controls using the most up to date MEDIPS package on my MacBook Air OS X El Capitan with 4GB of memory. I'm able to run the other quality controls: saturation analysis and sequence pattern coverage nicely. However this is the error I receive after running a CpG enrichment check. The file I ran it on is around 3GB and I only selected chr1:

    er = MEDIPS.CpGenrich(file=file,BSgenome=BSgenome, chr.select = chr.select, extend = extend,uniq = uniq)
    Reading bam alignment HCB3_45Low_merged_unique.bam
    Selecting chr1
    Total number of imported short reads: 3214352
    Extending reads...
    Creating GRange Object...
    Keep at most one 1 read mapping to the same genomic location.
    Number of remaining reads: 3211326
    Loading chromosome lengths for BSgenome.Hsapiens.UCSC.hg19...
    Calculating CpG density for given regions...

    *** caught segfault ***
    address 0x7fff4c68f380, cause 'memory not mapped'

    Traceback:
    1: .Call(.NAME, ..., PACKAGE = PACKAGE)
    2: .Call2("new_CHARACTER_from_XStringSet", x, xs_dec_lkup(x), PACKAGE = "Biostrings")
    3: .local(x, ...)
    4: as.character(ans)
    5: as.character(ans)
    6: .local(x, ...)
    7: getSeq(dataset, names = space(x), start = start(x), end = end(x), as.character = TRUE)
    8: getSeq(dataset, names = space(x), start = start(x), end = end(x), as.character = TRUE)

    I haven't found a solid answer online, but is this R's way of telling me I ran out of memory or a bug with Biostrings? I've updated Biostings prior to running this.

    Let me know what you think. Im very new to R programming, so simplifying would be appreciated.

    Thanks!

Latest Articles

Collapse

  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, Yesterday, 09:57 AM
0 responses
16 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, Yesterday, 09:51 AM
0 responses
11 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
66 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
47 views
0 reactions
Last Post SEQadmin2  
Working...