Re: [R-sig-phylo] MRM(), AICc, and other ways to compare effects

2012-09-18 Thread Roland Sookias
Cheers Luke and Alejandro Yes, I've been trying to use Rob's method. However I can't see how it allows incorporation of the other variables (communism, religion) I want as they are categorical/factors. How would you go about including these? I am interested in trying the lnam function. However I

[R-sig-phylo] R: variation in rates over time

2012-09-18 Thread Paolo Piras
Hi look also for motmot and auteur packages best Paolo Da: r-sig-phylo-boun...@r-project.org [r-sig-phylo-boun...@r-project.org] per conto di Matt Pennell [mwpenn...@gmail.com] Inviato: martedì 18 settembre 2012 1.46 A: Jason S Cc:

Re: [R-sig-phylo] MRM(), AICc, and other ways to compare effects

2012-09-18 Thread Alejandro Gonzalez V
Hi Roland, There are two functions for Rob's method to combine phylogenetic and spatial signal, although both - as far as I can tell - are based on gls methods, one uses independent contrasts, and thus I'm not sure if it is set up in such a way to allow you to include factors in the model.

Re: [R-sig-phylo] variation in rates over time, unexpected message when using Brownie.lite

2012-09-18 Thread Agus Camacho
Dear all, I tried to test something similar to Jason, but relating the change in rates to the acquisition of a trait, instead to a specific date. For that, I used a phylogeny with ten taxa, without outgroup. I used exactly the same script gently proposed by Liam but got the following error

Re: [R-sig-phylo] variation in rates over time, unexpected message when using Brownie.lite

2012-09-18 Thread Liam J. Revell
Hi Agus. It looks like in your example the Hessian matrix of partial second derivatives (the curvature) at the optimum is numerically singular. I'm not sure what this means - it could be that the surface is very flat, or alternatively that it is highly correlated. Is the error fatal, or do