gpa_fit: GPA

In psolymos/EDMAinR: Euclidean Distance Matrix Analysis in R

GPA

Description

Fit GPA or WGPA to landmark data.

Usage

gpa_fit(x, B = 0, ncores = getOption("Ncpus", 1L),
    weighted=FALSE, ...)
## S3 method for class 'gpa_fit'
print(x, truncate=40, ...)

Arguments

x	an EDMA data object of class 'edma_data'.
B	nonnegative integer, the number of bootstrap replicates.
weighted	logical, use shapes::procWGPA instead of shapes::procGPA.
ncores	positive integer, the number of cores to use when bootstrapping. Use options(Ncpus = 2) to set it to 2 globally.
truncate	numeric, number of characters to print for the object title.
...	arguments passed to shapes::procGPA or shapes::procWGPA.

Value

Returns only form matrix, SigmaKstar is NA.

Author(s)

Peter Solymos wrote the wrapper for shapes::procGPA.

References

Gower, J.C. (1975). Generalized Procrustes analysis, Psychometrika, 40, 33–50.

Examples

file <- system.file(
    "extdata/crouzon/Crouzon_P0_Global_MUT.xyz",
    package="EDMAinR")
x <- read_xyz(file)
x <- x[,,1:10] # 10 specimens

## nonparametric fit
fit <- gpa_fit(x, B=9)
fit
str(Meanform(fit))
str(SigmaKstar(fit))

psolymos/EDMAinR documentation built on Aug. 25, 2023, 7:54 p.m.