Q matrices now contain rates

All models to rTraitDisc are now passed as matrices, not as characters
(e.g., “ER”).

num_states and rate implicit with model as matrix

incorporate restricted ER model (transition only between adj states)
refactoring, adding restricted er simulation