pcamean: PCA of the mean matrix

Description Usage Arguments Value Note Author(s) References Examples

View source: R/pcamean.R

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 mariaantonietta.delferraro@yahoo.it
Henk A.L. Kiers h.a.l.kiers@rug.nl
Paolo Giordani paolo.giordani@uniroma1.it

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.

Related to pcamean in ThreeWay...