Library: | stats |
See also: | factoranalysis factorlikehood varimax factorrot |
Quantlet: | varimaxval | |
Description: | auxiliary quantlet for varimax, it calculates the value of the varimax negative criterion |
Usage: | V = varimaxval(x) | |
Input: | ||
x | n x 1 vector or n x c matrix where each column represents a set of all angles for rotation of a loadings matrix for evaluation of varimax criterion. n is a number of angles needed for a unique identification of some rotation. In k-dimensional space, to obtain k x k rotation matrix we need k*(k-1)/2 angles. | |
Output: | ||
V | scalar or 1 x c vector of resulting values of negative varimax. |
library("nummath") library("xplore") library("stats") y=read("decathlo") y=y[,1:10] b=factoranalysis(y,3,"PFM",10) globalQvarimax=b.loadings putglobal("globalQvarimax") k=cols(b.loadings) count=k*(k-1)/2 init=0*matrix(count) varimaxval(init) rot=varimax(b.loadings) globalQvarimax=rot putglobal("globalQvarimax") varimaxval(init)
The value of the varimax criterion before and after the varimax rotation: Contents of V [1,] -0.51824 Contents of V [1,] -0.83458