FactoMineR (version 1.34)

GPA: Generalised Procrustes Analysis

Description

Performs Generalised Procrustes Analysis (GPA) that takes into account missing values.

Usage

GPA(df, tolerance=10^-10, nbiteration=200, scale=TRUE,  group, name.group = NULL, graph = TRUE, axes = c(1,2))

Arguments

df

a data frame with n rows (individuals) and p columns (quantitative varaibles)

tolerance

a threshold with respect to which the algorithm stops, i.e. when the difference between the GPA loss function at step n and n+1 is less than tolerance

nbiteration

the maximum number of iterations until the algorithm stops

scale

a boolean, if TRUE (which is the default value) scaling is required

group

a vector indicating the number of variables in each group

name.group

a vector indicating the name of the groups (the groups are successively named group.1, group.2 and so on, by default)

graph

boolean, if TRUE a graph is displayed

axes

a length 2 vector specifying the components to plot

Value

A list containing the following components: A list containing the following components:

Details

Performs a Generalised Procrustes Analysis (GPA) that takes into account missing values: some data frames of df may have non described or non evaluated rows, i.e. rows with missing values only. The algorithm used here is the one developed by Commandeur.

References

Commandeur, J.J.F (1991) Matching configurations.DSWO press, Leiden University. Dijksterhuis, G. & Punter, P. (1990) Interpreting generalized procrustes analysis "Analysis of Variance" tables, Food Quality and Preference, 2, 255--265 Gower, J.C (1975) Generalized Procrustes analysis, Psychometrika, 40, 33--50 Kazi-Aoual, F., Hitier, S., Sabatier, R., Lebreton, J.-D., (1995) Refined approximations to permutations tests for multivariate inference. Computational Statistics and Data Analysis, 20, 643--656 Qannari, E.M., MacFie, H.J.H, Courcoux, P. (1999) Performance indices and isotropic scaling factors in sensory profiling, Food Quality and Preference, 10, 17--21

Examples

## Not run: 
# data(wine)
# res.gpa <- GPA(wine[,-(1:2)], group=c(5,3,10,9,2),
#     name.group=c("olf","vis","olfag","gust","ens"))
# 
# ### If you want to construct the partial points for some individuals only
# plotGPApartial (res.gpa)
# ## End(Not run)