| 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