fit_and_compare_bm_models function

Fit and compare Brownian Motion models for multivariate trait evolution between two data sets.