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
procD.lm(coords ~ Csize + species * site, data = gdf)
Run the code above in your browser using DataLab