Calculates Bro and Kiers's core consistency diagnostic (CORCONDIA) for a fit parafac or parafac2 model. For Parafac2, the diagnostic is calculated after transforming the data.

Usage

1

Arguments

X

Three-way data array with dim=c(I,J,K) or four-way data array with dim=c(I,J,K,L). Can also input a list of two-way or three-way arrays (for Parafac2).

object

Object of class "parafac" (output from parafac) or class "parafac2" (output from parafac2).

divisor

Divide by number of factors (default) or core sum of squares.

Details

The core consistency diagnostic is defined as

100 * ( 1 - sum( (G-S)^2 ) / divisor )

where G is the least squares estimate of the Tucker core array, S is a super-diagonal core array, and divisor is the sum of squares of either S ("nfac") or G ("core"). A value of 100 indiciates a perfect multilinear structure, and smaller values indicate greater violations of multilinear structure.