Zunächst muss der Datensatz eingelesen werden:
data(women)
Über das Kommando 'data' werden mitgelieferte Datensätze in eine R Sitzung geladen.
Über viele Datensätze gibt ausführliche Informationen, die als Hilfeseiten mitgeliefert werden.
help(women)
Einige Kommandos lassen den Datensatz genauer erkunden.
names(women)
## [1] "height" "weight"
str(women)
## 'data.frame': 15 obs. of 2 variables:
## $ height: num 58 59 60 61 62 63 64 65 66 67 ...
## $ weight: num 115 117 120 123 126 129 132 135 139 142 ...
Auch Bilder lassen sich in ein solches Dokument einbetten. als Beispiel soll der Zusammenhang zwischen Größe und Gewicht der Individuen angeschaut werden.
plot(women, xlab = "Height (in)", ylab = "Weight (lb)", main = "women data: American women aged 30-39")
Das Bild lässt einen linearen Zusammenhang vermuten. Das lässt sich leicht mit einer Regression überprüfen.
lineareRegression <- lm(women$weight ~ women$height)
summary(lineareRegression)
##
## Call:
## lm(formula = women$weight ~ women$height)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.733 -1.133 -0.383 0.742 3.117
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -87.5167 5.9369 -14.7 1.7e-09 ***
## women$height 3.4500 0.0911 37.9 1.1e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 1.53 on 13 degrees of freedom
## Multiple R-squared: 0.991, Adjusted R-squared: 0.99
## F-statistic: 1.43e+03 on 1 and 13 DF, p-value: 1.09e-14
Die Variable 'lineareRegression' enthält das vollständige Ergebnis einer Regression in R. Einige Kommandos wissen, wie mit Objekten verschiedenen Typs umzugehn ist, z.B. 'plot'.
plot(lineareRegression)
Ist die Regression angemessen? Na ja.
Jedenfalls bekommt man mit *.Rmd Dokumenten eine Möglichkeit leicht lesbare Dokumente zu erzeugen.