Question: bad restore file magic number with AnnotationHub
gravatar for jason.serviss
2.1 years ago by
jason.serviss0 wrote:

I am failing to load ENCODE RNAseq expression data with after downloading via AnnotationHub.

ah = AnnotationHub()
dm <- query(ah, c("EncodeDCC", "Homo sapiens"))


AnnotationHub with 1 record

# snapshotDate(): 2016-11-15 
# names(): AH1376
# $dataprovider: EncodeDCC
# $species: Homo sapiens
# $rdataclass: GRanges
# $title: wgEncodeCshlLongRnaSeqHelas3CellPapGeneGencV7
# $description: 
# $taxonomyid: 9606
# $genome: hg19
# $sourcetype: BED
# $sourceurl:
# $sourcelastmodifieddate: NA
# $sourcesize: NA
# $tags: c("011WC,012WC", "HeLa-S3", "wgEncodeCshlLongRnaSeq", "RnaSeq", "ENCODE Mar 2012 Freeze",
#   "2011-07-13", "2012-04-13", "wgEncodeEH000173", "GSM765402", "Gingeras", "CSHL",
#   "LID16633,LID16634", "iIDR", "cell", "9b6a25c6c8e577edf94e8a28f41f0f3a", "wgEncode", "2x76D",
#   "longPolyA", "Illumina_GA2x", "2e+06", "Nist14", "4583", "gtf", "GeneGencV7") 
# retrieve record with 'object[["AH1376"]]' 


out <- dm[["AH1376"]]


downloading from ‘’
retrieving 1 resource
  |======================================================================================================| 100%
Error: failed to load resource
  name: AH1376
  title: wgEncodeCshlLongRnaSeqHelas3CellPapGeneGencV7
  reason: bad restore file magic number (file may be corrupted) -- no data loaded
In addition: Warning message:
file ‘1376’ has magic number 'chr16'
  Use of save versions prior to 2 is deprecated


R version 3.3.2 (2016-10-31)

Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Sierra 10.12.1

[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] stats4    parallel  stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] GenomicRanges_1.26.1 GenomeInfoDb_1.10.1  IRanges_2.8.1        S4Vectors_0.12.0     AnnotationHub_2.6.1 
[6] BiocGenerics_0.20.0 

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.7                   digest_0.6.10                 mime_0.5                     
 [4] R6_2.2.0                      xtable_1.8-2                  DBI_0.5-1                    
 [7] RSQLite_1.0.0                 BiocInstaller_1.24.0          httr_1.2.1                   
[10] zlibbioc_1.20.0               XVector_0.14.0                curl_2.2                     
[13] Biobase_2.34.0                shiny_0.14.2                  httpuv_1.3.3                 
[16] AnnotationDbi_1.36.0          htmltools_0.3.5               interactiveDisplayBase_1.12.0


Is there a known problem here or have I missed something essential? Thanks in advance for any help you can offer!

ADD COMMENTlink modified 2.0 years ago by Valerie Obenchain ♦♦ 6.7k • written 2.1 years ago by jason.serviss0

Thanks for posting this bug. We are looking into it.


ADD REPLYlink written 2.1 years ago by Valerie Obenchain ♦♦ 6.7k
gravatar for Valerie Obenchain
2.0 years ago by
Valerie Obenchain ♦♦ 6.7k
United States
Valerie Obenchain ♦♦ 6.7k wrote:


I'm finally back with an answer. A similar question came up which I've responded to here:

A: AnnotationHub: two broad peaks from same url are not the same

The file you had problems with was saved with the wrong file extension (.Rda when it was really a gff). It was one of a group of files with several problems either in how the data were processed or how the metadata were specified. I gave a more complete answer on the other post but the gist is that these resources were too dammaged to fix and have been removed from the hub. Similar resources are available via query().

Thanks for reporting this.


ADD COMMENTlink written 2.0 years ago by Valerie Obenchain ♦♦ 6.7k
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: 415 users visited in the last hour