달력

52025  이전 다음

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31

library(biomaRt)

ensembl=useMart("ensembl")  # using ensembl database data

ensembl=useDataset("hsapiens_gene_ensembl",mart=ensembl)   # from ensembl using homosapien gene data


dat = as.matrix(read.table('/Volumes/Data1/data/TCGA-bresat/SAAV/4_20170808_breast105_PAV_ano.txt',sep='\t',header=TRUE))

for(i in 1:nrow(dat))

{

    

    key = strsplit(dat[i,12],'/')

    for(j in length(key))

    {

      id=getBM(attributes=c("hgnc_symbol") ,filters = "ensembl_gene_id",values= key[[1]][j], mart= ensembl)

    }

  

}



Posted by jjbang
|