# NOT RUN {
## What do the data look like
head(GarlandJanis.Original)
head(GarlandJanis.varcov)
## An example of a GLS fit
fit.gls.GJ <- with(GarlandJanis.Original,
phylog.gls.fit(cbind(body.mass,hind.l.length),
running.speed, GarlandJanis.varcov)
)
summary(fit.gls.GJ) # summary of the gls model; same as with IC
# }
# NOT RUN {
# This data set can be prepared from the original pdi files
# (in directory Examples) as:
GarlandJanis.Orig <- read.pdi.data(c("49ms.pdi","49hmt.pdi"),
variable.names = c("body.mass", "running.speed",
"hind.l.length","mtf.ratio"))
Garland.Janis.Orig$clade <- as.factor(c(rep("Carnivore",19),
rep("Herbivore",30)))
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab