07.01.2013 Aufrufe

Grafiken und Statistik in R

Grafiken und Statistik in R

Grafiken und Statistik in R

MEHR ANZEIGEN
WENIGER ANZEIGEN

Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.

YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.

y=<br />

β0+β1·x1+<br />

β2·x2+...+ɛ<br />

y<br />

160 170 180 190 200<br />

●<br />

●<br />

●<br />

●<br />

● ●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

20 30 40 50 60 70<br />

x<br />

y=<br />

β0+β1·x n<br />

1 +<br />

β2·x2+...+ɛ<br />

y<br />

160 170 180 190 200<br />

●<br />

●<br />

●<br />

●<br />

● ●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

20 30 40 50 60 70<br />

x<br />

4.3 Regressionsanalyse<br />

y<br />

160 170 180 190 200<br />

●<br />

●<br />

●<br />

●<br />

● ●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

●<br />

20 30 40 50 60 70<br />

x<br />

curve( # L<strong>in</strong>ie Konfidenz<strong>in</strong>tervall<br />

predict(modell,<br />

newdata = data.frame(x = x),<br />

<strong>in</strong>terval = "confidence")[ , "lwr"], # lower Intervall<br />

add = TRUE,<br />

lty="dashed") # L<strong>in</strong>ientyp s. 29<br />

curve( # L<strong>in</strong>ie Konfidenz<strong>in</strong>tervall<br />

predict(modell,<br />

newdata = data.frame(x = x),<br />

<strong>in</strong>terval = "confidence")[ , "upr"],<br />

add = TRUE,<br />

lty="dashed") # L<strong>in</strong>ientyp s. 29<br />

4 <strong>Statistik</strong><br />

� Vorhersage<strong>in</strong>tervalle analog mit <strong>in</strong>terval = "prediction"; predict() <strong>und</strong> curve() funktionieren auch mit anderen (nichtl<strong>in</strong>earen)<br />

Modellgleichungen<br />

4.3.2 L<strong>in</strong>ear, multipel<br />

E<strong>in</strong> Modell der Form y = β0 + β1 · x1 + β2 · x2 + ... + ɛ<br />

data(airquality) # Daten e<strong>in</strong>lesen<br />

?airquality # Hilfe zum Datensatz<br />

names(airquality) # welche Variablen?<br />

attach(airquality) # wenn attach(...) dann entfällt 'data=airquality'i.d.R.<br />

ozone.lm

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!