Implement SIMMAP in R
SIMMAP creates a stochastic simulation of discrete character trait evolution on a phylogenetic tree. This is a realization from the Bayesian posterior distribution of possible transitions given the tip data and the tree.
A) Stochastic simulation under a given model
B) Infer the stochastic model
a) Should implement an arbitrary size finite state Markov model
b) Should provide fast visualization or realizations
This should be incorporated into the ouch package as a way of generating possible paintings. Will require either a simmap that allows changes only at nodes or a modification of the ouch code to allow changes along a branch.