1 Statistische Auswertungen mit R Universität Kassel, FB 07 ...

Statistische Auswertungen mit R Universität Kassel, FB 07 Wirtschaftswissenschaften
Dipl.-Volksw. Markus Pullen Wintersemester 2012/13
Übung 5. Sitzung:
1. Plotten Sie den BMI (Hochachse) gegen die Körpergröße (Rechtsachse) und ergänzen Sie die
entsprechende Regressionsgerade.
Gestalten Sie die Graphik aussagekräftiger.
> attach(groe.gew)
> plot(groesse, BMI, main = "BMI und Körpergröße", sub = "ACHTUNG:
FIKTIVE DATEN", xlab="Größe", ylab="Body-Mass-Index", las = 1,
xlim=c(0,200) )
> # Es stellt sich allerdings die Frage, inwiefern es sinnvoll ist,
die Kennzahl BMI, in deren Berechnung die Größe „Größe“ eingeflossen
ist, gegen die Größe zu plotten. Dies liefert möglicherweise eher
Informationen über die mathematischen Eigenschaften der BMI-Funktion
als über die betrachteten Merkmalsträger. Auch die Ermittlung einer
Regressionsgeraden ist entsprechend kritisch zu sehen. Bitte
betrachten Sie diese Aufgabenstellung in erster Linie als
„Fingerübung“.
Body-Mass-Index
32
30
28
26
BMI und Körpergröße
0 50 100 150 200
Größe
ACHTUNG: FIKTIVE DATEN
> BMI.gerade = lm ( BMI ~ groesse )
> BMI.aab = BMI.gerade $ coeff [1]
> BMI.stg = BMI.gerade $ coeff [2]
1

plot(groesse, BMI, main = "BMI und Körpergröße", sub = "ACHTUNG:
FIKTIVE DATEN", xlab="Größe", ylab="Body-Mass-Index", las = 1,
xlim=c(0,200), ylim=c(0,35), pch=as.character(geschl) )
> curve ( BMI.aab + BMI.stg * x, 150:200, add=TRUE , col="GREEN")
Body-Mass-Index
35
30
25
20
15
10
5
0
BMI und Körpergröße
Erzeugen Sie entsprechende Graphiken, getrennt nach Frauen und Männern.
> attach (teildaten $ m)
The following object(s) are masked from 'groe.gew':
BMI, geschl, gewicht, groesse
> plot (groesse, BMI, main = "BMI und Körpergröße", sub = "ACHTUNG:
FIKTIVE DATEN", las=1, xlim = c(0,200), ylim = c(0, 35), pch =
as.character(geschl) )
> BMI.gerade.m = lm ( BMI ~ groesse)
2
w
m
w m
m
w m
0 50 100 150 200
Größe
ACHTUNG: FIKTIVE DATEN
w
w
m

abline (coef = BMI.gerade.m $ coeff, col = "blue")
BMI
> detach (teildaten $ m)
> attach (teildaten $ w)
The following object(s) are masked from 'groe.gew':
BMI, geschl, gewicht, groesse
> plot (groesse, BMI, main = "BMI und Körpergröße", sub = "ACHTUNG:
FIKTIVE DATEN", las=1, xlim = c(0,200), ylim = c(0, 35), pch =
as.character(geschl) )
> BMI.gerade.w = lm ( BMI ~ groesse)
> abline (coef = BMI.gerade.w $ coeff, col = "red")
BMI
35
30
25
20
15
10
5
0
35
30
25
20
15
10
5
0
BMI und Körpergröße
0 50 100 150 200
groesse
ACHTUNG: FIKTIVE DATEN
BMI und Körpergröße
0 50 100 150 200
groesse
ACHTUNG: FIKTIVE DATEN
3
w
w
m
w
w
m
w
m
m
m