Dear Bill, Thanks for pointing out that this functionality is already in the psych package. Shouldn't factor.residuals() avoid this computation for oblique rotations?
Regards, John ------------------------------ John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web: socserv.mcmaster.ca/jfox > -----Original Message----- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On > Behalf Of William Revelle > Sent: December-01-08 10:26 AM > To: John Fox; 'Don McNeil' > Cc: r-help@r-project.org > Subject: Re: [R] factanal question > > Don and John, > factor.residuals in the psych package does what you want (and > basically what John wrote). > > Bill > > > At 9:30 AM -0500 12/1/08, John Fox wrote: > >Dear Don, > > > >All long as you leave the factors unrotated or do an orthogonal rotation (as > >is the default), you can compute reproduced correlations among the variables > >from the factor loadings, and thus residual correlations given the loadings > >and the original correlation matrix, both of which are accessible in the > >object returned by factanal(); the following isn't carefully tested, but > >should work: > > > >repRes <- function(F, round=3){ > > A <- loadings(F) > > R <- F$correlation > > RR <- A %*% t(A) > > ResR <- R - RR > > list(reproduced.correlations=round(RR, round), > > residual.correlations=round(ResR, round)) > > } > > > >Here F is an object returned by factanal(). The diagonal elements of the > >reproduced correlations are the communalities, and of the residual > >correlations, the uniquenesses. > > > >To do this from an oblique rotation would require the factor-correlation > >matrix, which, as has been pointed out previously, factanal() oddly doesn't > >provide. In this case, that's not a real impediment, since reproduced and > >residual correlations are invariant with respect to rotation of the factors. > > > >I hope this helps, > > John > > > >------------------------------ > >John Fox, Professor > >Department of Sociology > >McMaster University > >Hamilton, Ontario, Canada > >web: socserv.mcmaster.ca/jfox > > > > > >> -----Original Message----- > >> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > >On > >> Behalf Of Don McNeil > >> Sent: November-30-08 11:39 PM > >> To: r-help@r-project.org > >> Subject: [R] factanal question > >> > >> Dear R users: > >> I'm wondering if it's possible to get the residual correlation matrix > when > >> using factanal. > >> Since factanal assumes that the errors are normally distributed and > >> independent (provided the factor model fits the data) this would be > >useful. > >> Of course you would need to submit the data to the function to get the > >> residuals (not just their correlation matrix), but it should be possible > >to > >> get the residual correlation matrix if only the data correlation matrix > is > >> provided. > >> Don McNeil > >> > >> ______________________________________________ > >> R-help@r-project.org mailing list > >> https://stat.ethz.ch/mailman/listinfo/r-help > >> PLEASE do read the posting guide > >http://www.R-project.org/posting-guide.html > >> and provide commented, minimal, self-contained, reproducible code. > > > >______________________________________________ > >R-help@r-project.org mailing list > >https://stat.ethz.ch/mailman/listinfo/r-help > >PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > >and provide commented, minimal, self-contained, reproducible code. > > > -- > William Revelle http://personality-project.org/revelle.html > Professor http://personality-project.org/personality.html > Department of Psychology http://www.wcas.northwestern.edu/psych/ > Northwestern University http://www.northwestern.edu/ > Attend ISSID/ARP:2009 http://issid.org/issid.2009/ > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
