Re: [R] how to get residuals in factanal

2008-01-07 Thread Yijun Zhao
1. ?factanal says: == The factor analysis model is x = Lambda f + e ... but the major restriction is that the scores be uncorrelated and of unit variance, and that the errors be independent with variances Phi, the uniquenesses. ===

Re: [R] how to get residuals in factanal

2008-01-06 Thread Prof Brian Ripley
On Sun, 6 Jan 2008, Yijun Zhao wrote: p.s. I tried to use both "regression" and "Bartlett" way to get the scores. In both cases, the scores are uncorrelated, but the errors are NOT uncorrelated to the scores, and are also NOT uncorrelated among themselves. What am I missing? factanal() are su

Re: [R] how to get residuals in factanal

2008-01-06 Thread Yijun Zhao
p.s. I tried to use both "regression" and "Bartlett" way to get the scores. In both cases, the scores are uncorrelated, but the errors are NOT uncorrelated to the scores, and are also NOT uncorrelated among themselves. What am I missing? factanal() are supposed to give independent errors vecto

Re: [R] how to get residuals in factanal

2008-01-06 Thread Yijun Zhao
The factanal was called with 'varimax' rotation. The factors scores are uncorrelated. But the residuals I got by using X - sum(loadings*factors-scores) is not uncorrelated to the factor scores. I thought the residuals should be independent to the factor scores as ?factanal says: =

Re: [R] how to get residuals in factanal

2008-01-06 Thread Prof Brian Ripley
On Sun, 6 Jan 2008, Yijun Zhao wrote: > In R factanal output, I can't find a function to give me residuals e. > > I mannually got it by using x -lamda1*f1 -lamda2*f2 - ... -lamdan*fn, but > the e > I got are not uncorrelated with all the f's. > > What did I do wrong? Please help. What did you u

[R] how to get residuals in factanal

2008-01-06 Thread Yijun Zhao
In R factanal output, I can't find a function to give me residuals e. I mannually got it by using x -lamda1*f1 -lamda2*f2 - ... -lamdan*fn, but the e I got are not uncorrelated with all the f's. What did I do wrong? Please help. Yijun __