For this recipe, we'll use the treespace package to compute distances and clusters. We'll use ape and adegraphics for accessory loading and visualization functions. The input data here will be all 20 files in datasets/ch4/gene_trees, each of which is a Newick-format tree for a single gene in each of 15 species.