[R] factanal question
William Revelle
lists at revelle.net
Tue Dec 2 18:35:48 CET 2008
At 11:04 AM -0500 12/1/08, John Fox wrote:
>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?
John,
Good suggestion. I will add that in the next revision.
Bill
>
>Regards,
> John
>
>------------------------------
>John Fox, Professor
>Department of Sociology
>McMaster University
>Hamilton, Ontario, Canada
>web: socserv.mcmaster.ca/jfox
>
>
>> -----Original Message-----
>> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
>On
>> Behalf Of William Revelle
>> Sent: December-01-08 10:26 AM
>> To: John Fox; 'Don McNeil'
>> Cc: r-help at 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: r-help-bounces at r-project.org
>[mailto:r-help-bounces at r-project.org]
>> >On
>> >> Behalf Of Don McNeil
>> >> Sent: November-30-08 11:39 PM
>> >> To: r-help at 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 at 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 at 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 at 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/
More information about the R-help
mailing list