Advertisement
favorov

RefSeq to HGNC names

Aug 26th, 2014
266
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 0.42 KB | None | 0 0
  1. if (!require('Homo.sapiens'))
  2. {
  3.   source("http://bioconductor.org/biocLite.R")
  4.   biocLite("Homo.sapiens")
  5.   library('Homo.sapiens')  
  6. }
  7.  
  8. isonames<-c('NM_001276435',
  9. 'NM_001276436',
  10. 'NM_001276437',
  11. 'NM_001276438',
  12. 'NM_001276439',
  13. 'NM_002243',
  14. 'NM_170736',
  15. 'NM_170737',
  16. 'NM_003966')
  17.  
  18.  
  19. geneSymbols <- select(
  20.     org.Hs.eg.db,
  21.     keys=isonames,
  22.     cols=c('SYMBOL'),
  23.     keytype='REFSEQ'
  24. )
  25.  
  26. print(geneSymbols$SYMBOL)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement