Question: biostrings report "XStringSet object is too big to be unlisted" error
gravatar for zwfanwz678
5 months ago by
zwfanwz6780 wrote:

hello everone,

I using Biostrings::readDNAStringSet to cut some sequence from hg19.fa, but an error occurs when I run the program(as follows).

seq <- Biostrings::readDNAStringSet("hg19.fa")
genome <- GenomeInfoDb::Seqinfo(genome = NA_character_)
gr_a <- BiocGenerics::unstrand(rtracklayer::import("file.bed", genome = genome))
a <- seq[gr_a]

gr_a is a GRanges object, when run a <- seq[gr_a], R report

"Error in .Call2("XStringSet_unlist", x, PACKAGE = "Biostrings") : 
  XStringSet object is too big to be unlisted (would result in an XString
  object of length 2^31 or more)"

I try these codes several times, but sometimes the code runs very well, sometimes it report the error above.

what's wrong with these code?


ADD COMMENTlink written 5 months ago by zwfanwz6780
Please log in to add an answer.


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