ThreeWay (version 1.1.3)

splithalfCP: Split-Half Analysis

Description

Performs split-half analysis for Candecomp/Parafac.

Usage

splithalfCP(X, n, m, p, r, centopt, normopt, scaleopt, addanal, conv,  maxit, ort1, ort2, ort3, 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

r

Number of extracted components

centopt

Centering option (see cent3)

normopt

Normalization option (see norm3)

scaleopt

Scaling option (see renormsolCP)

addanal

Number of additional runs

conv

Convergence criterion

maxit

Maximal number of iterations

ort1

Type of constraints on A (see CP)

ort2

Type of constraints on B (see CP)

ort3

Type of constraints on C (see CP)

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

Afull

Component matrix for the A-mode (full data)

As1

Component matrix for the A-mode (split n.1)

As2

Component matrix for the A-mode (split n.2)

Bfull

Component matrix for the B-mode (full data)

Bs1

Component matrix for the B-mode (split n.1)

Bs2

Component matrix for the B-mode (split n.2)

Cfull

Component matrix for the C-mode (full data)

Cs1

Component matrix for the C-mode (split n.1)

Cs2

Component matrix for the C-mode (split n.2)

References

P.M. Kroonenberg (2008). Applied Multiway Data Analysis. Wiley, New Jersey.

See Also

CP

Examples

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: 
# # Split-half analysis on CP solution
# splitCP <- splithalfCP(TVdata, 30, 16, 15, 2, 0, 0, 0, 5, 1e-6, 10000, 1, 1, 1, 
#  labSTUDENT, labSCALE, labPROGRAM)
# # Split-half analysis on CP solution (when labels are not available)
# splitCP <- splithalfCP(TVdata, 30, 16, 15, 2, 0, 0, 0, 5, 1e-6, 10000, 1, 1, 1)
# ## End(Not run)