Tried new code from Mahler's Continuos traits tutorial. Added a step to sort trait data based on the order of the tip labels before calculating the vcv matrices; I need to fix this in my other analyses