pcasup1: PCASup Analysis

In ThreeWay: Three-Way Component Analysis

Description

Computes PCASup analysis for the direction concerning the reduced mode.

Usage

 pcasup1(X, n, m, p, model)

Arguments

X	Matrix (or data.frame coerced to a matrix) of order (n x mp) containing the matricized array (frontal slices)
n	Number of A-mode entities
m	Number of B-mode entities
p	Number of C-mode entities
model	Tucker1 model choice (1 for T1-A, 2 for T1-B, 3 for T2-C)

Value

A list including the following components:

A	Matrix of the eingenvectors of the supermatrix containing the frontal slices of the array (A-mode)
B	Matrix of the eingenvectors of the supermatrix containing the horizontal slices of the array (B-mode)
C	Matrix of the eingenvectors of the supermatrix containing the lateral slices of the array (C-mode)
la	Vector of the eigenvalues of the supermatrix containing the frontal slices of the array (A-mode)
lb	Vector of the eigenvalues of the supermatrix containing the horizontal slices of the array (B-mode)
lc	Vector of the eigenvalues of the supermatrix containing the lateral slices of the array (C-mode)

Note

pcasup1 computes the Tucker1 solution.
Cumulative sum of eigenvalues and fits from PCAsup applied to the reduced mode are automatically printed.

Author(s)

Maria Antonietta Del Ferraro [e-mail]
Henk A.L. Kiers [e-mail]
Paolo Giordani [e-mail]

References

H.A.L. Kiers (1991). Hierarchical relations among three-way methods. Psychometrika 56: 449–470.
H.A.L. Kiers (2000). Towards a standardized notation and terminology in multiway analysis. Journal of Chemometrics 14:105–122.
L.R Tucker (1966). Some mathematical notes on three-mode factor analysis. Psychometrika 31: 279–311.

See Also

T1

Examples

data(Bus)
# PCA-sup for T1-B
pcasupBus <- pcasup1(Bus, 7, 5, 37, 2)

Example output

PCASUP: eigenvalues mode B
       Eigenvalue Fit(%)
Comp.1  396973.69  98.68
Comp.2    3403.74  99.52
Comp.3    1251.26  99.83
Comp.4     377.82  99.93
Comp.5     295.50 100.00

ThreeWay documentation built on May 2, 2019, 9:20 a.m.