powered by
Compute the Frobenius penalty for all the parameters in the model
matrix.penalty(U, penalty)