if (FALSE) {
data(plethodon)
Y.gpa <- gpagen(plethodon$land, PrinAxes = FALSE)
gdf <- geomorph.data.frame(Y.gpa)
attributes(gdf)
gdf <- geomorph.data.frame(Y.gpa, species = plethodon$species,
site = plethodon$site)
attributes(gdf)
# Using geomorph.data.frame to facilitate analysis
anova(procD.lm(coords ~ Csize + species * site, data = gdf))
}
Run the code above in your browser using DataLab