Code: Alles auswählen
meinreg <- function (X, Y) {
x <- X - mean(X); y <- Y - mean(Y)
Anstieg <- sum(x*y) / sum(x^2)
c(mean(Y) - Anstieg*mean(X), Anstieg)
}
with(swiss, meinreg(Education, Fertility))
lm(Fertility ~ Education, data=swiss)
Moderator: EDi
Code: Alles auswählen
meinreg <- function (X, Y) {
x <- X - mean(X); y <- Y - mean(Y)
Anstieg <- sum(x*y) / sum(x^2)
c(mean(Y) - Anstieg*mean(X), Anstieg)
}
with(swiss, meinreg(Education, Fertility))
lm(Fertility ~ Education, data=swiss)