回归：用户定义的函数

Question

这是一个简单的线性回归函数，用于计算β值，残差和拟合值。如何为此函数添加多个预测变量？

regression=function(x,y){
 x=cbind(1,x)
 beta=solve(t(x)%*%x)%*%t(x)%*%y
 residuals=y-x%*%beta
 fitted.values=x%*%beta
 res=list(coefficients=beta,residuals=residuals,fitted.values=fitted.values)
 class(res)="regression" 
 return(res)
}