Hi Jorge- I would just fit and compare the models manually, without relying on batch processing of models with fitdAICrc. Each model fitted individually gives you an AIC and a log-likelihood, so it's straightforward to do this. I would not trust comparisons between continuous change models (like SPVAR) and the discrete shift models, because (as in Rabosky 2006 Evolution) there is an overfitting issue that occurs with the multirate Yule shift models.
In general, I would suggest using more recently developed methods/implementations rather than those in laser, such as those in the TreePar or DDD packages, or our BAMM software (www.bamm-project.org) and associated BAMMtools R package. ~Dan On Apr 9, 2014, at 10:40 AM, Jorge Avaria Llautureo wrote: > Hi, > I compared variable rate models in LASER package with the fitdAICrc function, > but I cannot incorporate the yule4rate, SPVAR, EXVAR, and BOTHVAR models. How > I can incorporate these models to be evaluated? > > Thanks, > Jorge, > -- > Msc. Jorge Avaria Llautureo > Lab. Ecología Evolutiva y Filoinformática > Programa de Doctorado en Sistemática & Biodiversidad > Universidad de Concepción Chile. _____________________ Dan Rabosky Assistant Professor & Curator of Herpetology Museum of Zoology & Department of Ecology and Evolutionary Biology University of Michigan Ann Arbor, MI 48109-1079 USA [email protected] http://www-personal.umich.edu/~drabosky http://www.lsa.umich.edu/ummz/ [[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]/
