2011年5月26日星期四

get genetic location for gene names - using biomaRt library

From the biomaRt vignette:

library(biomaRt)

ensembl=useMart("ensembl", dataset="hsapiens_gene_ensembl")

getBM(c("ensembl_gene_id","hgnc_symbol","start_position","end_position"),
            filters = c("chromosome_name", "start","end"), 
            values = list(16, 1100000, 1250000), 
            mart = ensembl)

没有评论:

发表评论