[BioC] How to use datapackage homology

[Robert Gentleman]

Hi,
 And so, using version 1.6.2 of all data packages,

 v1=unlist(as.list(rae230aLOCUSID))

 ##drop those with no mappings
  v2=v1[!is.na(v1)]

 ##get the homology codes
 ##at most one per LL so we can unlist
 v3 = unlist(mget(as.character(v2), homologyLL2HGID, ifnotfound=NA))

 #see how many we got
 sum(is.na(v3))
 length(v3)

 #I get about 11000

 now you have them - note that there is a vignette on using the
 homology data in the annotate package; any comments, contributions
  would be welcome

 Robert

-- 
+---------------------------------------------------------------------------+
| Robert Gentleman                 phone : (617) 632-5250                   |
| Associate Professor              fax:   (617)  632-2444                   |
| Department of Biostatistics      office: M1B20                            |
| Harvard School of Public Health  email: rgentlem at jimmy.harvard.edu        |
+---------------------------------------------------------------------------+