Dear colleagues,
I am testing the impact of categorical binary characters (habitat and
presence/absence of symbionts) on a continuous variable (log of body size)
using PGLS...
I am not sure if I should remove the intercept from the formulae and the
biological interpretation of the absence
Dear Julien,
maybe I dont understand your rmodel...but IF your
model has one continuous dep. and one categorical
(binary) indep. it looks like an ANOVA model: in this
case phy.anova() [or phy.manova() if you have 1
dependents] in the R package geiger does it. IF the
model is different
Hi Jonathan.
The thing to do in this case is to estimate K with within-species
variability. This is described in Ives et al. (2007; Syst. Biol.) and
implemented in the phytools function phylosig. This will give an
unbiased estimate of K. (K estimated when within-species variability is