This function simulates the evolution of a group of traits evolving as independent Brownian motions with inverse gamma distributed rates.

1 | ```
norminvgamma.shift.sim.group(phy, alpha, beta, rates, n)
``` |

`phy` |
an ape format phylogeny on which to simulate |

`alpha` |
the shape parameter of the inverse gamma distribution |

`beta` |
the scale parameter of the inverse gamma distribution |

`rates` |
a vector of rates, with each entry corresponding to an edge of phy. Rates should be in the same order as edges in phy$edge |

`n` |
the number of traits to simulate |

A matrix with each row corresponding to a species and each column corresponding to a trait

1 2 | ```
data(yeast)
norminvgamma.shift.sim.group(yeast.tree,2,2,rep(1,6),10)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.