Factors = PARAFAC(X,1);
Where X is the data and 1 is the number of factors to extract. Remember that in PARAFAC the models are not nested, so if a two-component model is to be estimated one can not use the one-component model as a basis as in ordinary two-way PCA.