# NOT RUN {
# a GLS fit
data(Lacertid.varcov)
data(Lacertid.Original)
ex.gls.phylog <-
phylog.gls.fit(Lacertid.Original$svl,Lacertid.Original$clutch.size,Lacertid.varcov)
ex.gls.phylog
# }
# NOT RUN {
# This data set can also be obtained from the pdi files
# (see example in GarlandJanis.Original), or as:
LacertidSim <- read.sim.data(c("ifsmi.sim","ihshw.sim","iclag.sim","icfxx.sim"),
pdi.files=c("ifsmi.pdi","ihshw.pdi","iclag.pdi", "icfxx.pdi"),
variable.names = c("svl","svl.matur","hatsvl","hatweight",
"clutch.size", "age.mat","cl.freq", "xx"))
LacertidSim <- LacertidSim[,-10]
LacertidOriginal <- LacertidSim[LacertidSim$sim.counter==0,-1]
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab