Entering edit mode
5.8 years ago
naveen_gabriel
•
0
To fit the brownian motion for independent evolution of traits I used below code:
rownames(carni)<-gsub("_",".",rownames(carni)) #Renaming to match the tip of tree
#Independent Evolution
brown_indepen <- mvBM(carni_tree,carni,param=list(constraint="diagonal"))
mvBM is the function from mvMORPH package. It seems by giving constraint="diagonal" would fit the model for traits to be evolved independently. I would like to know the idea behind it. Also, I got the below output after running the code. I would also like to know what below output means.
LogLikelihood: -589.2373
AIC: 1186.475
AICc: 1186.771
4 parameters
Estimated rate matrix
______________________
size range
size 73.70335 0.00000
range 0.00000 13.54524
Estimated root state
______________________
size range
theta: 38.43947 13.78439
Size and range were the two traits.
