列印

Please login or register.
一小時 一天 一週 一個月 永遠
請輸入帳號, 密碼以及預計登入時間

[admin]

http://www.ensembl.org/info/data/biomart/biomart_r_package.html

程式碼: [Select]

source("http://bioconductor.org/biocLite.R")
biocLite("biomaRt")

程式碼: [Select]

library(biomaRt)
listEnsembl()
listEnsembl("GRCh=37")
listEnsembl(version=78)

[admin]

http://stackoverflow.com/questions/41490657/best-way-to-get-list-of-snps-by-gene-id/41496733

程式碼: [Select]

library(biomaRt)

#load the human variation data
variation = useEnsembl(biomart="snp", dataset="hsapiens_snp")

#look up a single gene and get SNP data
getBM(attributes = c(
  "ensembl_gene_stable_id",
  'refsnp_id',
  'chr_name',
  'chrom_start',
  'chrom_end',
  'minor_allele',
  'minor_allele_freq'),
  filters = 'ensembl_gene',
  values ="ENSG00000166813",
  mart = variation
)