bad restore file magic number with AnnotationHub
1
0
Entering edit mode
@jasonserviss-9818
Last seen 2.4 years ago

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

library(AnnotationHub)
ah = AnnotationHub()
dm <- query(ah, c("EncodeDCC", "Homo sapiens"))
dm["AH1376"]

 

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: http://hgdownload.cse.ucsc.edu/goldenpath/hg19/encodeDCC/wgEncodeCshlLongRnaSeq/wgEncodeCshlLo...
# $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"]]

require(“GenomicRanges”)

downloading from ‘https://annotationhub.bioconductor.org/fetch/1376’
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

sessionInfo()

R version 3.3.2 (2016-10-31)

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

locale:
[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!

annotationhub • 789 views
ADD COMMENT
0
Entering edit mode

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

Valerie

ADD REPLY
0
Entering edit mode
@valerie-obenchain-4275
Last seen 2.1 years ago
United States

Hi,

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.

Valerie

ADD COMMENT

Login before adding your answer.

Traffic: 368 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