[R-sig-phylo] PGLS, categorical data and regression through origin

2012-06-06 Thread Julien Lorion
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

Re: [R-sig-phylo] PGLS, categorical data and regression through origin

2012-06-06 Thread ppiras
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

Re: [R-sig-phylo] Question about Kcalc with 1 data points per species

2012-06-06 Thread Liam J. Revell
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