Tweede beweging · Verschillen koppelen

Oefening 3.1 · Correlatie en scatterplot

Samenhang tussen twee variabelen · correlation & scatterplot

Data bij deze oefening — diamantjes.csv de twaalf steentjes uit het boek · penguins.csv 342 pinguïns · bigfive.csv 2634 ingevulde vragenlijsten · wat elke kolom betekent

prikkel

Langere snavel, ondiepere snavel?

samenhang · negatieve r

Ergens op een kale rots bij Antarctica staat iemand met een schuifmaat en 342 pinguïns die daar niet om gevraagd hebben.1 Van elke pinguïn schrijft hij twee dingen op: hoe lang de snavel is, en hoe diep. Met diep bedoelen we de hóógte van de snavel — van boven naar onder, van opzij bekeken, niet de breedte en niet de lengte. Meer niet: twee getallen per beest, netjes onder elkaar.

De snavel in profiel: lengte loopt van voor naar achter, diepte is de hoogte van boven naar onder.

En nu wij. Want in die twee kolommetjes zit iets raars verstopt. Reken je de samenhang tussen snavellengte en snaveldiepte uit over alle 342 pinguïns, dan komt er een getal uit dat negatief is: r = −.24. In gewone taal: hoe langer de snavel, hoe minder diep — alsof een dier met een lange neus per se een dunne neus heeft.

Voelt dat kloppen? Mij niet. Trek even mee — we gaan het zelf uitzoeken, en onderweg blijkt dat getal iets heel anders te zeggen dan het lijkt.

speel het

Waar r vandaan komt — de twaalf steentjes eerst

covariantie · sum of products

Twee kolommen, twaalf steentjes: karaat en glans, netjes naast elkaar. En jij wilt één getal dat zegt hoe die twee samen bewegen. Het halve werk heb je in het boek al gedaan — van elk steentje ken je de afwijking van het gemiddelde, de blauwe streepjes. Alleen heeft elk steentje er hier twéé: eentje op karaat (het gemiddelde is 1.60) en eentje op glans (gemiddelde 50.00).

Dan de zet die je uit het boek kent: vermenigvuldig die twee afwijkingen per steentje met elkaar. Wijzen ze dezelfde kant op — allebei erboven, of allebei eronder — dan komt er plus uit; spreken ze elkaar tegen, dan min.

Doe er twee met de hand, dan zie je het gebeuren. De twaalf steentjes hebben in het boek namen, en twee ervan komen hier goed uit. Anna weegt 1.0 karaat en glanst 20, dus haar afwijkingen zijn −0.6 karaat en −30 glanspunten: allebei onder het gemiddelde, en min keer min is plus. Haar bijdrage is +18. Kees weegt 2.1 karaat en glanst 30 — zwaar, maar dof. Zijn afwijkingen zijn +0.5 en −20, die spreken elkaar tegen, en zijn bijdrage is −10.

afwijking op karaat afwijking op glans die twee maal elkaar
Anna −0.6 −30 +18
Kees +0.5 −20 −10

Bij onze twaalf duwen acht steentjes de som omhoog, drie liggen precies op een gemiddelde en tellen dus nul, en alleen Kees telt tegen. Tel alles op en er staat 54.00 — die som heet de sum of products. Deel door n − 1 = 11, en daar is hij: de covariantie (Sxy), 4.91. Het gemiddelde rechthoekje.

En nu de vraag die je bij élk getal hoort te stellen: 4.91 wát? Reken de eenheid maar gewoon mee — een afwijking in karaat keer een afwijking in glans is… karaat-glanzen. Is 4.91 karaat-glanzen veel samenhang? Weet je niet. Niemand denkt in karaat-glanzen. Ruw is ruk. De covariantie zit nog op de schaal van zijn twee variabelen; je kunt er de ríchting aan aflezen — plus, dus zwaarder gaat samen met glanzender — maar niet de sterkte.

De reparatie is één deling. Deel de covariantie door de twee standaardafwijkingen — 0.3908 karaat en 21.7423 glanspunten — en kijk wat er met de eenheden gebeurt: karaat streept weg tegen karaat, glans tegen glans. Wat overblijft is een kaal getal zonder eenheid, de gestandaardiseerde versie van dezelfde samenhang, en dat getal kan nooit buiten −1 en +1 komen: r = .58. Dát kun je in één keer lezen, bij diamantjes net zo goed als bij pinguïns.

OpmerkingEerst zoals in het boek — bouw de correlatie zelf · JASP 0.98.1

Open diamantjes.csv en maak twee kolommen: klik op de + rechts van de laatste kolom, geef een naam, kies type Scale en het icoon R.

S_xy     ->   cov(karaat, glans)
r_zelf   ->   cov(karaat, glans) / (sd(karaat) * sd(glans))

In S_xy staat in elke rij hetzelfde getal: de covariantie, 4.909 — dat is de 4.91 uit de tekst hierboven, met de derde decimaal erbij die JASP standaard toont. In r_zelf staat 0.578; daar is de eenheid uit weggedeeld.

Dan hetzelfde met één klik: Regression → Correlation, karaat en glans bij Variables, Pearson’s r aan (die staat al aan), en onder Plots de Scatter plots. Bij Pearson’s r staat 0.578.

En welke lijn daarbij hoort: Regression → Linear Regression, glans als Dependent Variable en karaat bij Covariates. In Coefficients staat (Intercept) −1.429 en karaat 32.143.

Correlation is dus geen zwarte doos: het is jouw covariantie, gedeeld door jouw twee standaardafwijkingen. JASP doet alleen het strepen sneller.

Let op het verschil tussen wat je scherm zegt en wat je zelf opschrijft. JASP drukt 0.578 af, met nul en al; in de tekst hierboven staat r = .58. Dat is geen slordigheid van een van de twee. Een klikpad citeert je scherm; jouw eigen zin volgt APA, en daar laat je de nul weg omdat het ontbréken ervan iets zegt: dit getal kan niet voorbij de 1 komen. De getallen in de JASP-kaders van dit boekje staan daarom zoals JASP ze toont, zodat je kunt nakijken of jij hetzelfde voor je hebt. Kwadrateer die .58 en je hebt r² = .33 — een derde van de verschillen in glans loopt langs het gewicht, de rest niet.

OpmerkingEerst de twaalf diamantjes

Twee vragen voor je verder leest, en dan mag je los op iets nieuws.

a) Een karaat is gewoon een gewicht: 200 milligram. Weeg de steentjes dus eens in milligram (maak een kolom mg met karaat * 200) en vraag de covariantie opnieuw op. (één getal, en één zin over wat er wél en niet veranderd is aan de twaalf steentjes)

b) Deel die nieuwe covariantie door de twee standaardafwijkingen die er nu bij horen. (één getal, en zeg wat je herkent)

a) 981.82 milligram-glanzen — tweehonderd keer de covariantie van daarnet. (Reken je 200 × 4.91 met de hand, dan krijg je 982.00; JASP vermenigvuldigt de onafgeronde 4.9091. Zelfde soort scheefte als overal waar je met een gedrukt getal verder rekent.) En er is niets aan de steentjes veranderd: geen enkel gewicht, geen enkele glans, geen enkel punt in de wolk. Alleen de maatstok waarmee je het gewicht opschrijft is anders. Een getal dat tweehonderd keer zo groot wordt zodra je een andere maat kiest, kan niet zeggen hoe sterk een verband is.

OpmerkingIn JASP · de steentjes in milligram
mg        ->   karaat * 200
S_xy_mg   ->   cov(mg, glans)

Er komt 981.818 uit.

b) .58 — hetzelfde getal als daarnet. Delen door de standaardafwijkingen haalt de maatstok eruit, en dat is precies waarom de correlatie bestaat.

OpmerkingIn JASP · en weer gestandaardiseerd
r_mg      ->   cov(mg, glans) / (sd(mg) * sd(glans))

Er staat 0.578, precies wat r_zelf ook gaf.

Nu jij — bij de pinguïns, die je inmiddels kent. Nieuwe vraag aan bekende beesten: hangt de lengte van een snavel samen met zijn diepte?

Haal de pinguïns binnen en gooi ze in een plaatje. Elke pinguïn één stip: snavellengte op de ene as, snaveldiepte op de andere.

OpmerkingIn JASP · trek eraan · JASP 0.98.1
  1. Open penguins.csv (menu linksboven → Open).
  2. Regression → Correlation, en sleep snavellengte en snaveldiepte naar Variables. Laat Pearson’s r aan staan.
  3. Vink onder Plots de Scatter plots aan: 342 stippen, een wolk.

Bij Pearson’s r staat −0.235.

Daar is-ie: een wolk die van linksboven naar rechtsonder hangt, en het getal −0.235 is één cijfer voor “hoe strak kantelt die wolk omláág?”. De schaal loopt van −1 (een perfecte dalende lijn) via 0 (geen rechtlijnig verband) naar +1 (een perfecte stijgende lijn). −.24 is een flauwe daling met veel rommel eromheen — zwak, maar het staat er. Langere snavel, ondiepere snavel. Zegt het getal.

snap het

Drie soorten in één wolk

confounding · derde speler · omgeklapt teken · r²

De verleiding is nu te schrijven: “bij pinguïns gaat een langere snavel samen met een ondiepere.” Het getal is netjes uitgerekend, de zin lijkt te volgen. En toch is-ie mis. Kleur de wolk eens per soort:

OpmerkingIn JASP · kleur de wolk · JASP 0.98.1
  1. Zet in dezelfde Correlation-analyse, onder Plots, bij de scatterplot de optie Color by / group op soort. JASP kleurt de stippen dan per soort.
  2. En reken r nu apart uit binnen één soort: klik op het trechtertje boven de data (Filter) en typ in het R-vak soort == "Gentoo". De Correlation rekent opnieuw, en bij Pearson’s r staat 0.643.
  3. Zet het filter daarna weer uit — anders rekent alles hierna over de Gentoo’s.

Ineens zie je het. Er zijn drie soorten in die wolk, elk als een eigen kluitje. En binnen elk kluitje loopt de samenhang de andere kant op: langere snavel, juist diepere snavel. Reken maar na — Adélie +.39, Chinstrap +.65, Gentoo +.64. Alle drie positief. Het teken is omgeklapt.

Wat gebeurde er? De Gentoo’s (rechtsonder) hebben lange, ondiepe snavels; de andere twee soorten (linksboven) korte, diepe. Gooi je ze op één hoop, dan trekt dat soortverschil de hele wolk schuin naar beneden — en dáár komt die −.24 vandaan. Niet uit een verband tussen snavellengte en -diepte, maar uit het verschil tússen de soorten. De soort is een derde speler die zich als het verband vermomt. Het boek noemt dat straks confounding; jij zag het hier gebeuren, in de kleuren van een scatterplot.

En dat is de dubbele les van dit blok, mooier dan één van de twee alleen:

  • De rekenmachine zei: over deze 342 dieren is r precies −.24. Daar valt niets op af te dingen; reken het maar na.
  • De scatterplot zei: klopt — maar niet wát je denkt. De richting is een luchtspiegeling van de soort.

Uitgerekend en waar zijn dus twee verschillende dingen. Een correlatie kan tot achter de komma kloppen en je tóch de verkeerde kant op wijzen. Daarom lees je r nooit zonder de wolk ernaast — en zeker niet zonder je af te vragen wie er nog meer in die wolk zit. Zou je het in één eerlijke zin moeten opschrijven, dan niet “langere snavels zijn ondieper”, maar: over drie soorten heen hangt de wolk omlaag, maar binnen elke soort omhoog — het min-teken is een soorteffect.

(En dat kleine getal r² dat je weleens ziet — hier .06 — is het “aandeel gedeelde variantie”: van alle verschil in snaveldiepte deelt maar een paar procent iets met snavellengte. Ook dat is geen “verklaring”, alleen boekhouding.)

jouw beurt

Van snavel naar vin

correlatie berekenen · kleuren per soort · soortverschil

Van snavel naar vin. In penguins.csv zitten ook de kolommen vinlengte (in mm) en gewicht (in gram). Zelfde recept, andere vraag: hangen die twee samen?

OpmerkingIn JASP · zet klaar · JASP 0.98.1

Open penguins.csv. De covarianties maak je in het R-vak, zoals bij de diamantjes; de correlatie en de wolk komen uit Regression → Correlation, en de gemiddelden per soort uit Descriptives.

OpmerkingJouw beurt

Beantwoord de vragen op volgorde; ze bouwen op elkaar. De eerste drie zijn de route die je hierboven op de twaalf steentjes hebt zien lopen — ruw, en dan pas gestandaardiseerd. De laatste heet met opzet Z) en niet k): de wolk-vraag staat altijd achteraan.

Het zijn er elf, en dat is te veel voor één zitting. Na f) heb je een hele correlatie van begin tot eind gemaakt — dat is een goed moment om te stoppen en er morgen weer in te vliegen. Bij g) begint iets anders: daar gaat het niet meer over het getal maar over wie er in de wolk zit.

Tussen haakjes staat telkens wat er in een compleet antwoord hoort. Staat er meer antwoorden zijn goed, of dat een deel aan jou is, dan is de verwoording vrij — bij de rest ligt het antwoord vast. En onder Tentamenvragen, aan het eind van deze oefening, komen dezelfde onderwerpen nog een keer langs in de vorm waarin een tentamen ze stelt.

a) Reken de covariantie uit tussen de vinlengte en het gewicht. (rapporteer Sxy mét zijn eenheid — en die eenheid is twee woorden lang)

b) Doe hetzelfde nog eens, maar zet het gewicht eerst in kilo’s (maak een kolom kg met gewicht / 1000). (het nieuwe getal, en één zin over waarom dit een probleem is voor een maat die sterkte zou moeten uitdrukken)

c) Deel allebei je covarianties door de twee standaardafwijkingen die erbij horen. (twee getallen, en zeg wat er met de eenheid gebeurd is)

d) Teken de wolk. (één zin over welke kant hij op hangt en hoe strak — de richting ligt vast, hoe strak is jouw oordeel)

e) Zeg in gewone taal wat het teken hier betekent. (één zin, over pinguïns — de richting ligt vast, de verwoording is aan jou; een getal mag erbij als je zegt wat het meet)

f) Kwadrateer de correlatie. (rapporteer r², en één zin over wat dat getal wél zegt en wat het niet zegt)

g) Vertrouw dat getal niet voor je de soorten hebt gezien. Kleur de wolk per soort en reken r apart binnen elke soort. (drie correlaties, elk met zijn n)

h) En toch is de r over alle 342 heen hóger dan binnen élke soort afzonderlijk. Hoe kan dat?

Doe het in twee stappen. Vraag eerst per soort het gemiddelde gewicht en de gemiddelde vinlengte op — zes getallen, één Descriptives met soort bij Split. Kijk daarna naar die zes en probeer te zeggen wat de Gentoo’s anders maakt. (twee of drie zinnen; kom je er niet uit, dan is het antwoord hieronder geen nederlaag maar de bedoeling)

i) Bij de snavels hierboven klapte het teken om; hier niet. Wat is het verschil? (twee zinnen)

j) Schrijf je bevinding op in één eerlijke zin. (N, allebei de variabelen, de r over alles én het bereik binnen de soorten — die vier moeten erin staan, de woorden eromheen zijn aan jou; meer antwoorden zijn goed)

Z) Is de wolk recht, en zitten er geen punten tussen die de boel verneuken? Een r vat één ding samen: hoe strak de punten om een rechte lijn liggen. Kijk dus terug naar je wolk uit d) — buigt hij ergens, en hangt er één punt zo ver weg dat het in zijn eentje aan het getal trekt? (twee korte uitspraken, en zeg welke van je antwoorden hierboven eraan hangen)

a) Sxy = 9824.42 millimeter-grammen. Zes cijfers voor de komma, en een eenheid die je nergens anders tegenkomt.

OpmerkingIn JASP
S_vg      ->   cov(vinlengte, gewicht)

Er komt 9824.416 uit.

b) In kilo’s is het 9.8244 millimeter-kilogrammen. Precies duizend keer kleiner, want je hebt duizend keer kleinere getallen in de ene kolom gezet — en verder niets. Geen enkele pinguïn is van gewicht veranderd, geen enkel punt is in de wolk verschoven, en toch staat er nu een ander getal. Een maat die van je weegschaal afhangt, kan niet zeggen hóé sterk een verband is.

OpmerkingIn JASP
kg        ->   gewicht / 1000
S_vkg     ->   cov(vinlengte, kg)

Er komt 9.824 uit — op de vier decimalen van het antwoord hierboven is dat 9.8244.

c) Allebei .87. De standaardafwijking van de vinlengte is 14.0617 mm, die van het gewicht 801.9545 g, en die van het gewicht in kilo’s 0.8020 kg. Deel de covariantie door het product van de twee die erbij horen en je krijgt twee keer hetzelfde kale getal.

En dat kómt doordat de eenheid weg is. In de teller stonden millimeter-grammen; in de noemer staat millimeter maal gram. Die vallen tegen elkaar weg, en wat overblijft kan niet meer weten in welke maat je gewogen hebt.

OpmerkingIn JASP
r_vg      ->   cov(vinlengte, gewicht) / (sd(vinlengte) * sd(gewicht))
r_vkg     ->   cov(vinlengte, kg) / (sd(vinlengte) * sd(kg))

Allebei 0.871. Dezelfde 0.871 staat bij Regression → Correlation onder Pearson’s r.

d) Een strak stijgende wolk, van linksonder naar rechtsboven, en veel minder rommel eromheen dan bij de snavels.

OpmerkingIn JASP

De wolk: Regression → Correlation met vinlengte en gewicht, en onder Plots de Scatter plots aan.

e) Pinguïns met langere vinnen zijn zwaarder. Meer niet — er staat niet dat de vin het gewicht maakt.

f) r² = .76. Drie kwart van alle verschil in gewicht loopt langs de vinlengte mee. Wat het níét is: een verklaring. Het is boekhouding — welk deel van de verschillen de twee maten delen, niet waaróm ze dat doen.

g) Adélie r = .47 (n = 151), Chinstrap r = .64 (n = 68), Gentoo r = .70 (n = 123).

OpmerkingIn JASP

Kleur de wolk met Color by / group op soort, en zet met het trechtertje (Filter) telkens één soort aan: soort == "Adelie", en zo verder. Bij Pearson’s r staat 0.468, 0.642 en 0.703. Zet het filter daarna weer uit.

h) Houd je gekleurde plaatje uit g) erbij — dáár staat het antwoord in, en de zes getallen hieronder zeggen alleen wat je erop ziet. Vraag ze op:

OpmerkingIn JASP

Descriptives, met gewicht en vinlengte bij Variables en soort bij Split. Er staat 3700.662 · 3733.088 · 5076.016 gram en 189.954 · 195.824 · 217.187 mm.

Kijk nu naar wat de Gentoo’s zijn: veruit de zwaarste dieren (5076 gram tegen 3701 en 3733) én de langste vinnen (217.2 mm tegen 190.0 en 195.8). Zwaar én lange vin — dat is dezelfde richting als het verband binnen elke soort. Het soortverschil duwt de wolk dus mee langs die diagonaal in plaats van er dwars op te staan.

En dan tel je over alle 342 heen twéé dingen bij elkaar op: het verband binnen elke soort, én de sprong van de ene soort naar de andere. Twee kluitjes die al schuin liggen én bovendien lángs diezelfde diagonaal van elkaar af liggen, vormen samen één lange strakke wolk — strakker dan elk kluitje afzonderlijk.

Wil je zien hoe groot die sprong is, zet hem dan naast de rommel bínnen een soort — maar dat is een verdieping bovenop je antwoord, niet iets wat je uit de zes getallen had moeten halen:

OpmerkingIn JASP
sprong    ->   sd(tapply(gewicht, soort, mean))
rommel    ->   sqrt(mean(tapply(gewicht, soort, var)))

Er komt 784.868 en 451.711 uit — de 785 en 452 uit het antwoord hierboven, op hele grammen afgerond.

De sprong tussen de soorten is dus ruim anderhalf keer zo groot als de rommel binnen één soort (785 gedeeld door 452 is 1.74), en dat hele stuk telt in die .87 volledig mee.

i) Bij de snavels lag het soortverschil dwars op het verband binnen de soorten: de Gentoo’s hebben lange én ondiepe snavels, terwijl binnen elke soort een langere snavel juist een diepere is. Daar duwde de soort de wolk de andere kant op, en klapte het teken om. Hier ligt het soortverschil mee, en dan versterkt het wat er binnen de soorten al gebeurt. Zelfde derde speler, twee tegengestelde uitkomsten — en dat is precies waarom je aan het teken alleen niet genoeg hebt.

j) Bijvoorbeeld:

Bij 342 pinguïns hingen de vinlengte en het gewicht sterk positief samen, r = .87; binnen elk van de drie soorten was die samenhang zwakker, r = .47 tot .70.

Let op dat die tweede helft er staat. Zonder haar leest je lezer .87 als het verband bij één soort pinguïn, en dat is het niet.

Z) Ja — allebei komen ze hier goed uit. Maar het loont om te weten wát eraan zou hangen als dat anders was.

De wolk is recht. Je ziet in het plaatje geen bocht, en je kunt het narekenen. Trek er een lijn doorheen en kijk hoeveel die dekt: dat is R², en met één voorspeller is dat exact hetzelfde getal als de r² uit f). Op vier decimalen .7590. Sta de lijn nu een knik toe — een extra term met de vinlengte in het kwadraat — en er staat .7753.

Ruim anderhalve procent erbij dus, en toch is dat géén bewijs van een bocht. Een extra term kan de R² namelijk alleen maar omhoog duwen en nooit omlaag — ook een term die niets te melden heeft. (Dat is geen bevinding maar rekenkunde, en je komt het bij de meervoudige regressie nog eens tegen, waar een kolom willekeurige getallen precies datzelfde doet.) Wat je hier ziet is dus de gewone winst van een parameter erbij, niet een kromming die om aandacht vraagt.

OpmerkingIn JASP

Maak een kolom vin2 met vinlengte^2 en draai Regression → Linear Regression twee keer: eerst met alleen vinlengte bij Covariates, daarna met vinlengte én vin2. In Model Summary staat R² = 0.759 en daarna 0.775.

En geen punt trekt in zijn eentje aan het getal. Laat de zwaarste pinguïn weg en de correlatie beweegt pas in de vierde decimaal:

OpmerkingIn JASP

Zet in het Filter gewicht < max(gewicht) — dan valt precies het zwaarste dier weg, en je hoeft zelf geen gewicht in te tikken. Bij Pearson’s r staat dan 0.872 in plaats van 0.871; op vier decimalen is dat .8716 tegen .8712.

Doe je dat om de beurt voor alle 342, dan is de grootste verschuiving die één enkel dier veroorzaakt nog altijd maar .004. Dat is de derde decimaal, en die haalt geen enkele conclusie onderuit.

OpmerkingIn JASP · en wat hier niet kan

Dat om de beurt voor alle 342 doen is geen knop in JASP. Het R-broertje rekent het in twee regels na, en daar komt .0042 uit.

Wat eraan hángt, als het wél mis was gegaan:

  • Blijft staan: de drie correlaties binnen de soorten bij g), de zes soortgemiddelden bij h), en de redenering bij i). Die beschrijven wat er in de tabel staat, en daar verandert geen kromme wolk iets aan.
  • Doe je over: de covariantie bij a) en b) niet — die is gewoon een som — maar wél alles wat de sterkte uitdrukt: de r bij c), de r² bij f) die er het kwadraat van is, de uitspraak in pinguïn-taal bij e) die erop leunt, en de zin bij j). Een r meet rechtlijnigheid; ligt het verband krom, dan is het getal geen leugen maar wel een slechte samenvatting.

jouw beurt

Van vin naar vragenlijst

schaalscore · totaalscore · dezelfde route, ander onderwerp

Pinguïns hebben één groot voordeel: je kunt ze wegen. Bij mensen meet je zelden iets wat op een weegschaal past, en dan bouw je je variabele zélf uit een vragenlijst.

In bigfive.csv staan 2634 mensen die vijfentwintig uitspraken over zichzelf hebben beoordeeld op een schaal van 1 tot 6. Je gebruikt er hier tien: vijf gaan over neuroticisme — hoe snel raak je van slag — en vijf over openheid, de smaak voor nieuwe dingen en nieuwe ideeën. Eén rij is één mens.

Tien losse antwoorden zijn nog geen variabele. Je maakt er twee getallen van door per persoon het gemiddelde van zijn vijf antwoorden te nemen; dan staat er voor iedereen één neuroticisme-score en één openheid-score, allebei nog steeds op die schaal van 1 tot 6.

Zo’n getal heet een schaalscore: het gemiddelde van de vragen die bij één schaal horen. Het boek telt de vragen liever op en noemt dat een totaalscore (ook somscore, in het Engels sum score). Het is hetzelfde idee, alleen met grotere getallen, en de correlatie merkt daar niets van — bij b) en c) reken je zelf na waarom. Allebei zijn goed, als je er maar bij zegt welke je nam. Dit oefenboek middelt, omdat een gemiddelde op dezelfde 1 tot 6 blijft staan als de antwoorden zelf; zo lees je een score meteen terug als een antwoord op de vragenlijst.

En dan loop je precies hetzelfde rondje als hierboven. Dat is de bedoeling: de vinlengte is nu de openheid, en het gewicht is het neuroticisme, maar de route is regel voor regel dezelfde.

OpmerkingIn JASP · bouw de twee schalen · JASP 0.98.1

Open bigfive.csv en maak twee kolommen in het R-vak (type Scale):

neuroticisme   ->   rowMeans(data.frame(N1, N2, N3, N4, N5))
openheid       ->   rowMeans(data.frame(O1, O2, O3, O4, O5))

Let op wélk gemiddelde je hier neemt: één dwars over vijf kolommen, per persoon. Het mean(y)-blokje uit de sleep-route middelt juist één kolom over alle rijen — zie de gereedschapskist.

Kijk of je twee kolommen erbij hebt en of ze binnen de schaal blijven: zet openheid in Descriptives en kijk naar Minimum en Maximum. Er hoort iets tussen 1 en 6 te staan. Staat er een 0 of een 7, dan heb je een verkeerde kolom te pakken.

OpmerkingDezelfde route op vragenlijsten

Vier stappen, en dan de zin waar dit hoofdstuk op uitloopt.

a) Vraag van allebei de schalen het gemiddelde en de standaardafwijking op. (vier getallen, met de schaal erbij waarop ze staan)

b) Reken de covariantie uit tussen openheid en neuroticisme. Doe het daarna nog eens met de rijsóm in plaats van het rijgemiddelde (een kolom met N1 + N2 + N3 + N4 + N5, en net zo voor de O’s). (twee getallen, en zeg hoeveel keer ze schelen en waarom)

c) Deel allebei door de bijbehorende standaardafwijkingen. (twee getallen, en één zin over wat je hier herkent van de pinguïns)

d) Teken de wolk, en kwadrateer de correlatie. (rapporteer r², en één zin over hoe zo’n wolk eruitziet vergeleken met die van vin en gewicht)

e) Schrijf op wat je gevonden hebt, in twee zinnen en in deze vaste vorm:

Binnen deze n mensen geldt: … (r = …). Of dat ook buiten deze n opgaat, weet je nog niet.

(vul de n, de bevinding in gewone taal en de r in — en laat die tweede zin staan zoals hij er staat; voor de bevinding zijn meer antwoorden goed)

a) Neuroticisme: M = 3.16, SD = 1.19 schaalpunten. Openheid: M = 3.87, SD = 0.55 schaalpunten. Allebei op de schaal van 1 tot 6 waarop de vragen beantwoord zijn — dus de gemiddelde deelnemer zit bij neuroticisme net onder het midden en bij openheid er duidelijk boven.

OpmerkingIn JASP

Descriptives, met neuroticisme en openheid bij Variables: Mean 3.163 en 3.867, Std. deviation 1.195 en 0.546.

b) Met de rijgemiddelden is de covariantie 0.1050, met de rijsommen 2.6253. Dat scheelt precies een factor 25, en niet toevallig: een som is vijf keer een gemiddelde, en in een covariantie zitten twee variabelen, dus vijf maal vijf. Weer hetzelfde als bij de kilo’s van de pinguïns — je hebt de maat veranderd en niet de mensen.

OpmerkingIn JASP
N_som     ->   N1 + N2 + N3 + N4 + N5
O_som     ->   O1 + O2 + O3 + O4 + O5
cov_gem   ->   cov(neuroticisme, openheid)
cov_som   ->   cov(N_som, O_som)

Er komt 0.105 en 2.625 uit.

c) Allebei .16. Dat is wat je bij de pinguïns ook zag: zodra je door de twee standaardafwijkingen deelt, vallen de eenheden weg en maakt het niet meer uit of je in gemiddelden of in sommen rekende.

OpmerkingIn JASP

De standaardafwijkingen van de sommen staan in Descriptives: 5.973 voor N_som en 2.728 voor O_som.

r_gem     ->   cov(neuroticisme, openheid) / (sd(neuroticisme) * sd(openheid))
r_som     ->   cov(N_som, O_som) / (sd(N_som) * sd(O_som))

Allebei 0.161.

d) r² = .03. En de wolk zíet er ook zo uit: een dikke, ronde vlek zonder duidelijke richting, waar je met moeite een flauwe helling in ontdekt. Zet dat naast de vin-en-gewicht-wolk van hiervoor — daar lagen de punten in een strakke sliert. Hier delen de twee maten nog geen drie procent van hun verschillen.

OpmerkingIn JASP

Regression → Correlation met openheid en neuroticisme, Scatter plots aan. Bij Pearson’s r staat 0.161; kwadrateer dat zelf tot .03.

Dat de punten in blokjes lijken te vallen is geen fout: een gemiddelde van vijf hele getallen kan maar een beperkt aantal waarden aannemen, dus veel mensen landen op precies hetzelfde plekje.

e) Bijvoorbeeld:

Binnen deze 2634 mensen geldt: wie hoger scoort op openheid, scoort ook iets hoger op neuroticisme (r = .16). Of dat ook buiten deze 2634 opgaat, weet je nog niet.

Die tweede zin is geen slag om de arm en geen bescheidenheid. Alles wat je in dit hoofdstuk hebt uitgerekend — de covariantie, de correlatie, de r² — beschrijft de mensen die je vóór je hebt, en verder niemand. Wat er nodig is om die stap naar buiten wél te zetten, is precies waar de volgende oefening mee begint.

tentamen

Tentamenvragen

covariantie tegenover correlatie · derde speler · r²

Deze zes hebben de vorm waarin een tentamen dit onderwerp stelt: vier opties, en precies één is de beste — de andere drie kunnen best waar zijn en toch niet het beste antwoord zijn. Rekenen hoeft alleen bij vraag 2, en daar nog grof ook.

OpmerkingTentamenvragen

1. Een onderzoeker meldt over deze 342 pinguïns dat de covariantie tussen de vinlengte en het gewicht 9824.42 is, en schrijft erbij: “dat is een sterk verband.” Welke kritiek op die zin is de beste?

  1. Het getal klopt, maar hij had eerst moeten uitzoeken of het verband ook buiten deze 342 dieren geldt.
  2. Onjuist, want 9824.42 ligt buiten het bereik van −1 tot +1.
  3. De covariantie zegt de richting wél en de sterkte niet: zet je het gewicht in kilo’s, dan staat er 9.8244 terwijl er aan de dieren niets veranderd is.
  4. Onjuist, want de samenhang tussen de vinlengte en het gewicht van 342 dieren laat zich niet in één getal samenvatten.

2. De covariantie tussen de vinlengte en het gewicht is 9824.42 millimeter-grammen, en de twee standaardafwijkingen zijn 14.06 millimeter en 801.95 gram — allebei op twee decimalen afgerond. Wat is de correlatie tussen de vinlengte en het gewicht? Je hoeft niet precies te rekenen: twee van de vier vallen af zodra je weet waar een correlatie kán liggen.

  1. .76
  2. .87
  3. 12.25
  4. 698.75

3. Over alle 342 pinguïns is de correlatie tussen de snavellengte en de snaveldiepte r = −.24. Binnen Adélie is hij +.39, binnen Chinstrap +.65 en binnen Gentoo +.64. Welke uitspraak is de beste?

  1. Bij pinguïns gaat een langere snavel samen met een ondiepere snavel; de drie soorten apart zijn te klein om iets tegen die −.24 in te brengen.
  2. Een van de vier getallen moet verkeerd uitgerekend zijn — correlaties met een tegengesteld teken kunnen niet allebei over dezelfde dieren gaan.
  3. Binnen elke soort geldt het verband wél, dus die −.24 is een fout in het bestand.
  4. Het min-teken komt uit het verschil tússen de soorten en niet uit een verband tussen de snavellengte en de snaveldiepte.

4. Kwadrateer je de correlatie tussen de vinlengte en het gewicht, dan staat er r² = .76. Wat zegt die .76?

  1. Ruim driekwart van de verschillen in gewicht loopt langs de vinlengte mee; waaróm staat er niet bij.
  2. Van de 342 pinguïns ligt 76 procent op of vlak bij de lijn.
  3. De vinlengte verklaart 76 procent van het gewicht van een pinguïn.
  4. De samenhang tussen de vinlengte en het gewicht is 76 procent van wat hij maximaal kan zijn.

5. Tussen de vinlengte en het gewicht is de correlatie over alle 342 pinguïns hóger dan binnen elke soort afzonderlijk. Bij de snavels gebeurde iets anders: daar klapte het teken óm. Waar zit het verschil?

  1. Bij de vinnen zijn de drie soorten ongeveer even groot, bij de snavels niet.
  2. Bij de snavels meet je twee dingen aan hetzelfde lichaamsdeel, en dan werken de twee maten elkaar tegen.
  3. Bij de snavels ligt het verschil tussen de soorten dwars op het verband binnen één soort; bij de vinnen ligt het in dezelfde richting.
  4. Bij de vinnen is het verband zo sterk dat het teken niet meer kan omklappen.

6. Bij de 2634 ingevulde vragenlijsten maak je van vijf antwoorden één score. Neem je per persoon het rijgemiddelde, dan is de covariantie tussen de openheid en het neuroticisme 0.1050; neem je de rijsom, dan is hij 2.6253. Wat gebeurt er met de correlatie tussen die twee?

  1. Die wordt 25 keer zo groot, net als de covariantie.
  2. Die blijft .16, want delen door de twee standaardafwijkingen haalt de maatstok eruit.
  3. Die wordt 5 keer zo groot, want een som van vijf items is vijf keer een gemiddelde.
  4. Dat kun je niet zeggen zonder de twee standaardafwijkingen erbij.

1 c · 2 b · 3 d · 4 a · 5 c · 6 b

Hieronder staat per vraag waaróm, en waarom de drie andere opties aantrekkelijk zijn. Heb je er een goed, dan kun je die overslaan.

1. c. Dit is de kilo-proef van Jouw beurt, nu als keuze: hetzelfde stel dieren, een ander getal, alleen doordat je in kilo’s weegt. a is waar en geeft antwoord op een andere vraag — of het verband buiten deze 342 dieren geldt, weet je inderdaad niet, maar de kritiek gaat over het woord sterk en niet over de stap naar buiten. b heeft het goede bereik bij de verkeerde maat: −1 tot +1 is het bereik van de correlatie, en juist de covariantie kent zo’n bereik niet. d is de veilige reflex: één getal mág een samenhang samenvatten, dat is precies waar de correlatie voor bestaat.

2. b. Delen door de twee standaardafwijkingen: 9824.42 / (14.06 × 801.95) = .87. En je hoefde het niet uit te rekenen: een correlatie ligt tussen −1 en +1, dus c en d vallen meteen af. c is wat je krijgt als je alleen door de standaardafwijking van het gewicht deelt: 9824.42 / 801.95 = 12.25. d is de andere helft van diezelfde vergissing, alleen door de standaardafwijking van de vinlengte: 9824.42 / 14.06 = 698.75. Dat allebei zo ver buiten het bereik vallen ís de aanwijzing: er ontbreekt een deling. Blijft a over, en dat is r² — het kwadraat, niet de correlatie zelf.

3. d. De Gentoo’s hebben lange, ondiepe snavels en de andere twee soorten korte, diepe; dat soortverschil kantelt de hele wolk omlaag. a neemt het getal over de hele hoop als de uitspraak over een pinguïn, en dat is precies de zin die dit blok verbiedt. b is de reactie waar je jezelf de schuld geeft: de vier getallen kloppen alle vier, ze gaan alleen over verschillende hopen. c heeft de goede helft en de verkeerde conclusie — de −.24 is netjes uitgerekend en geen fout; hij vat iets anders samen dan je dacht.

4. a. Boekhouding: welk deel van de verschillen de twee maten delen, niet waaróm. b verwart een aandeel in de verschillen met een aandeel in de dieren; r² zegt niets over hoeveel pinguïns er dichtbij liggen. c is de verleidelijkste, en er staat één woord te weinig: het gaat over de verschillen in gewicht en niet over het gewicht. Een pinguïn van 4200 gram bestaat niet voor 76 procent uit vinlengte. d maakt van r² een percentage van het maximum; dat zou betekenen dat r = .87 en r² = .76 twee manieren zijn om hetzelfde te zeggen, en dan had je er maar één nodig.

5. c. Zelfde derde speler, twee uitkomsten: ligt het soortverschil dwars op het verband binnen een soort, dan klapt het teken om; ligt het mee, dan wordt de samenhang over alles héén juist sterker dan binnen één soort. a is na te rekenen en niet waar: de drie soorten zijn 151, 68 en 123 dieren, bij de vinnen net zo goed als bij de snavels. b klinkt als een meetkundig argument maar zegt niets — binnen elke soort gaat een langere snavel wél samen met een diepere. d maakt van sterkte een slot op het teken; de snavels laten zien dat een soortverschil een teken kan omzetten hoe sterk het verband binnen de soorten ook is.

6. b. De correlatie blijft .16. Een som is vijf keer een gemiddelde, en dat zit in de teller én in de noemer, dus het valt weg — dezelfde wet als bij de kilo’s van de pinguïns. a trekt de factor 25 van de covariantie door naar de correlatie, en dat is juist wat de correlatie níét doet. c heeft het goede getal bij de verkeerde maat: vijf keer geldt voor één variabele, en in een covariantie zitten er twee, dus vijf maal vijf. d is te bescheiden — je hoeft de standaardafwijkingen niet te kennen om te weten dat ze evenredig meeschuiven.

Alle getallen in deze zes vragen heb je eerder op deze bladzij zelf uitgerekend, bij Jouw beurt. De enige die nieuw zijn, zijn de foute afleiders — en die staan hierboven allemaal als deling voluit, zodat je ze kunt narekenen zonder ze te geloven.

◠ Zelfde vorm in het boek

Hier: één getal voor de samenhang, en de wolk die verklapt wat het getal verzwijgt.

In het boek: W4 — Hangen ze samen?

Daar begint het bij de wolk, niet bij het getal — en de rechte lijn komt pas als je hebt gezien wie er allemaal in de wolk zitten.

Voetnoten

  1. Palmer Penguins — Gorman, Williams & Fraser (2014), Palmer Station Antarctica LTER; via het R-pakket palmerpenguins (CC0). Alle getallen in dit blok zijn op die échte data nagerekend.↩︎