x1<-sample(1:100) x2<-sample(1:100) y<-2*x1+0.01*x2+rnorm(100) p<-data.frame(y,x1,x2) lm(y~x1+x2,data=p) f<-function(r){ sum((y-r[1]*x1-r[2]*x2)^2) } optim(c(1,3),f) ##Lasso f2<-function(r){ sum((y-r[1]*x1-r[2]*x2)^2)+20000*sum(abs(r)) } o…