compounds <- c("C00022", "C00153", "C00216", "C01167") organism <- c("M00008") for (compound in compounds) { compound_info <- keggFind("compound", compounds) if(length(compound_info$organism) > 0) { cat(paste(compound, "belongs to", organism, "\n")) } else { cat(paste(compound, "has no associated organism(s) in KEGG\n")) } }
code to explanationWed, 07 Jun 2023