Regressie en correlatie
Oefening verband lengte – gewicht van volwassen persoon
#reg en cor: verband lengte - gewicht
rm(list=ls(all=TRUE))
library(units)
lengte <- set_units(c(170, 179, 185, 171, 188, 175, 190, 166, 186,
192), cm)
gewicht <- set_units(c(73, 74, 81, 74, 82, 70, 82, 60, 90, 85), kg)
yonx <- lm(gewicht ~ lengte)
xony <- lm(lengte ~ gewicht)
print(coef(yonx))
print(coef(xony))
plot(x = lengte, y = gewicht, main= "Verband lengte & gewicht")
,abline(yonx, lwd = 2, col = "red")
abline(a = -coef(xony)[1]/coef(xony)[2], b = 1/coef(xony)[2], lwd =
2, col = "blue")
Geen leerstof (te verwarrend)
print(paste("y=", signif(coef(yonx)[1], 4), "+",
signif(coef(yonx)[2], 4), "*x"))
print(paste("x=", signif(coef(xony)[1], 4), "+",
signif(coef(xony)[2], 4), "*y"))
print(cor(lengte, gewicht))
=> r = correlatiecoëfficiënt
Als |r| > 0,95 =aangetoond, bruikbaar verband
#Afwerken van plot met legende
legend("topleft", legend = c("regressielijn van y op x: gewicht
i.f.v. lengte",
"regressielijn van x op y: lengte i.f.v. gewicht"),
lwd = 2, col = c("red", "blue"))
,