Stochastic modeling of phylogenies raises five questions that have received varying levels of attention from quantitatively inclined biologists. 1) How large do we expect (from the model) the ratio of maximum historical diversity to current diversity to be? 2) From a correct phylogeny of the extant species of a clade, what can we deduce about past speciation and extinction rates? 3) What proportion of extant species are in fact descendants of still-extant ancestral species, and how does this compare with predictions of models? 4) When one moves from trees on species to trees on sets of species (whether traditional higher order taxa or clades within PhyloCode), does one expect trees to become more unbalanced as a purely logical consequence of tree structure, without signifying any real biological phenomenon? 5) How do we expect that fluctuation rates for counts of higher order taxa should compare with fluctuation rates for number of species? We present a mathematician's view based on an oversimplified modeling framework in which all these questions can be studied coherently.

Simple stochastic models for phylogenetic trees on species have been well studied. But much paleontology data concerns time series or trees on higher-order taxa, and any broad picture of relationships between extant groups requires use of higher-order taxa. A coherent model for trees on (say) genera should involve both a species-level model and a model for the classification scheme by which species are assigned to genera. We present a general framework for such models, and describe three alternate classification schemes. Combining with the species-level model of Aldous and Popovic (Adv Appl Probab 37:1094-1115, 2005), one gets models for higher-order trees, and we initiate analytic study of such models. In particular we derive formulas for the lifetime of genera, for the distribution of number of species per genus, and for the offspring structure of the tree on genera.