Build data and AI skills | 50% off

Last chance! 50% off unlimited learning

Sale ends in


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

Run this code
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)

Run the code above in your browser using DataLab