Hi
When trying to map from a vector of refseq accessions to uniprot id's I get an
error:
refseq <- c("YP_139402", "YP_141320", "YP_820357", "YP_006002448",
"YP_006040776",
"YP_006340074", "YP_006340075", "YP_138830", "YP_819836", "YP_006001858"
)
select(UniProt.ws, keys=refseq, keytype="REFSEQ_PROTEIN", cols=c('UNIPROTKB'))
I have traced the error back to .mergeList. The problem seems to be that when
mapping to UNIPROTKB a data frame with two columns named ACC+ID is produced
making the by='ACC+ID' argument in merge() ambiguous.
The easy fix seems to be to only use 'ACC+ID' for internally keeping track of
the mapping and use another name for the UNIPROTKB id?
best
Thomas Pedersen, PhD student at Technical University of Denmark
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel
