Re: [R-sig-phylo] Ancestral state estimates of continuous traits

2013-03-16 Thread Emmanuel Paradis
: Fri, 15 Mar 2013 14:23:20 To: Alejandro Gonzalezalejandro.gonza...@ebd.csic.es Cc: R-phylo Mailing-listr-sig-phylo@r-project.org Subject: Re: [R-sig-phylo] Ancestral state estimates of continuous traits Hi Alejandro. That must be a bug in ace(...,method=GLS). I would also suggest you check out

Re: [R-sig-phylo] Ancestral state estimates of continuous traits

2013-03-16 Thread Theodore Garland Jr
Cc: R-phylo Mailing-list Subject: Re: [R-sig-phylo] Ancestral state estimates of continuous traits Alejandro and Liam, You can also try method=REML which gives much better estimates of sigma^2 than ML. I think I'll make it the default for continuous characters. Best, Emmanuel -Original

[R-sig-phylo] Ancestral state estimates of continuous traits

2013-03-15 Thread Alejandro Gonzalez
Hello, I am using ape to obtain ancestral state estimates for continuous traits. Two options are available, either maximum likelihood or a GLS method. I am comparing the results of both methods and one difference between the two puzzles me, I hope someone can enlighten me. Under GLS the