[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

From |
"Cathy L. Antonakos" <cathya@umich.edu> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: rotated components in pca, order of |

Date |
Thu, 5 Jul 2007 12:57:59 -0400 (EDT) |

Okay. that's great. Thanks I understand now! Cathy On Thu, 5 Jul 2007, khigbee@stata.com wrote:

Cathy L. Antonakos <cathya@umich.edu> asks:I'm running a principal components analysis in Stata, and am puzzled by the order of the components in the rotated solution. Prior to rotation, the components are listed as "comp 1" "comp 2" "comp 3" "comp 4" in the table of component loadings. After rot ation, the order of the components changes to: 3, 2, 1, 4. The loadings also change, which I expected. I can figure out which of the rotated components is what, based on running a correlation between the rotated component scores and selected variables that I expect to load highly on particular components. But I don't know why the order of the components has changed. I did look in Statalist and the manual and found references on how to generate component scores, but nothing more. My question is - are the rotated components reordered because due to rotation the percentage of variance explained by a particular component has changed? After rotation, in this example, does component 3 now explain more variance than component 1 did in the unrotated solution?If you had Stata 10 instead of an earlier version of Stata, the rotated components would be labeled in the output from 1 to 4 instead of 3, 2, 1, 4. The issue is only cosmetic--the rotated components are the same either way. Rotation is creating new components that are linear combinations of (all of) the original components. Afterwards the resulting columns are reordered from highest to lowest percent of variance explained (just as you guessed in your question above). Previous to Stata 10, when we reordered the rotated components, the column labels ended up reordered at the same time. There is no good reason to do that, and so with Stata 10 the columns are labeled from 1 up.I'm also wondering if anyone can explain what the "component rotation matrix" is. I've seen this in factor analysis output as "factor rotation matrix" but despite quite a bit of searching, I haven't been able to find information about what this is.Maybe the best way to understand the "factor rotation matrix" (found in the e(r_T) returned matrix result) is by example. Borrowing an example from "[MV] pca postestimation" we have . webuse audiometric . pca l* r*, comp(4) . rotate . mat L = e(L) . mat rL = e(r_L) . mat rT = e(r_T) matrix L has the unrotated components, matrix rL has the rotated components, and rT has the component rotation matrix. These matrices are related by rL = L * rT You can verify by . mat chk = L * rT . mat list rL . mat list chk Ken Higbee khigbee@stata.com StataCorp 1-800-STATAPC * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

* * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**References**:**Re: st: rotated components in pca, order of***From:*khigbee@stata.com

- Prev by Date:
**st: RE: do-file editor - difference between Stata 9 and 10** - Next by Date:
**st: generate under bylist** - Previous by thread:
**Re: st: rotated components in pca, order of** - Next by thread:
**st: RE: SSC & search - thanks** - Index(es):

© Copyright 1996–2016 StataCorp LP | Terms of use | Privacy | Contact us | What's new | Site index |