Я хотел бы закодировать уравнение в R, которое рассчитывает критерии регрессии хребта

Question

Критерий определяется следующей функцией на прилагаемом изображении. Я работал над этим некоторое время, но это то, что у меня есть. Мне нужно

ridge.crit <- function(y, X, lambda, beta){
n <- length(y) #number of rows
p <- ncol(X) #number of predictor variables
  for(i in 1:n){
    (y[i] - beta[1] - sum(beta[2:p+1] * X[i,]))^2
##WE GOTTA ADD THEM TOGETHER AND STORE FINAL ANSWER AND RETURN IT!
    }
  return(final_answer)
}