BoxPlot.
Een 'box- en whisker-plot' geeft een grafische samenvatting van een aantal eigenschappen van een steekproef. Je ziet de mediaan, het eerste en derde kwartiel (de grenzen van de 'box': daarbinnen ligt 50% van de waarnemingen), en de maximum en minimumwaarde (de 'whiskers'), en, als ze voorkomen ook de uitschieters. Ook kun je met deze plot beoordelen of de verdeling al of niet scheef is [ligt de mediaan in het midden van de box?]
> bron <- "http://www.mzandee.net/~zandee/statistiek/data/gegevens.txt"
> cohort <- read.table(bron, header=T)
> attach(cohort)
> names(cohort)
[1] "lichaam" "arm" "pols" "geslacht" "hand" "ogen"
> boxplot(lichaam)
Lichaamslengte van cohort studenten
Een boxplot is vooral handig als je twee of meer steekproeven met elkaar wilt vergelijken.
> boxplot(lichaam ~ geslacht)
Lichaamslengte van cohort studenten, gegroepeerd op sexe (m/v)
Een dergelijk plaatje houdt echter nog geen formele toetsing in op een eventueel verschil in gemiddelde. Daarvoor moeten we een t-toets gebruiken. Wel kunnen we in een boxplot zg inkepingen ('notches') aanbrengen die de grenzen van het 95% vertrouwens interval van het gemiddelde afbeelden. We krijgen daarmee wel een duidelijke hint dat bij formele toetsing het verschil significant zal blijken te zijn (geen overlap in vertrouwensintervallen).
> boxplot(lichaam ~ geslacht, notch=T)
We kunnen ook het effect van meer dan twee onafhankelijke groeperingen tegelijk afbeelden, zoals sexe en links/rechts-handigheid.
> boxplot(lichaam ~ geslacht + hand, notch=T)
We zien dat het betrouwbaarheidsinterval van de gemiddelde lichaamslengte van linkshandige vrouwen zo groot is (door de geringe steekproefomvang) dat het de buiten de waargenomen hoogste en laagste waarneming in de steekproef valt.
BoxPlot mbv lattice pakket:
> library(lattice)
> bwplot(lichaam ~ geslacht + hand, data = cohort,
+ ylab = "Lengte (cm)",
+ panel = function(...) {
+ panel.grid(v = 0, h = -1)
+ panel.bwplot(...)
+ },
+ par.settings = list(plot.symbol = list(pch = 4)))
Bron:
Website van John Maindonald en John Braun's boek : Data Analysis and Graphics Using R - An Example-Based Approach, met plot voorbeelden en R-code.
