> Ef = X - Mf;
Use these residuals for plotting purposes as was done with the raw data in the introductory chapter. Use the residuals for calculating the sum-squared residuals (or obtain these directly in the output from the algorithm) for scree-plots.
For calculating the core consistency use the m-file corcondia. For the model with loadings in the vector Factorf calculate the core consistency as
> cf = corcond(X,Factorf);
where f is the number of components. Compare the core consistencies for various numbers of components:
> plot([1:5],[c1 c2 c3 c4 c5])
> xlabel('Number of components')
> ylabel('Core consistency')