# 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 obtained from the original dsc file as:
GarlandJanis.varcov <- read.phylog.matrix("49ms.dsc")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab