Kwantielplot.
Kwantielplots, of qq-plots, zijn handig om, bijvoorbeeld, een indruk te krijgen van de normaliteit van een verdeling [van waarden in een steekproef].
Eerst de data binnen halen:
> bron <- "http://www.mzandee.net/~zandee/statistiek/data/gegevens.txt"
> cohort <- read.tabel(bron, header=T)
> attach(cohort)
> names(cohort)
[1] "lichaam" "arm" "pols" "geslacht" "hand" "ogen"
Vervolgens een histogram maken van de variabele lichaam, en de bijbehorende dichtsheidfunctie:
> dens <- density(lichaam) # bereken de dichtheidsfunctie
> xlim <- range(dens$x) # gebruik de dichtheid om grenswaarden voor x-as te bepalen
> ylim <- range(dens$y) # idem, y-as
> hist(lichaam, breaks=150+(0:10)*5, xlab="Lichaamslengte", probability=T, xlim=xlim, ylim=ylim)
> lines(dens)
Ziet er redelijk normaal uit. Misschien ietsje scheef met een staartje aan de linker kant.
> qqnorm(lichaam)
> qqline(lichaam, col=2)
