We need the input .biom file in datasets/ch5/rich_high_count_otu.biom and the metacoder and RColorBrewer packages.