Re: [R-sig-phylo] Normality requirement for assessment of lambda withphylosig (phytools) and fitContinuous (geiger)

2012-04-25 Thread Enrico Rezende
Complementing Ted's comment, the randomization test in Blomberg et al. (2003) and the K statistics is available the package picante. Now, caution is warranted because having many species with the same phenotypic value in the dataset (i.e., many zeros) can dramatically decrease the statistical

Re: [R-sig-phylo] Comparing means under a phylogenetic perspective?

2011-03-23 Thread Enrico Rezende
Rafael, a phylogenetic paired t-test should do the job. Please check this reference: Lindenfors P, Revell LJ, and Nunn CL. Sexual dimorphism in primate aerobic capacity: a phylogenetic test. J Evolution Biol 23: 1183-1194, 2010. Liam, any comments on this? Enrico Rafael Maia escribió:

Re: [R-sig-phylo] Dealing with Bounded Trait Measures

2011-03-05 Thread Enrico Rezende
David, on the top of my head, if no species measurement strictly corresponds to zero, you may log-transform the data. You may then simulate Brownian motion in log-transformed values, which will correspond to a boundary of zero in a linear scale (i.e., the more negative the log number, the

Re: [R-sig-phylo] K for discrete unordered traits

2011-03-02 Thread Enrico Rezende
Hi Stacey, you can try this function, which is analogous to Blomberg's K for discrete traits and follows Maddison Slatkin 1991. An example is available at the bottom of the syntax. Cheers, Enrico ## Function phylo.signal.disc

Re: [R-sig-phylo] PICs in non-ultrametric trees

2011-02-02 Thread Enrico Rezende
Hi Peter, the non-standardized contrasts tells you how much sister lineages differ phenotypically (i.e., phenotype lineage A - phenotype lineage B) during the computation of contrasts (assumptions of brownian motion, etc). When you standardize them, you are dividing this difference by the

Re: [R-sig-phylo] fitContinuous running for days

2010-11-15 Thread Enrico Rezende
Manuela and David, perhaps you want to try the OU-transformation syntax I wrote to translate the MATLAB code developed by Blomberg, Ives Garland. It works well with small phylogenies, I never tried it in large ones and am actually curious to see what would be the outcome. You will find the

Re: [R-sig-phylo] PIC vs. PGLS

2010-09-27 Thread Enrico Rezende
Hi everybody, I am able replicate the results of independent contrasts and PGLS employing R (mapping the intercept to the phylogenetic means). I also recall comparing results of PDAP employing a star phylogeny against a regular regression in SPSS for Windows. And as far as I can tell, this

Re: [R-sig-phylo] Simulate continuous character evolution on phylogenetictrees using different models of evolution

2009-11-06 Thread Enrico Rezende
Hi Christian, try the function shown below (an example follows). Everybody else is welcome to give it a try and diagnose if everything works fine (which seems to be the case as far as I am aware). Cheers, Enrico 'cor.brown' -

Re: [R-sig-phylo] test correlated evolution using two phylogenies?(consumerresource)

2009-06-07 Thread Enrico Rezende
Hi Mick, I am not entirely sure of the best approach to analyze your data, but here are a few comments (IN CAPS): Dear r sig phylo users, I am doing a phylogenetic comparative analysis of (quantitative) foraging behaviour and body size ratio between consumers resources, and am not sure how