# pcamean: PCA of the mean matrix In ThreeWay: Three-Way Component Analysis

## Description

Performs Principal Component Analysis (PCA) of the mean matrix aggregated over mode number indicated by `aggregmode`.

## Usage

 `1` ```pcamean(X, n, m, p, laba, labb, labc) ```

## 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 `laba` Optional vector of length `n` containing the labels of the `A`-mode entities `labb` Optional vector of length `m` containing the labels of the `B`-mode entities `labc` Optional vector of length `p` containing the labels of the `C`-mode entities

## Value

A list including the following components:

 `Y` An object of class `matrix` containing the mean matrix `ev` A vector containing the eigenvalues of `Y` `A1` Component matrix for the `A` mode based on varimax rotation of loadings `B1` Component matrix for the `B` mode based on varimax rotation of loadings `C1` Component matrix for the `C` mode based on varimax rotation of loadings `A2` Component matrix for the `A` mode based on oblique ‘HKIC’ (Harris-Kaiser Independent Cluster) orthomax rotation of loadings `B2` Component matrix for the `B` mode based on oblique ‘HKIC’ (Harris-Kaiser Independent Cluster) orthomax rotation of loadings `C2` Component matrix for the `C` mode based on oblique ‘HKIC’ (Harris-Kaiser Independent Cluster) orthomax rotation of loadings

## Note

`aggregmode` denotes the mode over which means are computed (1 for `A`-mode, 2 for `B`-mode, 3 for `C`-mode).
`aggregmode` is provided interactively.

## Author(s)

Maria Antonietta Del Ferraro [email protected]
Henk A.L. Kiers [email protected]
Paolo Giordani [email protected]

## References

H. Kaiser (1958). The varimax criterion for analytic rotation in factor analysis. Psychometrika 23:187–200.
C. Harris \& H. Kaiser (1964). Some mathematical notes on three-mode factor analysis. Psychometrika 29:347–362.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15``` ```data(TV) TVdata=TV[[1]] labSCALE=TV[[2]] labPROGRAM=TV[[3]] labSTUDENT=TV[[4]] # permutation of the modes so that the A-mode refers to students TVdata <- permnew(TVdata, 16, 15, 30) TVdata <- permnew(TVdata, 15, 30, 16) ## Not run: # PCA on the mean matrix TVpcamean <- pcamean(TVdata, 30, 16, 15, labSTUDENT, labSCALE, labPROGRAM) # PCA on the mean matrix (when labels are not available) TVpcamean <- pcamean(TVdata, 30, 16, 15) ## End(Not run) ```

ThreeWay documentation built on May 29, 2017, 11:52 p.m.