[BioC] KEGGSOAP questions: how to handle multiple annotations in R data frames, and why does a
Nianhua Li
nli at fhcrc.org
Tue Nov 28 02:37:56 CET 2006
Hi, Alan,
I just had a quick look of your email. Not sure whether the following script is
helpful or not:
## your code
cpdID<-c(1,2,3,4,5,6)
mass<-c(129.0426, 147.0532, 208.0848, 220.0848, 204.0899, 777.0317)
RT<-c(1,2,3,4,5,6)
ppmerror<-c(4,11,75,7,21,55)
floatmass=NULL
for (i in 1:length(cpdID)) {
floatmass[i]<-if(ppmerror[i]<10) {1e-5*mass[i]} else{(ppmerror[i]/10^6)*mass[i]}
}
## something else
res <- mapply(search.compounds.by.mass, mass, floatmass)
final_res <- cbind(cpdID, mass, floatmass, res)
sessionInfo()
...
KEGGSOAP SSOAP RCurl XML
"1.9.1" "0.4-1" "0.8-0" "1.2-0"
More information about the Bioconductor
mailing list