Statistiek-1 + R

Staafdiagram

» Start

Staafdiagram.


Voor het maken van een staafdiagram moeten we een kruistabel [met frekwenties] aanbieden aan de functie barplot. Zo'n tabel maak je met de functie table. Bijvoorbeeld voor de frekwenties van voorkomen van de verschillende oogkleuren in de variabele ogen uit de datamatrix die we zelf hebben samengesteld.:


> 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" 

   

    > table(ogen):

    ogen

    blauw bruin grijs groen zwart

    26    13     5    21     1 


Dat is dus een kruistabel met 1 rij, met frekwenties voor oogkleur, waarmee de functie barplot dus direct een staafdiagram kan maken.


> barplot(table(ogen), beside=TRUE)

bar-ogen.pdf

Als je naar het voorbeeld kijkt in de Rweb graphics tutorial, dan wordt daarin de datamatrix VADeaths uit het pakket datasets gebruikt.

Als je VADeaths uitprint zie je dit:


    > VADeaths

        Rural Male Rural Female Urban Male Urban Female

    50-54       11.7          8.7       15.4          8.4

    55-59       18.1         11.7       24.3         13.6

    60-64       26.9         20.3       37.0         19.3

    65-69       41.0         30.9       54.6         35.1

    70-74       66.0         54.3       71.1         50.0


Dus een kruistabel met een klasse-indeling in de rijen, en 4 kolommen met frekwenties voor 2 x 2 combinaties van Rural, Urban, Male, Female.


Ook in dit voorbeeld is de input voor de functie barplot dus een kruistabel, maar dan wel een met 4 rijen. Dat kan dus 4 staafdiagrammen opleveren.


Gelet op dit voorbeeld, moeten we dus eerst een kruistabel zien te maken voor oogkleur-frekwenties, uitgesplitst op man/vrouw, dus een kruistabel met 2 rijen [ipv 1], voor dat we het staafdiagram kunnen plotten. Hoe doen we dat?


    > table(cohort[,c(4,6)])

             ogen

    geslacht blauw bruin grijs groen zwart

        m    15     6     2    11     0

        v    11     7     3    10     1


Als R opdracht staat daar dus:

Maak een frekwentietabel voor de gegevens in de 4e en 6e kolom van het data-frame met de naam data.


En dit commando levert dan het bijbehorende staafdiagram op:


    > barplot(table(cohort[,c(4,6)]), beside=TRUE)

staaf-beside.pdf


Bron:

Buijs, A. - Statistiek om mee te werken. Stenfert Kroese, Groningen (2003)

Maindonald, J. & J. Brown. - Data Analysis and Graphics Using R. Cambridge University Press, Cambridge (2003)