Hello everyone,
I am trying to implement Liam Revell's suggestion on the evaluation of
Pagel's lambda simultaneous to fitting PGLS to minimize the effects of
wrong model selection (OLS vs. PGLS) on species data (i.e. Revell 2010,
Methods in Ecology and Evolution 1: 319-329).
Does anyone
Hi Antigoni.
The other responders are correct, but just to be a little more concrete:
# using ape nlme
library(ape); library(nlme)
# assuming your data are contained in named vectors x y
y-y[names(x)]
fit-gls(y~x,data.frame(x,y),correlation=corPagel(1,tree,fixed=FALSE),method=ML)
# assuming
