Skip to content

Instantly share code, notes, and snippets.

@freesci
Created March 25, 2013 20:24
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save freesci/5240374 to your computer and use it in GitHub Desktop.
Save freesci/5240374 to your computer and use it in GitHub Desktop.
Creating an annotation database in Bioconductor for non-model organism
source("http://bioconductor.org/biocLite.R")
library(AnnotationForge)
makeOrgPackageFromNCBI(version="0.0.1",author="PS",maintainer="PS <ps@pawelszczesny.org>",outputDir=".",tax_id=4081,genus="Solanum", species="lycopersicum")
#here you need to wait a bit, few hours maybe
install.packages("org.Slycopersicum.eg.db", repos=NULL, type="source")
library("org.Slycopersicum.eg.db/")
keys<-head(keys(org.Slycopersicum.eg.db))
head(select(org.Slycopersicum.eg.db, keys=keys,cols=c("ENTREZID","ACCNUM","ALIAS","CHR","PMID","REFSEQ","SYMBOL","UNIGENE" ,"GENENAME"))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment