Hi Antigoni You can try ape (and nlme) with the functions gls and corPagel (the latter to define the covariance structure, ie the evolutionary model) and also caper (I don't know the exact function there).
Cheers Alejandro _____________________________________ Dr. Alejandro Gonzalez Voyer Estación Biológica de Doñana (CSIC) Avenida Américo Vespucio s/n 41092 Sevilla Spain Tel: +34- 954 466700, ext 1749 E-mail: [email protected] Web-site (Under Construction): Group page: http://consevol.org/index.html Personal web-page: http://consevol.org/members/alejandro.html On 28/11/2012, at 18:02, Antigoni Kaliontzopoulou <[email protected]> wrote: > Hello everyone, > > I am trying to implement Liam Revell's suggestion on the evaluation of > Pagel's lambda simultaneous to fitting PGLS to minimize the effects of wrong > model selection (OLS vs. PGLS) on species data (i.e. Revell 2010, Methods in > Ecology and Evolution 1: 319-329). > > Does anyone know if this kind of simultaneous optimization of lambda and PGLS > parameters is presently implemented in R? The obvious place would be > phytools, but I could not find anything similar in that package. > > I would be grateful if you can provide any suggestions. > > Antigoni > > -- > ATTENTION: New e-mail address, please change to [email protected] > > ------------------- > Antigoni Kaliontzopoulou > > Post-doctoral researcher > > CIBIO, Centro de Investigacao em Biodiversidade e Recursos Geneticos > > Campus Agrario de Vairao, 4485-661 Vairao > > PORTUGAL > > Department of Ecology, Evolution, and Organismal Biology > > Iowa State University, Ames, > > Iowa 50011, USA > > _______________________________________________ > R-sig-phylo mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo _______________________________________________ R-sig-phylo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
