# Load the reference data
data(ko_to_kegg_reference)
# View structure
str(ko_to_kegg_reference)
# Get unique pathways
unique_pathways <- unique(ko_to_kegg_reference$pathway_id)
length(unique_pathways)
# Find KOs for a specific pathway
glycolysis_kos <- ko_to_kegg_reference[ko_to_kegg_reference$pathway_id == "ko00010", ]
head(glycolysis_kos)
Run the code above in your browser using DataLab