# Lade den Datensatz "raq.csv" aus dem Ordner discovr_csv, den man hier
# herunterladen kann: https://www.discovr.rocks/csv/discovr_csv.zip
<- readr::read_csv("data/discovr_csv/raq.csv")
raq.tib ## Rows: 2571 Columns: 24
## ── Column specification ──────────────────────────────────────────────────────────────────
## Delimiter: ","
## chr (1): id
## dbl (23): raq_01, raq_02, raq_03, raq_04, raq_05, raq_06, raq_07, raq_08, raq_09, raq_...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
# Lösche die Variable "id", die ganz vorne im Datensatz steht
<- raq.tib |>
raq_items_tib select(-id)
# Berechne die Korrelationen für alle Variablen mit allen Variablen (items)
<- raq_items_tib |>
raq_cor cor()
# Gebe einen Korrelationsplot mit dem Paket "psych" raus.
#psych::corPlot(raq_cor, upper = FALSE)
# Sortiere die Variablen danach, wie stark sie miteinander korrelieren
<- corrplot::corrMatOrder(raq_cor, order = 'FPC')
order.FPC <- corrplot::corrMatOrder(raq_cor, order = 'hclust')
order.hc
# Speichere die Ordnung als Matrix
<- raq_cor[order.hc, order.hc]
raq_cor.FPC
# Schönere Korrelationsplots gibt es mit dem Paket "corrplot" und dem Befehl "corrplot"
::corrplot(raq_cor.FPC, tl.col='black', tl.cex=.75) corrplot
8 Faktorenanalysen
Auch als (explorative) Faktorenanalyse bekannt.
8.1 Die Folien zur Sitzung
Der Vorlesungsmitschnitt
8.2 Multikollinearität und Dimensionsreduktion
Wenn zwei oder mehr Variablen stark miteinander zusammenhängen, ergibt sich eine hohe Multikollinearität. Bei zwei Variablen ergib sich also keine Punktwolke in einem zweidimensionalen Koordinatensystem, sondern eher eine Reihe Punkte die sehr nahe um eine (Regressions)Gerade liegen. Im Grunde kann man sie statt auf zweidimensionalen Ebene auf einer eindimensionalen Gerade darstellen. Wenn wir zB 5 UV’s hätten, die alle sehr hoch miteinander korrelieren (die r’s > .7), dann können wir diese fünf Dimensionen ohne grosse Verluste auf einer Dimension darstellen. Wir könnten einen Index bauen, der diese fünf UV’s abbildet, ohne das wir viele Informationen verlieren würden (wir hoffen, dass wir nur Messrauschen «verlieren» und kaum substantielle Varianz).
8.3 Indices
Indices sind Variablen, die viele andere zusammenfassen. Das ist bei Aktienindices so und bei Mittelwert- oder Summenindices in der Statistik auch. Die Faktoren als Ergebnis einer explorativen Faktorenanalyse sind auch Indices, die die zugrundeliegenden (latenten) Variablen optimal abbilden, und zwar in dem Sinne, dass sie die Varianz aller Variablen optimal erfassen (also nicht alle gleichberechtigt, sondern die bevorzugt, die untereinander stark korrelieren).
8.4 Faktorenanalyse in Worten
In einem ersten Schritt wird eine Regressionsgerade so in die «Punktwolke» aller Dimensionen gelegt, dass diese Gerade alle Variablen möglichst gut abbildet, also ihre Varianz maximal abbildet. In der Regel bleibt dabei einige Varianz übrig, da die Punkte der Wolke nicht alle beziehungsweise mehrheiltich nicht auf der Geraden liegen. Diese ganze Varianz kann nun eine zweite Gerade haben, die so in die übrigen Punkte gelegt wird, dass sie nicht (steht senkrecht auf der ersten Geraden) oder kaum mit der ersten Gerade korreliert (Winkel gleich oder nahe 90°).
8.5 Explorative Faktorenanalyse (Principal Component Analysis – PCA)
Merkmale von Fällen wie Personen, Gruppen oder Inhalten usw. werden gemessen und als Variablen gespeichert. Die möglichen Merkmale einer zu messenden Population (GG) spannen einen Merkmalsraum mit vielen Merkmalsdimensionen auf. Das kann man sich tatsächlich auch räumlich vorstellen. Wenn zum Beispiel allen Befragten mehrere Fragen zu ihrer Medienzufriedenheit auf einer 5er-Skala gestellt werden, dann kann jeder bei jedem zu bewertenden Medium. Zum Beispiel könnte Axel bei der SRG1 eine 2 vergeben, weil er sie nicht so gut findet, Bernd eine 4 für eine bessere Bewertung usw. Bei der Bewertung von 3+ gibt Axel eine 4 und Bernd eine 3. Die Bewertungen für die SRG1 könnten wir auf eine X-Achse legen und die Bewertungen für 3+ auf eine Y-Achse. Damit hätten wir einen Merkmalsraum von zwei Dimensionen (da klingt es mit dem Raum noch komisch, wird aber einfach immer so genannt). Hätten wir noch ein drittes Medium, wie zum Beispiel RTL+, dann könnte man das auf eine Z-Achse packen, die dann die dritte Dimension wäre und damit auch im alltagssprachlichen ein schöner Merkmalsraum mit drei Merkmalsdimensionen. Allerdings gibt es da noch TV-Ostschweiz, tele bärn, 4+, Pro7, ARD, SRG2 usw. Das heisst, wir haben in der Normalität einen viel grösseren Merkmalsraum mit etlichen Merkmalsdimensionen. Nun kann es interessant sein, ob hinter den unterschiedlichen Bewertungen der Medien unterschiedliche zugrundeliegende Vorlieben stecken. Es könnte doch sein, dass viele Leute aus einem bildungsbürgerlichen Anspruch heraus eher Arte, 3SAT, öffentlich rechtlichen wie die Sender der SRG oder ARD und ZDF besser bewertet als die privaten kommerziellen bzw. regionalen Sender. Andere finden vielleicht generell die öffentlich rechtlichen Sender blöd, weil sie ihnen Staatsnähe unterstellen oder sie Gebühren zahlen müssen. Es könnte also zugrundeliegende bzw. latente Merkmale geben, die zu den Einzelbewertungen auf den Dimensionen führen. Wenn dem so ist, dann müssten die Variablen, die zu einer latenten Dimension gehören, stark miteinander korrelieren, also z.B. alle Bewertungen zu den öffentlich-rechtlichen Sendern. Das bedeutet, wir könnten diese Bewertungen auch auf diese latente Dimension reduzieren. Zack fertig: Dimensionsreduktion!
Faktorenanalysen dienen genau dieser Dimensionsreduktion. Sie werden eingesetzt, um latente Konstrukte zu identifizieren, die die Ausprägungen der gemessenen (manifesten) Variablen bestimmen bzw. determinieren. Wir setzen die Faktorenanalyse aber auch ein, um z.B. das Problem der Multikolliniearität bei Regressionsanalysen in den Griff zu kriegen. Das Ziel der Faktorenanalyse ist es daher, eien Vielzahl an Variablen auf wenige zugrundeliegende Faktoren zu reduzieren, die man dann gut und klar unterscheiden kann, die also nicht oder nur wenig miteinander korrelieren.
8.6 Ablauf einer Faktorenanalyse
Die Schritte der Faktorenanalyse sind: 1. Voranalyse über Korrelationstabellen (Ausschluss von Variablen, die mit keiner anderen korrelieren) 2. Extraktion der Faktoren 3. Bestimmung der Anzahl Faktoren (Scree Plots der Eigenwerte) 4. Rotation der Faktoren (bessere Verteilung der Varianzaufklärung) 5. Eignung für die Variablen über Kommunalitäten (ggf. Ausschluss gering abgebildeter Variablen mit Kommunalitäten < .4) 6. Interpretation und Benennung der Faktoren 7. Speichern der Faktoren (eine Form der Indizes) für weitere Analysen, wie Regressionen
Um eine explorative Faktorenanalyse durchzuführen, müssen:
- die Variablen ausgewählt werden, für die Faktoren extrahiert werden sollen.
- Vor der eigentlichen Faktorenanalyse sollte man sich die Korrelationsmatrix anschauen, sowie
- KMO und Bartlett-Test durchlaufen lassen, um festzustellen, ob die Variablen für eine FA geeignet sind.
- Im Anschluss wählt man die Extraktionsmethode und Rotationsmethode aus, die verwendet werden sollen.
- Dann werden die Faktoren extrahiert und man erstellt einen Screeplot, um anhand des Ellenbogenkriteriums zu bestimmen, wie viele Faktoren eine gute Faktorlösung ergeben. Ein gängiges Kriterium ist, dass die Eigenwerte der Faktoren grösser als 1 sein sollen.
- Die Rotation der Faktoren ist nicht zwingend, aber so üblich und empfohlen, dass sie im Grunde dazugehört. Die Faktoren können orthogonal (unkorreliert) oder oblique (leichte Korrelationen zugelassen) rotiert werden.
- Dann kann man anhand der Faktorladungen jeder Variable die Faktoren interpretieren. Dabei charakterisieren hohe Faktorladungen die Faktoren.
- Haben einzelne Variablen kleine Kommunalitäten (<.4) und eine hohe Uniqueness (>.6), entfernt man diese Variablen aus den Faktorenanalysemodell und führt die Schritte 2. bis 6. nochmals ohne die entsprechenden Variablen aus.
- Am Ende schaut man sich noch die Gesamtvarianzaufklärung der Faktorlösung an, die Auskunft darüber gibt, wie gut die Faktoren insgesamt die Variablen abbilden.
- Will man die Faktorwerte als Indizes weiterverarbeiten, speichert man sie am Ende in seinen Daten ab.
8.7 The R Anxiety
Als Beispiel wird hierfür die R-Angst-Skala von Andy Field [@Field2022] verwendet. Die Fragen und zugehörigen Variablen sind:
- raq_01: Statistics make me cry
- raq_02: My friends will think I’m stupid for not being able to cope with R
- raq_03: Standard deviations excite me
- raq_04: I dream that Pearson is attacking me with correlation coefficients
- raq_05: I don’t understand statistics
- raq_06: I have little experience of computers
- raq_07: All computers hate me
- raq_08: I have never been good at mathematics
- raq_09: My friends are better at statistics than me
- raq_10: Computers are useful only for playing games
- raq_11: I did badly at mathematics at school
- raq_12: People try to tell you that R makes statistics easier to understand but it doesn’t
- raq_13: I worry that I will cause irreparable damage because of my incompetence with computers
- raq_14: Computers have minds of their own and deliberately go wrong whenever I use them
- raq_15: Computers are out to get me
- raq_16: I weep openly at the mention of central tendency
- raq_17: I slip into a coma whenever I see an equation
- raq_18: R always crashes when I try to use it
- raq_19: Everybody looks at me when I use R
- raq_20: I can’t sleep for thoughts of eigenvectors
- raq_21: I wake up under my duvet thinking that I am trapped under a normal distribution
- raq_22: My friends are better at R than I am
- raq_23: If I am good at statistics people will think I am a nerd
8.8 Korrelationsmatrix
Die Korrelationsmatrix ist die Basis für Faktorenanalysen (im Grunde braucht man nur die Korrelationsmatrix (+ Fallzahl) und die ursprünglichen Daten nicht). Mit dem folgenden Befehlen kann man sich die Korrelationsmatrix rausgeben lassen.
Man sieht hier schon, dass die Korrelationen nicht wahnsinnig gross sind, aber sich wie Haufen bilden. Die Variablengruppen, die untereinander hoch korrelieren, gehen vermutlich auf ein gemeinsames latentes Konstrukt zurück. Diese latenten Konstrukte werden im Folgenden auch als Faktoren bezeichnet.
8.9 Anzahl Faktoren bestimmen
Wenn wir eine Faktorenlösung suchen, müssen wir erstmal die Anzahl sinnvoller Faktoren bestimmen. Das geht mit dem «psych»-Paket und der Analyse «fa.parallel». Dort werden die Eigenwerte (eigen values) der Faktoren angezeigt. Die Eigenwerte sind der Anteil der Varianzaufklärung eines Faktors relativ zur Anzahl der Variablen in der Faktorenanalyse. Wenn also ein Eigenvalue bei 1 ist, erklärt ein Faktor so viel wie eine einzelne Variable.
## Parallel analysis suggests that the number of factors = 4 and the number of components = NA
Im Plot kann man sehen, dass der erste Faktor einen Eigenwert von knapp 6 hat, also so viel Varianz aufnimmt, wie sechs Variablen im Ursprung. Der zweite Faktor ist noch über 1. Das bedeutet, er erklärt etwas mehr als eine Ursprungsvariable. Der dritte und der vierte Faktor erklären etwas weniger als eine Variable. Da das alte Kaiser-Kriterium (Eigenwerte müssen über 1 sein) etwas sehr holzschnittartig ist, haben sich findige Statistiker ausgedacht, dass man die FA simulieren könnte, unter der Annahme, dass die Faktoren nichts erklären. Diese Simulation durch mehrfaches ziehen von Stichproben aus den Daten (FA Resampled Data) ergibt, dass 4 Faktoren mehr besser sind als die informationslose Simulation. Also ist die Faktorlösung 4.
8.10 Faktorladungen und Uniqueness
Mit dieser Analyse können wir jetzt die Faktorenanalyse rechnen. Als Anzahl «n» der Faktoren (nfactors) geben wir die 4 aus der Analyse von oben ein (siehe Abbildung @ref(fig-Parallelanalyse)).
Die Uniqueness ist der Varianzanteil, den eine Variable ganz alleine hat, also nicht mit den anderen teilt. Die Uniqueness ist das Gegenteil von Kommunalität (Communality), also der gemeinsamen Varianz mit der Faktorenlösung. Rechnerisch ergibt sich je Variable die Uniqueness aus 1 - Kommunalität. Hohe Uniqueness bedeutet, dass eine Variable nicht gut in die Faktorenanalyse passt, weil sie eben nicht gut durch die Faktoren abgebildet wird, sondern einzigartig (unique) ist. Für die Variable selbst und ggf. für ihre Integration in ein Modell ist eine hohe Uniqueness gut, da sie auch bedeutet, dass es keine Probleme mit Multikollinearität gibt. Die Variable kann also getrost aus der Faktorenanalyse entfernt und als eigenständige Variable in eine Modell aufgenommen werden.
Die Complexity gibt an, wie viele Faktoren gebraucht werden, um die Variable abzubilden. Wenn sie 1 ist, dann wird eine Variable von einem Faktor abgebildet. Ist sie zum Beispiel 1.97 braucht es zwei Faktoren, um die Variable darzustellen. Geringe Komplexität ist in dem Fall gut, da sie zu einer klaren Faktorenlösung führt.
## Loading required namespace: GPArotation
## # Rotated loadings from Factor Analysis (oblimin-rotation)
##
## Variable | MR1 | MR2 | MR4 | MR3 | Complexity | Uniqueness
## ---------------------------------------------------------------
## raq_06 | 0.84 | | | | 1.00 | 0.27
## raq_18 | 0.63 | | | | 1.03 | 0.57
## raq_13 | 0.57 | | | | 1.02 | 0.66
## raq_07 | 0.56 | | | | 1.02 | 0.65
## raq_10 | 0.49 | | | | 1.08 | 0.74
## raq_15 | 0.48 | | | | 1.04 | 0.71
## raq_05 | 0.45 | | 0.39 | | 1.97 | 0.46
## raq_14 | 0.42 | | | | 1.06 | 0.78
## raq_09 | | 0.81 | | | 1.02 | 0.38
## raq_23 | | 0.79 | | | 1.02 | 0.41
## raq_19 | 0.26 | 0.56 | | | 1.41 | 0.50
## raq_22 | | 0.52 | | | 1.29 | 0.59
## raq_02 | 0.25 | 0.48 | | | 1.54 | 0.62
## raq_21 | | | 0.59 | | 1.04 | 0.60
## raq_04 | | | 0.56 | | 1.01 | 0.67
## raq_20 | | | 0.54 | | 1.02 | 0.68
## raq_16 | | | 0.51 | | 1.02 | 0.77
## raq_03 | | | -0.43 | | 1.01 | 0.80
## raq_01 | | | 0.39 | | 1.06 | 0.83
## raq_12 | | | 0.37 | | 1.07 | 0.89
## raq_08 | | | | 0.88 | 1.00 | 0.25
## raq_11 | | | | 0.72 | 1.00 | 0.45
## raq_17 | | | | 0.68 | 1.00 | 0.51
##
## The 4 latent factors (oblimin rotation) accounted for 40.12% of the total variance of the original data (MR1 = 13.23%, MR2 = 9.76%, MR4 = 8.83%, MR3 = 8.29%).
8.11 Interpretation der Faktorenanalyse
Mit dieser Faktorenlösung können wir jetzt die Faktoren interpretieren.
Der erste Faktor lädt hoch auf folgenden Items. Wir können diesen Faktor als l «Probleme mit Computern»** labeln:
- raq_05: I don’t understand statistics
- raq_06: I have little experience of computers
- raq_07: All computers hate me
- raq_10: Computers are useful only for playing games
- raq_13: I worry that I will cause irreparable damage because of my incompetence with computers
- raq_14: Computers have minds of their own and deliberately go wrong whenever I use them
- raq_15: Computers are out to get me
- raq_18: R always crashes when I try to use it
Beachte: Das Item «raq_05» lädt auch hoch auf dem zweiten Faktor MR2.
Wenn man die Fragen anschaut, die hoch auf dem zweiten Faktor MR2 laden, deuten darauf hin, dass es die Befragten Angst haben, von ihren Peers komisch angesehen zu werden. Nennen wir diesen Faktor «Angst vor sozialer Bewertung»:
- raq_02: My friends will think I’m stupid for not being able to cope with R
- raq_09: My friends are better at statistics than me
- raq_19: Everybody looks at me when I use R
- raq_22: My friends are better at R than I am
- raq_23: If I am good at statistics people will think I am a nerd
Beim Faktor MR3 wird deutlich, dass es hier eine Angst vor Statistik gibt. Nennen wir den Faktor «Angst vor Stastik»:
- raq_01: Statistics make me cry
- raq_03: Standard deviations excite me
- raq_04: I dream that Pearson is attacking me with correlation coefficients
- raq_05: I don’t understand statistics
- raq_12: People try to tell you that R makes statistics easier to understand but it doesn’t
- raq_16: I weep openly at the mention of central tendency
- raq_20: I can’t sleep for thoughts of eigenvectors
- raq_21: I wake up under my duvet thinking that I am trapped under a normal distribution
Bei den übrigen Fragen, die auf dem Faktor MR4 laden, geht es eher um Mathematik. Wir könnten also sagen der Faktor MR4 ist «Angst vor Mathe».
- raq_08: I have never been good at mathematics
- raq_11: I did badly at mathematics at school
- raq_17: I slip into a coma whenever I see an equation
Wir können noch schauen, ob die Variablen mit Doppelladungen plausibel sind. Also schauen wir zum Beispiel auf das Item raq_05 «I don’t understand statistics». Das scheint mit einer geringen Selbstwirksamkeit in Bezug auf Computer und Statistik zusammenzuhängen. Es spiegelt die Selbsteinschätzung wieder, dass man Statistik und Computer «nicht kann».
8.12 Faktorendiagramm
Faktorenanalysen kann man mit solchen Diagrammen darstellen. Hier sieht man auch, wie stark die einzelnen Faktoren miteinander korrelieren, wenn man die Faktoren nicht gezwungen hat, orthogonal zu sein, also unkorrliert.
## Saving 6 x 4 in image
LEF 8
Essayfragen 8
E8.1 Bitte schauen Sie sich die folgenden drei (mistigen) Tabellen an, die aus fiktiven Publikationen stammen. Jedes mal werden die Koeffizienten als «beta» bezeichnet, aber es könnten nur die unstandardisierten Regressionskoeffizienten b sein oder die standardisierten Regressionskoeffizienten std. b. Was ist wohl was? a) Begründen Sie Ihre Entscheidung! b) Eklären Sie, warum es nicht \(\beta\)s sein können.
Essayfragen 8.2-8.3
E8.2 Schreiben Sie zu jedem der folgenden Schlagworte in ganzen Sätzen auf, zu welchen Zwecken Faktorenanalysen alles eingesetzt werden können: a) latente Variablen entdecken b) Multikollinearität c) Indices d) Messung latenter Konstrukte
E8.3 a) Erklären Sie, was Faktorladungen sind. b) Wie gehen Sie vor, wenn Sie Faktorladungen interpretieren wollen?
Essayfragen 8.4-8.10
Lesen Sie @song2004 und beantworten Sie dann folgende Fragen und die anschliessenden MCs.
E8.4 Was wird mit Cronbachs \(\alpha\) beschrieben?
E8.5 Was für eine Faktorenanalyse haben Song et al durchgeführt?
E8.6 Mit welcher Methoden wurden die Faktoren rotiert?
E8.7 Wie viele initiale Faktoren hat die Faktorenanalyse herausgegeben?
E8.8 a) Nach welchen Kriterien wurden die Faktoren ausgewählt? b) Wie bewerten Sie das Vorgehen?
E8.9 a) Wie viele Items (Fragen) wurden mit der Faktorenanalyse ausgewertet? b) Auf wie viele Faktoren wurden die Items reduziert?
E8.10 Wie viel Varianz erklären die Faktoren, die von Song et al ausgewählt wurden?
MC-Fragen 8
MC 8.1.
MC 8.2.
MC 8.3.
MC 8.4.
MC 8.5.
MC 8.6.
Insgesamt (nur MCs) von Punkten, was % und etwa einer entspricht.