this piece of code will return coefficients :intercept , slop1 , slop2
set.seed(1) n=10 y=rnorm(n) x1=rnorm(n) x2=rnorm(n) lm.ft=function(y,x1,x2) return
You need predict -
predict
set.seed(1) n=10 y=rnorm(n) x1=rnorm(n) x2=rnorm(n) lm.ft=function(y,x1,x2) # return(lm(y~x1+x2)$coef) return(lm(y~x1+x2)) res=lm.ft(y,x1,x2) ypredicted <- predict(res) residuals <- y - ypredicted