This function is used in the main function CA3variants when the input parameter is catype="NSCA3".
It decomposes the Marcotorchino index, computes principal axes, coordinates, weights of rows and columns,
total inertia (equal to the Marcotorchino index) and the rank of the matrix.
nsca3basic(x, p, q, r, test =10^-6, ctr = T, std = T, sign =TRUE)
Arguments
x: The three-way contingency table.
p: The number of components of the first mode.
q: The number of components of the second mode.
r: The number of components of the third mode.
test: The treshold used in the algorithm.
ctr: The flag parameter to center the data (T or F), if F the data are not centered.
std: The flag parameter to weight the data (T or F), if F the data are not weighted.
sign: The input parameter for changing the sign to the components according to the core sign.
Returns
x: The original three-way contingency table.
xs: The weighted three-way contingency table.
xhat: The three-way contingency table reconstructed after Tuckals3 by means of the principal components and core array.
nxhat2: The inertia of the three-way non-symmetrical correspondence analysis
for one response (the three-way Marcotorchino index).
prp: The proportion of inertia reconstructed using the principal components and the core array to the total inertia.
To select the model dimensions (number of principal components), we examine the inertia explained by the p, q, r principal components with respect to the overall fit.
a: The row principal components.
b: The column principal components.
cc: The tube principal components.
g: The core array (generalized singular values) calculated by using the Tuckals3 algorithm.
They help to explain the strength of the association among the three principal components.
iteration: The number of iterations that are required for the TUCKALS3 algorithm to converge.
References
Beh EJ and Lombardo R (2014) Correspondence Analysis, Theory, Practice and New Strategies. John Wiley & Sons.