Create a gist now

Instantly share code, notes, and snippets.

What would you like to do?
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