Learn R Programming

FactoMineR (version 1.01)

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, 
    coord=c(1,2), group, name.group = NULL, graph = TRUE)

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
coord
a length 2 vector specifying the components to plot
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

Value

  • A list containing the following components:
  • RVa matrix of RV coefficients between partial configurations
  • RVsa matrix of standardized RV coefficients between partial configurations
  • simia matrix of Procrustes similarity indexes between partial configurations
  • scalinga vector of isotropic scaling factors
  • depan array of initial partial configurations
  • consensusa matrix of consensus configuration
  • Xfinan array of partial configurations after transformations
  • correlationscorrelation matrix between initial partial configurations and consensus dimensions
  • PANOVAa list of "Procrustes Analysis of Variance" tables, per assesor (config), per product(objet), per dimension (dimension)

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

Run this code
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
plot.GPApartial (res.gpa)

Run the code above in your browser using DataLab