fitmultiBM function

Fits multi-rate (or multi-trend) discrete character dependent continuous trait evolution