Skip to content
Prev 7316 / 21312 Next

[Bioc-devel] AnnotationHubData Error: Access denied: 530

Hi Johannes, Sonali,
On 04/10/2015 09:40 AM, Arora, Sonali wrote:
A GRanges object is not the same as a GTF file and I guess Johannes
wants access to the GTF file. Are these GTF files available on
AnnotationHub?

@Johannes - Here is one alternative: You could take a different approach
and implement some equivalent of makeTxDbFromGRanges() for EnsDb
objects. So people could just do:

   library(ensembldb)
   ensdb <- makeEnsDbFromGRanges(gtf[[1]])

like they can do right now with makeTxDbFromGRanges():

   library(GenomicFeatures)
   txdb <- makeTxDbFromGRanges(gtf[[1]])

That way you don't need a recipe or try to add things to AnnotationHub 
at all.

@Sonali - These GRanges objects I get from AnnotationHub have no genome
information and their seqlevels are not sorted:

   > seqinfo(gtf[[1]])
   Seqinfo object with 22 sequences from an unspecified genome; no 
seqlengths:
     seqnames seqlengths isCircular genome
     X              <NA>       <NA>   <NA>
     9              <NA>       <NA>   <NA>
     8              <NA>       <NA>   <NA>
     7              <NA>       <NA>   <NA>
     6              <NA>       <NA>   <NA>
     ...             ...        ...    ...
     12             <NA>       <NA>   <NA>
     11             <NA>       <NA>   <NA>
     10             <NA>       <NA>   <NA>
     1              <NA>       <NA>   <NA>
     MT             <NA>       <NA>   <NA>

I know it's easy enough to sort the seqlevels with sortSeqlevels() but
what about having these things done by the recipe instead?

Thanks,
H.