fitmultiBM function

Fits multi-rate discrete character dependent continuous trait evolution