Deutsches R-Forum

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)