Aufgaben vollständig gelöst
-
Sechs Punkte mit 99,49% der Varianz auf PC1 und 0,51% auf PC2 8 Schritte
Sechs Punkte: (2, 1), (3, 2), (4, 2,2), (5, 3,2), (6, 3,9), (7, 5). Das Bedienfeld weist 99,49% der Varianz auf PC1 und 0,51% auf PC2 aus, was wie eine Punktwolke erscheint, die im Wesentlichen eine Gerade ist. Berechnen Sie, wie weit die Punkte tatsächlich von dieser Geraden entfernt liegen.
-
Zuerst zentrieren, da es im Folgenden nur um die Streuung und nicht um die Position geht. Die x-Werte verlaufen von 2 bis 7 in Einzelschritten, sodass ihr Mittelwert genau in der Mitte des Bereichs liegt und die zentrierte x-Spalte symmetrisch ist: ±0,5, ±1,5, ±2,5. Nur für y ist Rechenaufwand erforderlich.
-
Bei der Kovarianz werden die summierten Produkte durch n − 1 = 5 dividiert, nicht durch 6. Ein Freiheitsgrad wurde bereits für die Schätzung des gerade subtrahierten Mittelwerts verbraucht, und bei nur sechs Punkten macht diese Korrektur 20% aus.
-
Zwei Zahlen genügen, um beide Eigenwerte zu erhalten, und keine von beiden erfordert eine Richtung. Die Spur ist die Gesamtvarianz; die Determinante misst, wie weit die Matrix von der Singularität entfernt ist. Gegenüber einer Spur von 5,581667 entspricht eine Determinante von 0,156933 einer Matrix, die sehr nahe an Rang 1 liegt.
-
Lösen Sie die charakteristische quadratische Gleichung. λ1 ergibt sich als etwa das 200-Fache von λ2 — die Form dieser Punktwolke verdichtet in einem einzigen Vergleich.
-
Die Zeilen für die erklärte Varianz stellen genau dieses Verhältnis dar und nichts anderes. 99,49% ist somit eine Aussage über λ1 und λ2 und erbt deren Natur: Varianzen, die die Einheiten der Daten im Quadrat tragen. Behalten Sie diesen Gedanken im Hinterkopf.
-
Die Richtung ist der Eigenvektor, und für eine 2×2-Kovarianzmatrix mit einer Nebendiagonale ungleich null kann man ihn ohne weiteres Lösen direkt hinschreiben: (λ1 − Syy, Sxy) zeigt bereits entlang von PC1. Skaliert auf Einheitslänge ergibt sich der ausgegebene Vektor.
-
Drehen Sie diesen Vektor um eine Vierteldrehung, um v2 zu erhalten, und projizieren Sie den ersten zentrierten Punkt darauf. Das Ergebnis, 0,031179, ist die zweite Zahl in der Score-Tabelle. Diese gesamte Spalte — die das Bedienfeld ausgibt, ohne je zu erklären, was sie bedeutet — ist also der vorzeichenbehaftete Abstand jedes Punkts von der PC1-Geraden, und der schlechteste der sechs Werte beträgt −0,236855.
-
Nun zur Frage, die von den Prozentangaben verborgen wurde. λ2 ist eine Varianz, also ein quadratischer Abstand, und die typische Abweichung ist ihre Quadratwurzel. Vergleichen Sie das mit der Quadratwurzel aus λ1, welche die Streuung entlang der Geraden darstellt.
Antwort
Die Punkte liegen etwa 7,13% der Länge der Wolke von der PC1-Geraden entfernt, nicht 0,51%. Beide Zahlen beschreiben dieselben sechs Punkte; die eine ist ein Verhältnis von Quadraten, die andere ein Verhältnis von Längen, und die Quadratwurzel dazwischen macht hier einen Faktor von 14 aus. Diese Lücke ist die Falle in jeder Regel der Form „behalte 99% der Varianz“: Um die rechtwinklige Abweichung unter 1% der Streuung zu halten, muss PC2 unter 0,01% der Varianz liegen — einhundertmal strenger, als die Regel klingt. In Dateneinheiten beträgt die Abweichung √λ2 = 0,168104, und sie ist nicht verborgen — sie ist im Diagramm als graues Segment von jedem grünen Punkt zu seiner blauen Projektion eingezeichnet, und diese Segmente haben offensichtlich keine Länge von null.
-
-
PC1 mit Steigung 0,769 und eine Anpassung der kleinsten Quadrate von 0,763 6 Schritte
PC1 wird als (0,792687, 0,609629) ausgegeben, eine Richtung mit der Steigung 0,769. Die Anpassung von y an x nach der Methode der kleinsten Quadrate aus genau denselben Kovarianzeinträgen ergibt 2,67/3,5 = 0,763. Zwei Geraden durch eine Wolke aus sechs Punkten. Zeigen Sie, welche davon die steilere ist — und dass es immer diese eine ist.
-
Die Methode der kleinsten Quadrate bezüglich y minimiert die vertikalen Abstände, und der Minimierer ist Sxy/Sxx. Beachten Sie, was im Nenner steht: die Streuung von x allein. Die Anpassung nimmt x als gegeben an und fragt nur, was y als Reaktion darauf tut.
-
Die Steigung von PC1 wird direkt aus dem ausgegebenen Einheitsvektor abgelesen, da eine Richtung (a, b) die Steigung b/a hat. Sie liegt etwa 6 Tausendstel über der Regressionssteigung — klein, aber nicht null.
-
Um die beiden richtig zu vergleichen, schreiben Sie die Steigung von PC1 so auf, wie sie durch den Eigenvektor definiert ist, statt als Verhältnis von Dezimalzahlen. Der Zähler ist bei beiden dasselbe Sxy, sodass der gesamte Unterschied im Nenner liegt: Sxx für die Regression, λ1 − Syy für PC1.
-
Drücken Sie λ1 aus. Ob PCA den kleineren Nenner hat, ist nun zu einer Frage nach einer einzelnen Quadratwurzel geworden.
-
Und diese Quadratwurzel ist genau dann kleiner als die Spur, wenn die Determinante positiv ist, was für jede Punktwolke gilt, die keine perfekte Gerade ist. Der Nenner von PC1 ist daher immer der kleinere, sodass PC1 immer die steilere der beiden ist. Nichts an diesen sechs Punkten wurde hierfür verwendet.
-
Tauschen Sie die Rollen und führen Sie stattdessen eine Regression von x auf y durch. Das ist eine dritte Gerade, und wenn man ihre Steigung für dasselbe y-gegen-x-Bild umschreibt, liegt sie bei 0,779650. PC1 verändert sich bei diesem Tausch nicht — die Kovarianzmatrix ist symmetrisch in x und y —, sodass das Argument des vorherigen Schritts nun in die andere Richtung verläuft und PC1 die flachere dieses Paares ist.
Antwort
Drei Geraden durch eine Wolke: 0,762857 < 0,769067 < 0,779650, wobei PC1 strikt in der Mitte liegt. Die beiden Regressionen weichen um 2,2% bei einer Wolke voneinander ab, die das Bedienfeld als zu 99,49% eindimensional bezeichnet, und sie weichen voneinander ab, weil sie unterschiedliche Fragen beantworten — sage y aus x voraus oder sage x aus y voraus. PC1 beantwortet keine von beiden. Sie minimiert den rechtwinkligen Abstand, sodass eine Umbenennung der Achsen das Bild dreht und die Gerade genau dort belässt, wo sie war, und ihre Residuen sind die zweite Score-Spalte statt eines vertikalen Abstands. Die drei fallen nur dann zu einer einzigen zusammen, wenn det C = 0 ist, das heißt, wenn die Punkte wirklich kollinear sind und es nichts mehr gibt, worüber man uneinig sein könnte. Zwischen diesen Extremen ist PCA die geeignete Gerade, wenn beide Koordinaten einen Messfehler tragen, und die falsche Gerade, wenn Sie tatsächlich eine Vorhersage wünschen: Die Regressionssteigung ist absichtlich flacher, weil das Hinziehen zur Waagerechten den quadratischen Vorhersagefehler am kleinsten macht.
-
Lernpfad
Matrizen, die Dinge bewegen
Quellen (2)
- Insight block 3 — the perpendicular line of closest fit, thirty years before the name: K. Pearson, "On lines and planes of closest fit to systems of points in space." Philosophical Magazine 2(11), 559–572, 1901.
- Where "principal component" was coined: H. Hotelling, "Analysis of a complex of statistical variables into principal components." Journal of Educational Psychology 24(6), 417–441, 1933.