Re: [R-sig-phylo] Cophenetic repeats

2012-05-30 Thread Rob Lanfear
Hi Ben, Just to follow up - if your tree is ultrametric then any pair of tips which span the root will have the same distance (which will be equal to twice the total depth of the tree). Other nodes in the tree will also create equal distances in the same way, but for a reasonably balanced tree

[R-sig-phylo] mccr test for non ultrametric trees

2012-05-30 Thread pasquale.r...@libero.it
Hi all, I would enquire the list about a simple issue. Is there any method implemented to test for changes in diversification rate as applied to fossil (non- ultrametric) trees?. As far as I understand, the methods so far available work on ultrametric trees (I've inspected laser's

[R-sig-phylo] Non-parametric alternative to phylogenetic ANOVA?

2012-05-30 Thread Karin Schneeberger
Dear all I'm trying to compare one trait across three (unordered categorical) groups including 25 species (let's say for example basal metabolic rate of aquatic, terrestrial and aerial mammals). If the data would be normally distributed, I would simply use a phylogenetic ANOVA including a

Re: [R-sig-phylo] Non-parametric alternative to phylogenetic ANOVA?

2012-05-30 Thread Liam J. Revell
Hi Karin. GLS with x as a factor is a generalized ANOVA which assumes [in the case of gls(...,correlation=corBrownian)] that the residual error in the ANOVA model has evolved by Brownian evolution. If you read your data into data frame Z with row names as species names, for instance:

Re: [R-sig-phylo] mccr test for non ultrametric trees

2012-05-30 Thread David Bapst
Pasquale- This isn't a feasible solution, because the branch lengths of a paleo-tree are a function of birth, death and sampling rates. I'll be discussing work that Matt Pennell, Emily King and I have been doing relating to this issue next month at Evolution. -Dave On Wed, May 30, 2012 at 8:23

[R-sig-phylo] brunch

2012-05-30 Thread Yanthe Pearson
Hello, I have been trying to figure out a few details in the caper functions brunch and crunch. Say I want to run PIC for Y as a function of X (as seen in Garland and Ives 2000), my understanding is I would have to find the IC's for Y and X independently and apply the regression formulas.

Re: [R-sig-phylo] brunch

2012-05-30 Thread Graham Slater
Hi Yanthe, It's been a while since I've used these functions, but if I remember correctly, you can do what you ask with crunch but not brunch. What kinds of models can I use in crunch? Can I do Y as a function of X1 and X2, all of the variables being continuous? Are the contrasts for