Computation of fit contributions by combinations of modes in case of renormalization .
T3fitpartitioning(Xprep, n, m, p, AS, BT, CU, K, renormmode, laba, labb, labc)
Arguments
Xprep: Matrix (or data.frame coerced to a matrix) of order (nxmp) containing the matricized array (frontal slices)
n: Number of A-mode entities
m: Number of B-mode entities
p: Number of C-mode entities
AS: Component matrix for the A-mode
BT: Component matrix for the B-mode
CU: Component matrix for the C-mode
K: Matricized core array (frontal slices)
renormmode: Renormalization option (0 for no renormalization, 1 for fit contribution to total fit of each B- and C-mode component combination, 2 for fit contribution to total fit of each A- and C-mode component combination, 3 for fit contribution to total fit of each A- and B-mode component combination)
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
Returns
A list including the following components: - fitA: Fit contribution for the A-mode entities
fitB: Fit contribution for the B-mode entities
fitC: Fit contribution for the C-mode entities
ABcontr: Contribution to the goodness of fit contributions by combinations of A- and B-modes in case of renormalization
BCcontr: Contribution to the goodness of fit contributions by combinations of B- and C-modes in case of renormalization
ACcontr: Contribution to the goodness of fit contributions by combinations of A- and C-modes in case of renormalization
Note
The computation of the fit contributions by combinations of modes is done in case of renormalization .
data(Bus)# labels for Bus datalaba <- rownames(Bus)labb <- substr(colnames(Bus)[1:5],1,1)labc <- substr(colnames(Bus)[seq(1,ncol(Bus),5)],3,8)# T3 solutionBusT3 <- T3funcrep(Bus,7,5,37,2,2,2,0,1e-6)# Fitpartitioning of the T3 solutionFitT3 <- T3fitpartitioning(Bus,7,5,37, BusT3$A, BusT3$B, BusT3$C, BusT3$H,0, laba, labb, labc)# Fitpartitioning of the T3 solution (when labels are not available)FitT3 <- T3fitpartitioning(Bus,7,5,37, BusT3$A, BusT3$B, BusT3$C, BusT3$H,0)