Library: | stats |
See also: | factoranalysis varimax varimaxval factorrot |
Quantlet: | factorlikehood | |
Description: | auxiliary quantlet for factoranalysis.xpl, which computes the likelihood function |
Usage: | like = faclikehood(x) | |
Input: | ||
x | n x 1 vector or n x c matrix. Each column represents one vector of specifics. | |
Output: | ||
like | scalar or 1 x c vector as results of the likelihood function for the columns of x |
library("nummath") library("xplore") library("stats") y=read("geopol") y=y[,2:11] b=factoranalysis(y,4,"MLM",10) PSI=b.specifics factorlikehood(log(PSI./(1-PSI)))
The minimized likelihood of the FA model Contents of like [1,] 0.41388