| align.procrustes | R Documentation |
GPA is performed by a call to gpagen.
If outlier.analysis = TRUE then the function will display the results of
plotOutliers, then warp grids for the most extreme shapes.
The user is prompted to select a certain number of extreme shapes to exclude from
a repeated Procrustes alignment. This processes is iterative, until the user indicates that
it should stop or no shapes are requested for removal.
align.procrustes( A, curves = NULL, show.plot.gpa = TRUE, outlier.analysis = FALSE, display.outliers = 4, output.gdf = TRUE, ... )
A |
A list or 3-dimensional array containing XY shape coordinates for multiple specimens. |
curves |
An optional matrix defining which landmarks should be treated as semilandmarks by |
show.plot.gpa |
A logical value specifying whether to plot the Procrustes-aligned landmarks. |
outlier.analysis |
A logical value specifying whether to perform interactive outlier analysis. Alternatively, a numeric vector can be provided indicating the index positions of specimens to remove. A character vector of specimen names can also be provided, indicating the specimens to remove. |
display.outliers |
The number of outlier shapes to display in the interactive outlier analysis. |
output.gdf |
A logical value specifying whether to return output with
a list structure that includes an element, |
Returns a list with Procrustes-aligned coordinates in coords, centroid size in Csize.
All other elements from gpagen are in the the element gpagen.
Any other list elements in the input are retained.
Data provenance is also retained in provenance and updated with details of the GPA..
Dave Angelini david.r.angelini@gmail.com [aut, cre]
Adams and M.L. Collyer and A. Kaliontzopoulou. 2020. Geomorph: Software for geometric morphometric analyses. R package version 3.2.1. D.C. (Link)
data("Bombus.forewings", package = "borealis")
fw.gpa <- align.procrustes(Bombus.forewings)
fw.gpa <- align.procrustes(Bombus.forewings, show.plot.gpa = FALSE)
plot(fw.gpa$gpagen)
fw.gpa <- align.procrustes(Bombus.forewings, outlier.analysis = TRUE)
fw.gpa <- align.procrustes(Bombus.forewings, outlier.analysis = 57)
fw.gpa <- align.procrustes(Bombus.forewings,
outlier.analysis = c("FJ190828-002","DRA190718-001"),
print.progress = FALSE )
names(fw.gpa$provenance)
cat(fw.gpa$provenance$align.procrustes)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.