Dear all,
I'm using the mvMORPH package to fit the BM, EB, OU1, and OUM
macroevolutionary models to the PCA axes of the log-shape ratios of
whip-spiders. The simmap tree contains 69 species and 4 habitat states
(cave = 24, cave/forest = 6, city = 5, forest = 34).
At first, I used mvgls but couldn't find a way of extracting model
parameters (e.g., sigma, beta, r, and alfa), so I moved to the separate
functions for each model (mvBM, mvEB, and mvOU).
The two approaches produced different results, with the best mvgls model
(according to GIC) being OU1, and the best mvOU model (according to AIC)
being OUM.
Could someone please help me with the following questions?
1) How to extract model parameters from models fitted using mvgls?
2) I understand that mvgls incorporates penalized likelihood, whereas mvOU
does not. Does that explain the discrepancy in the results? Which approach
would be more appropriate in my case?
3) When I fit the mvOU (OU1) model, it says it has converged, but it also
says "Unreliable solution (Likelihood at a saddle point)". Does anyone know
how serious this is and how to deal with it?
Thank you!
Paulo
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list - [email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/[email protected]/
