Fit models for the different number of components
Perhaps decrease the maximal number of iterations to 500 if you are in a hurry, but don’t do that in “real life”. For each fitted model, the model of the data being
(f the number of components), can be obtained using the file
Calculate residuals as
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
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')