The AnnotationHub
package provides a client interface to resources stored at the AnnotationHub web service.
library(AnnotationHub)
Start by creating an AnnotationHub
object
ah = AnnotationHub()
ah
object to find a TxDb
object for hg38. How many such records are there?library(AnnotationHub)
ah = AnnotationHub()
query(ah, c('hg38','TxDb'))
txdb = ah[["AH52260"]]
transcripts
method to get a set of GRanges
associated with transcripts.tx = transcripts(txdb)
tx
ah
to find the narrow peak RoadMapEpigenome objects. How many records are available?ah_roadmap = query(ah, c("RoadMap", "narrowpeak"))
length(ah_roadmap)
dnase = query(ah_roadmap, 'dnase')
dnase
dnase_track = dnase[['AH45295']]
dnase_track
class(dnase_track)