Thanks everyone, exactly what I needed! Went ahead and used threshBayes
since I only have two characters, and indeed 95% HPD of r was > 0 (ESS >
1,000).
Thanks so much for all the advice!
Sean
On Sat, Apr 9, 2016 at 11:23 AM, Joe Felsenstein
wrote:
> Liam --
>
> Thanks,
Just to clarify Joe's email - Threshml can be called using my package
with Scott Chamberlain called Rphylip (not phytools). It is on CRAN;
however, Rphylip is just a 'wrapper' for PHYLIP, which still needs to be
installed locally. In the future we may change Rphylip so that it is
packaged with
Sean --
... or, if you want to do it *really correctly*, you can use the threshold
model of Sewall Wright for the discrete character and use the MCMC approach
that I proposed in 2012:
Felsenstein, J. 2012. A comparative method for both discrete and
continuous characters using the threshold
Alejandro is correct. You can also do it with phylogenetically independent
contrasts or computer simulations:
Garland Jr., T., P. H. Harvey, and A. R. Ives. 1992. Procedures for the
analysis of comparative data using phylogenetically independent contrasts.
Systematic Biology 41:18–32.
Hello Sean,
If the continuous variable is the “response” and the “independent” variable the
discrete one, you can use PGLS, this would be akin to an ANOVA and you can do
it accounting for phylogenetic non-independence.
Is this what you were after?
Cheers
Alejandro
Hello, I have a two traits, one categorical (binary) and one continuous,
and I want to test for a relationship between them accounting for
phylogenetic signal. I have found a plethora of sources for examining
relationships between multiple categorical traits and many others for
examining multiple