I knew there would be a mistake somewhere. Here is the correction:
## C <- 1+(lambda-1)*(1-mat) ## Lambda correction WRONG!!
C <- (lambda*mat)+(1-lambda)*diag(diag(mat)) ## correct, I think..
I tested it on a toy example of a tree with 5 species and it works fine..
I hope there are no other
Dear All,
Okay, so I tried the code that Dr. Blomberg provided, and it doesn't seem
to work. I rewrote the code slightly in order for it to try and predict the
body mass of the taxa used to calculate the line, rather than some new
taxon, because I need to do that in order to calculate prediction
Hi Russell,
And, for a fully PGLS and Bayesian model, see this:
https://pubmed.ncbi.nlm.nih.gov/17344851/
Best, Chris
On Thu, Jul 1, 2021 at 1:10 AM Theodore Garland
wrote:
>
> Russell,
> Please read this paper:
> https://pubmed.ncbi.nlm.nih.gov/10718731/
> Cheers
> Ted
>
>
> On Wed, Jun 30,
