Script 4

in script 4 :

Remove this line of code metrics_list$unique_taxa <- metrics_list$unique_taxa %>% select(-unique_taxa)

It does not cause issues to me, try to understand why it does for you.

also can you look at a way to produce an excel table with a sheet for each crop (beans, rice etc), and each sheet with 2 columns one with the taxa and the second with how many accessions are recorded in the combined dataset?