Skip to content

Instantly share code, notes, and snippets.

@dtenenba

dtenenba/ahprob.R

Created September 10, 2015 18:09
Show Gist options
  • Save dtenenba/ba93066fe8385590b6d9 to your computer and use it in GitHub Desktop.
Save dtenenba/ba93066fe8385590b6d9 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
ahprob.R
source("https://bioconductor.org/biocLite.R")
biocLite(c("AnnotationHub", "rtracklayer"))
library(AnnotationHub)
ah <- AnnotationHub()
ah <- subset(ah, species == "Homo sapiens")
qhs <- query(ah, "RefSeq")
genes <- qhs[qhs$genome == "hg19" & qhs$title == "RefSeq Genes"]
genes <- qhs[[1]]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment