我有物种的分类 ID,我可以从 NCBI ( https://www.ncbi.nlm.nih.gov/Taxonomy/TaxIdentifier/tax_identifier.cgi ) 获得物种和属名。但我想要这些数据中的门、类、顺序。
我尝试过taxize package,但不适用于大型数据集。
specieslist <- c("Clostridium", "Clostridium", "Achromobacter", "Achromobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Acinetobacter", "Actinomyces", "Actinomyces", "Aeromonas", "Agrococcus", "Alcanivorax", "Alkalihalobacillus", "Alloprevotella", "Aminobacterium", "Amniculibacterium", "Anaerocolumna", "Anaerocolumna", "Anaerocolumna", "Asticcacaulis", "Atopobium", "Bacillus", "Bacillus", "Bacteroidales", "Bacteroides", "Bacteroides", "Bacteroides", "Bacteroides", "Bacteroides", "Bacteroides", "Barnesiella", "Bifidobacterium", "Blochmannia", "Bordetella", "Brevibacillus", "Buchnera", "Burkholderia", "Butyricimonas", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Campylobacter", "Capnocytophaga", "Capnocytophaga", "Capnocytophaga", "Chroococcidiopsis", "Citrobacter", "Clostridium", "Clostridium", "Clostridium", "Clostridium", "Corynebacterium", "Corynebacterium", "Corynebacterium", "Corynebacterium", "Cutibacterium", "Dialister", "Dolosigranulum", "Enterobacter", "Enterococcus", "Entomoplasma", "Escherichia", "Escherichia", "Escherichia", "Eubacterium", "Fermentimonas", "Frankia", "Fusobacterium", "Fusobacterium", "Fusobacterium", "Fusobacterium", "Fusobacterium", "Fusobacterium", "Gemella", "Haemophilus", "Haemophilus", "Halomonas", "Hydrogenophaga", "Ilyobacter", "Klebsiella", "Klebsiella", "Klebsiella", "Klebsiella", "Klebsiella", "Kocuria", "Kytococcus", "Lachnoanaerobaculum", "Lachnospira", "Lachnospiraceae", "Lachnospiraceae", "Lacrimispora", "Lactobacillus", "Lactobacillus", "Lactobacillus", "Lactobacillus", "Lactobacillus", "Lancefieldella", "Lautropia", "Leptotrichia", "Leptotrichia", "Leptotrichia", "Leptotrichia", "Leptotrichia", "Leptotrichia", "Leptotrichia", "Leptotrichia", "Ligilactobacillus", "Limosilactobacillus", "Luteimonas", "Lysinibacillus", "Lysobacter", "Lysobacter", "Lysobacter", "Magnetospirillum", "Marivirga", "Megasphaera", "Megasphaera", "Meiothermus", "Methylobacterium", "Methylobacterium", "Methylobacterium", "Methylobacterium", "Methylobacterium", "Microbacterium", "Microbacterium", "Microbacterium", "Microbacterium", "Micrococcus", "Muribaculaceae", "Muribaculum", "Muribaculum", "Neisseria", "Neisseria", "Neisseria", "Neisseria", "Neisseria", "Neisseria", "Neisseria", "Nocardioides", "Nocardioides", "Paludibacter", "Pantoea", "Paracoccus", "Paracoccus", "Paraprevotella", "Pasteurella", "Petrimonas", "Phenylobacterium", "Phocaeicola", "Phocaeicola", "Phocaeicola", "Phyllobacterium", "Polaribacter", "Pontibacter", "Pontibacter", "Porphyromonas", "Porphyromonas", "Porphyromonas", "Porphyromonas", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Prevotella", "Proteus", "Pseudoleptotrichia", "Pseudomonas", "Pseudonocardia", "Pseudonocardia", "Raoultella", "Rheinheimera", "Romboutsia", "Roseivirga", "Roseococcus", "Rothia", "Rothia", "Rubrobacter", "Rubrobacter", "Rufibacter", "Saccharomonospora", "Saccharopolyspora", "Saccharopolyspora", "Salinivirga", "Salmonella", "Schaalia", "Sedimentisphaera", "Selenomonas", "Selenomonas", "Selenomonas", "Selenomonas", "Selenomonas", "Selenomonas", "Shigella", "Skermanella", "Sphingosinicella", "Spirosoma", "Staphylococcus", "Staphylococcus", "Stenotrophomonas", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptococcus", "Streptomyces", "Tannerella", "Tannerella", "Thermovirga", "Treponema", "Treponema", "Treponema", "Treponema", "Treponema", "Treponema", "Veillonella", "Veillonella", "Veillonella", "Veillonella", "Veillonella")
t <- tax_name(query = c(specieslist), get = c("phylum","class", "order", "family", "genus"), db = "ncbi")
请问有什么建议吗?