I know some of the function has imbedded in them some code to obtained the right matrices to perform the analyses.
It would be great to write these bit of code into functions that can be used by other function. An example of this would be a getTraitMatrix function that could replace code on L58-59 of fitDMC.