Skip to content

Group lasso

How a block-norm penalty selects predefined predictor groups, with exact group soft thresholding, a Lasso comparison and the limits of orthogonal-design arithmetic.

Group lasso is a regression penalty that can set predefined groups of coefficients to zero together. For a quant modelling a return series with several lags of each economic signal, it offers a way to retain or discard a whole signal family. Ordinary Lasso regression selects individual coefficients instead.

The unit of selection is a block

Partition the predictors into non-overlapping groups g=1,…,Gg=1,\ldots,G. For centred responses and predictors, one convention is

β^=arg min⁡β{12n∥y−Xβ∥22+λ∑g=1Gwg∥βg∥2},λ≥0,wg>0.\widehat\beta=\operatorname*{arg\,min}_{\beta} \left\{\frac{1}{2n}\lVert y-X\beta\rVert_2^2 +\lambda\sum_{g=1}^G w_g\lVert\beta_g\rVert_2\right\}, \qquad \lambda\geq0,\quad w_g>0.
Here nn is the number of observations, yy the response vector and XX the predictor matrix. The vector βg\beta_g contains the coefficients assigned to group gg; an intercept, if fitted separately, is normally unpenalised. The penalty sums Euclidean norms, not squared norms. Its kink when a whole block reaches zero creates group sparsity. Single-coefficient groups recover the L1 penalty. Yuan and Lin, §2.

Weights matter: wg=pgw_g=\sqrt{p_g}, where pgp_g is group size, is a common convention. Unit weights are another choice. Changing weights changes the fit. So does changing the loss normalisation without adjusting λ\lambda: replacing 1/(2n)1/(2n) by 1/21/2 requires multiplying λ\lambda by nn to preserve this objective's minimiser.

Two groups, one exact calculation

Consider the following synthetic, dimensionless proximal problem, with unit weights:
min⁡β  12∥z−β∥22+2(∥β1∥2+∥β2∥2),z1=(3,4),z2=(0.5,0.5).\min_\beta\;\frac12\lVert z-\beta\rVert_2^2 +2\bigl(\lVert\beta_1\rVert_2+\lVert\beta_2\rVert_2\bigr), \qquad z_1=(3,4),\quad z_2=(0.5,0.5).

This is the orthogonal-design calculation, rather than a general correlated regression fit. With loss 12∥y−Xβ∥22\frac12\lVert y-X\beta\rVert_2^2, it follows when X⊤X=IX^\top X=I and z=X⊤yz=X^\top y, after dropping a constant. It also appears as a proximal-gradient update, with the step size incorporated into the penalty parameter.

For any block, the exact rule is block soft thresholding, applied separately because the groups do not overlap:
β^g={0,∥zg∥2≤λwg,(1−λwg∥zg∥2)zg,∥zg∥2>λwg.\widehat\beta_g= \begin{cases} 0, & \lVert z_g\rVert_2\leq\lambda w_g,\\ \left(1-\dfrac{\lambda w_g}{\lVert z_g\rVert_2}\right)z_g, & \lVert z_g\rVert_2>\lambda w_g. \end{cases}

The first branch includes zg=0z_g=0, even when λ=0\lambda=0; never divide by a zero norm. To see the mechanism, fix the output norm rr. The closest vector of norm rr points along zgz_g, leaving the scalar problem 12(∥zg∥2−r)2+λwgr\frac12(\lVert z_g\rVert_2-r)^2+\lambda w_gr for r≥0r\geq0. Its minimiser is r=max⁡(∥zg∥2−λwg,0)r=\max(\lVert z_g\rVert_2-\lambda w_g,0).

BlockInput normCalculationOutput
Group 155(1−2/5)(3,4)(1-2/5)(3,4)(1.8,2.4)(1.8,2.4)
Group 20.5≈0.707\sqrt{0.5}\approx0.707Norm is below 22(0,0)(0,0)
Group 1's output norm is 33; both its coordinates shrink by the same factor, 0.60.6. With the same quadratic loss and λ=2\lambda=2, coordinatewise Lasso instead returns (1,2)(1,2) for its first two coefficients. These are different penalties; the common numerical λ\lambda does not make their regularisation strengths universally equivalent.
Group 1 moves from input (3,4) to (1.8,2.4) along the same ray, while coordinatewise Lasso gives (1,2). Group 1 input norm 5 exceeds threshold 2; Group 2 input norm approximately 0.707 is below 2 and gives the zero vector.
Group 1 moves from input (3,4) to (1.8,2.4) along the same ray, while coordinatewise Lasso gives (1,2). Group 1 input norm 5 exceeds threshold 2; Group 2 input norm approximately 0.707 is below 2 and gives the zero vector.
Original synthetic illustration with unit group weights and dimensionless coefficients. The left panel compares directions; the right compares input group norms with the threshold.

What the grouping commits you to

Groups must be supplied before fitting; the method does not discover economic families. A lag block, categorical encoding or spline basis can be a sensible unit, but a nonzero group does not establish that every component is informative. Ordinary group lasso does not encourage sparsity within a surviving group. Sparse-group lasso adds an L1 term for that separate task. Simon et al., §§1–2.
The norm is unchanged by an orthogonal rotation within a group, so a corresponding rotation of the design and coefficients preserves the criterion. Arbitrary rescaling does not. Choose a scaling convention, estimate standardisation from training data only, and preserve it in validation. Standardising columns does not make them orthogonal: the one-shot formula above cannot simply be applied to marginal coefficients from a correlated design.
For the ordinary-penalty foundations and a worked example of training-only scaling, see Ridge, Lasso, and Elastic Net (public preview; full lesson requires membership).

Sources

Last updated
  • group lasso
  • regularisation
  • block sparsity
  • feature selection
← Back to the glossaryBack to top ↑