REVIEW 4 major objections 5 minor 29 references
A Short Note of Comparison between Convex and Non-convex Penalized Likelihood
T0 review · 4 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read The paper argues that non-convex penalties such as SCAD and MCP can reach ground-truth sparsity without sacrificing estimation accuracy, while Lasso trades one against the other.
desk verdict A correct but unoriginal teaching note whose own experiments don't back the headline claim. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the first derivative of the penalty function $p_\lambda(t)$: Lasso has $p_\lambda'(t)=\lambda$ constant, SCAD has $p_\lambda'(t)=\lambda[I(t\le\lambda)+((a\lambda-t)_+/((a-1)\lambda))I(t>\lambda)]$, and MCP has $p_\lambda'(t)=(\lambda-t/a)_+$. This derivative determines the shrinkage pattern: a constant derivative means equal shrinkage of all coefficients, while a derivative that declines to zero means only small coefficients are penalized. The paper reads the bias, sparsity, and accuracy comparison off this derivative, then uses coordinate descent with closed-form soft-thresholding-type updates for each penalty to test the behavior computationally.
What would settle it
Run the n=200, p=1000 simulation with SCAD and MCP over a grid of shape parameters (e.g., $a=2.1,3.7,10$) and several random initializations; if some runs stop at visibly higher loss or sparsity levels far from 10, the claimed sparsity-without-accuracy-loss is an optimizer artifact rather than a penalty property. A cleaner check is to solve small instances to certified global optimality and compare the coordinate-descent solution against the global solution.
Extended reading notes
Core claim
The central claim is that non-convex penalties are not merely a theoretical alternative to Lasso: in a high-dimensional linear model with n=200, p=1000 and 10 nonzero coefficients, SCAD and MCP produce estimators with ground-truth sparsity level without sacrificing estimation accuracy, while Lasso requires a much larger penalty to reach that sparsity and then suffers high MSE. The mechanism is the derivative of the penalty: Lasso applies a constant penalty to every coefficient regardless of size, shrinking large coefficients and creating bias, whereas SCAD and MCP keep a positive penalty only for small values and let it drop to zero for large values, so large coefficients are left nearly unpenalized. The note derives closed-form coordinate-descent updates for all three penalties and reports that SCAD and MCP converge about three times slower than Lasso with unstable loss variation across iterations. On the rat eye gene-expression data, under cross-validated tuning, SCAD yields a sparser model with lower residual sum of squares than Lasso.
Load-bearing premise
The load-bearing premise is that the coordinate-descent solutions for SCAD and MCP show the penalties' true behavior rather than a quirk of the optimizer: Section 3 shows unstable loss variation and does not report the shape parameter, the convergence tolerance, or the starting values, so if the optimizer settles in a poor local optimum the claimed sparsity without accuracy loss could be an algorithmic accident.
Editorial extensions
If this is right
- In high-dimensional sparse linear models, Lasso forces a trade-off: a penalty level that recovers true sparsity over-shrinks the nonzero coefficients and inflates MSE.
- With SCAD or MCP, a single penalty level can deliver both the ground-truth sparsity level and low estimation error in the simulated setting.
- Non-convex penalties carry a computational cost: coordinate descent takes roughly three times as long on the toy problem and shows unstable loss variation.
- On the rat gene-expression data, SCAD gives a sparser model with lower residual sum of squares than Lasso under cross-validated tuning.
- The practical guide is to prefer Lasso when fast convex optimization matters and SCAD or MCP when bias on large coefficients matters.
- The choice between convex and non-convex penalties is not about one being universally better; it depends on whether accuracy or computational feasibility is the binding constraint.
Reading between the lines
- A natural extension of the derivative-based explanation is that any penalty whose derivative reaches zero for large coefficients should reproduce the sparsity-without-bias pattern, while any constant-derivative penalty should inherit Lasso's bias.
- Because Section 3 omits the shape parameter $a$, convergence tolerance, and initialization, rerunning the simulation across $a$ values and random starts would separate penalty behavior from optimizer behavior.
- The real-data conclusion rests on a single data set and a single cross-validation run; repeating the comparison with repeated CV splits or multiple data sets would show whether the sparser SCAD fit is stable.
- A practical extension the paper does not test is initializing SCAD or MCP from a Lasso solution, which could reduce the reported convergence time while keeping the favorable sparsity-accuracy profile.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This short note compares Lasso with the non-convex SCAD and MCP penalties in high-dimensional penalized likelihood estimation. The paper reviews the penalties, states their known bias/oracle properties, presents a small simulation study (Section 2), derives closed-form coordinate descent updates (Section 3), and reports a real-data comparison on a rat eye eQTL dataset (Section 4). The central empirical assertion is that SCAD and MCP can achieve the ground-truth sparsity level without sacrificing estimation accuracy, while Lasso cannot simultaneously satisfy both goals.
Significance. If the claims were fully supported, the paper would offer a useful, concise comparison with practical guidance on choosing between convex and non-convex penalties. The qualitative statements accord with the existing literature, and the coordinate descent formulas are standard. However, the paper provides no new theory, and its own empirical evidence is too incomplete to establish the central claim: the simulation lacks paired sparsity/MSE statistics, the non-convex optimization settings are unreported, and the real-data comparison relies on in-sample RSS. The manuscript is a reasonable pedagogical note, but the load-bearing quantitative assertions need substantial strengthening before the paper can be accepted as a standalone research contribution.
major comments (4)
- [Section 2, after Figure 3] The claim that 'the SCAD and the MCP penalties can give us an estimator with ground truth sparsity level without sacrificing estimation accuracy' is not established by the displayed results. Figure 3 plots MSE and sparsity level separately against lambda; for Lasso the text explicitly identifies lambda=0.12 (minimum MSE) and lambda>0.35 (sparsity 10), but no lambda is identified for SCAD or MCP at which the average sparsity equals 10 and the MSE is simultaneously low. Because the two criteria are not paired at a single lambda, the curves only suggest that some intermediate lambda may look favorable. Please report a paired summary: for each method and each lambda, the average sparsity and MSE together, or identify a concrete lambda where both objectives are met, with standard errors across the 100 replications.
- [Section 3, coordinate descent formulas] The closed-form updates for Lasso, SCAD, and MCP are not correctly normalized with respect to the objective in Eq. (1.1), which includes the factor 1/(2n). With L(beta)=||Y-Xbeta||^2, the coordinate-wise solution for Lasso is beta_j = S(X_j^T r, n lambda)/X_j^T X_j (or, if columns are standardized to unit L2 norm, S(X_j^T r, lambda)). The manuscript writes beta_j = S(X_j^T r, lambda) without stating the column normalization. For the real data, variables are standardized to zero mean and unit standard deviation, so X_j^T X_j is approximately n, making the displayed threshold inconsistent with Eq. (1.1). Please clarify the normalization used in the simulations and correct the formulas, since the lambda grid in Figure 3 depends on this scaling.
- [Section 3, Figure 4 and Table 1] The empirical comparison of convergence is not reproducible because the optimization settings are unreported. The paper never states the SCAD shape parameter a, the MCP shape parameter a, the initialization, whether warm starts were used, the convergence tolerance, the maximum number of iterations, or the lambda grid. For non-convex objectives, coordinate descent can converge to different local minima depending on these choices; the unstable loss variation in Figure 4 may therefore reflect the algorithm's terminal points rather than a property of the penalty family. Please report all of these settings and consider a sensitivity check with multiple starting values or a comparison against a well-tested implementation such as ncvreg.
- [Section 4, Figure 5] The real-data comparison uses the residual sum of squares (RSS) to conclude that SCAD 'performs better' than Lasso. Since RSS is an in-sample measure and the two methods have different sparsity levels under the selected lambdas, a lower RSS does not establish better predictive performance or better model selection. Please use out-of-sample prediction error, cross-validated deviance, or another criterion that accounts for model complexity. In addition, the real-data section omits MCP entirely, so it does not support the paper's stated three-way comparison.
minor comments (5)
- [Section 2] The sentence 'the Lasso penalty always introduces a bias on its large components without sacrificing model selection consistency' is ambiguous and potentially misleading: model selection consistency of Lasso is known to require additional conditions such as the irrepresentable condition. Please rephrase to state the precise theoretical result.
- [Section 1, Eq. (1.1)] The scaling 1/(2n) is introduced for the least-squares case, but L(beta) is defined generally as the negative log-likelihood. Please clarify whether Eq. (1.1) is meant to be generic or specific to the linear model.
- [Section 3, Figure 4] The quantity 'Loss Variation' plotted in Figure 4 is not defined. Specify whether it is |Q(beta^(k)) - Q(beta^(k-1))|, a relative change, or something else.
- [Throughout] There are several typographical issues: 'coordinate decent' should be 'coordinate descent', 'right penal' should be 'right panel', 'BardetBiedl' should be 'Bardet-Biedl', and 'Not surprised' should be 'It is not surprising'.
- [References] The paper should cite relevant algorithmic references for coordinate descent with SCAD and MCP, such as Breheny and Huang (2011) and Mazumder et al. (2011), and the oracle property results of Fan and Li (2001) and Zhang (2010).
Circularity Check
No significant circularity: the comparison is expository, the penalty formulas are standard definitions, and no fitted parameter, self-citation chain, or definitional reduction forces the paper's conclusion.
full rationale
The paper's derivation chain is self-contained in the relevant sense. The penalized likelihood objective (1.1) and the Lasso/SCAD/MCP penalty functions are standard definitions taken from the cited literature; the coordinate-descent update formulas in Section 3 are derived algebraically from those definitions, not fitted to the conclusion. The qualitative claims in Section 2 (Lasso shrinks large coefficients, SCAD/MCP flatten the penalty for large |t|) follow directly from the plotted derivative functions and are consistent with the externally cited results of Fan and Li (2001) and Zhang et al. (2010). The toy simulation in Model 1 is an illustration: no parameter is fitted to make SCAD/MCP look better, and no 'prediction' is produced from a fit that was calibrated on the same target. The central empirical assertion—that SCAD/MCP can give an estimator with ground-truth sparsity without sacrificing estimation accuracy—rests on visual comparison of the MSE and sparsity curves in Figure 3 rather than a paired statistic, and Section 3 omits the shape parameter a, initialization, and convergence tolerance for the non-convex coordinate-descent runs. These are transparency and evidentiary weaknesses, not circularity, because the curves are generated from the model and penalties rather than being defined to equal the conclusion. There are no load-bearing self-citations, no uniqueness theorem imported from the authors' prior work, and no ansatz smuggled in via citation. The comparison is therefore not circular; the relevant risk is correctness and rigor of the empirical support, which is outside the circularity score.
Assumptions & free parameters
free parameters (4)
- SCAD shape parameter a =
not reported
- MCP shape parameter a =
not reported
- Regularization parameter lambda =
simulation path 0.05 to 0.4; real data Lasso 0.146, SCAD 0.158
- Convergence tolerance and maximum iterations =
not reported
assumptions (4)
- domain assumption The true coefficient vector is exactly sparse (10 non-zero coefficients in simulation)
- domain assumption The data follow a linear model with i.i.d. normal errors and standardized predictors
- domain assumption Coordinate descent on the non-convex objectives reaches a relevant optimum
- domain assumption The oracle properties of SCAD/MCP cited from Fan and Li (2001) and Zhang et al. (2010) apply to the simulation and real-data settings
Cite this review
Pith. "Pith review of A Short Note of Comparison between Convex and Non-convex Penalized Likelihood." pith.science (2026). https://pith.science/paper/I3YJARQZ
@misc{pith2026250207655,
author = {Pith},
title = {Pith review of: A Short Note of Comparison between Convex and Non-convex Penalized Likelihood},
year = {2026},
howpublished = {\url{https://pith.science/paper/I3YJARQZ}},
note = {Machine review of arXiv:2502.07655}
}
read the original abstract
This paper compares convex and non-convex penalized likelihood methods in high-dimensional statistical modeling, focusing on their strengths and limitations. Convex penalties, like LASSO, offer computational efficiency and strong theoretical guarantees but often introduce bias in parameter estimation. Non-convex penalties, such as SCAD and MCP, reduce bias and achieve oracle properties but pose optimization challenges due to non-convexity. The paper highlights key differences in bias-variance trade-offs, computational complexity, and robustness, offering practical guidance for method selection. It concludes that the choice depends on the problem context, balancing accuracy
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Candes, E., Tao, T., et al. (2007). The dantzig selector: Statistical estimation when p is much larger than n. The annals of Statistics , 35(6):2313--2351
work page 2007
-
[2]
Chen, J., de Hoogh, K., Gulliver, J., Hoffmann, B., Hertel, O., Ketzel, M., Bauwelinck, M., Van Donkelaar, A., Hvidtfeldt, U. A., Katsouyanni, K., et al. (2019). A comparison of linear regression, regularization, and machine learning algorithms to develop europe-wide spatial models of fine particles and nitrogen dioxide. Environment international , 130:104934
work page 2019
-
[3]
Ci, H., Liu, M., Pan, X., Zhong, F., and Wang, Y. (2023a). Proactive multi-camera collaboration for 3d human pose estimation. arXiv preprint arXiv:2303.03767
work page Pith review arXiv 2023
-
[4]
Ci, H., Wang, C., Ma, X., and Wang, Y. (2019). Optimizing network structure for 3d human pose estimation. In Proceedings of the IEEE/CVF international conference on computer vision , pages 2262--2271
work page 2019
-
[5]
Ci, H., Wu, M., Zhu, W., Ma, X., Dong, H., Zhong, F., and Wang, Y. (2023b). Gfpose: Learning 3d human pose prior with gradient fields. In Proceedings of the IEEE/CVF conference on computer vision and pattern recognition , pages 4800--4810
work page 2023
-
[6]
Cui, S., Guo, X., Li, R., Yang, S., and Zhang, Z. (2024a). Double-estimation-friendly inference for high dimensional misspecified measurement error models. arXiv preprint arXiv:2409.16463
work page Pith review arXiv 2024
-
[7]
Cui, S., Guo, X., and Zhang, Z. (2024b). Estimation and inference in ultrahigh dimensional partially linear single-index models. arXiv preprint arXiv:2404.04471
work page Pith review arXiv 2024
-
[8]
Cui, S., Li, D., Li, R., and Xue, L. (2024c). Hypothesis testing for high-dimensional matrix-valued data. arXiv preprint arXiv:2412.07987
work page Pith review arXiv 2024
Show all 29 references
-
[9]
Cui, S., Sudjianto, A., Zhang, A., and Li, R. (2023). Enhancing robustness of gradient-boosted decision trees through one-hot encoding and regularization. arXiv preprint arXiv:2304.13761
2023 arXiv
-
[10]
and Li, R
Fan, J. and Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American statistical Association , 96(456):1348--1360
2001
-
[11]
Friedrich, S., Groll, A., Ickstadt, K., Kneib, T., Pauly, M., Rahnenf \"u hrer, J., and Friede, T. (2023). Regularization approaches in clinical biostatistics: A review of methods and their applications. Statistical Methods in Medical Research , 32(2):425--440
2023
-
[12]
and Satija, R
Hafemeister, C. and Satija, R. (2019). Normalization and variance stabilization of single-cell rna-seq data using regularized negative binomial regression. Genome biology , 20(1):296
2019
-
[13]
and Tibshirani, R
Hastie, T. and Tibshirani, R. (2004). Efficient quadratic regularization for expression arrays. Biostatistics , 5(3):329--340
2004
-
[14]
Hastie, T., Tibshirani, R., Friedman, J., and Franklin, J. (2005). The elements of statistical learning: data mining, inference and prediction. The Mathematical Intelligencer , 27(2):83--85
2005
-
[15]
Hoerl, A. E. and Kennard, R. W. (1970). Ridge regression: Biased estimation for nonorthogonal problems. Technometrics , 12(1):55--67
1970
-
[16]
Kotsilieris, T., Anagnostopoulos, I., and Livieris, I. E. (2022). Regularization techniques for machine learning and their applications
2022
-
[17]
Liang, H., Liu, X., Li, R., and Tsai, C.-L. (2010). Estimation and testing for partially linear single-index models. Annals of statistics , 38(6):3836
2010
-
[18]
Ma, L., Qiao, Y., Wang, R., Chen, H., Liu, G., Xiao, H., and Dai, R. (2024). Machine learning models decoding the association between urinary stone diseases and metabolic urinary profiles. Metabolites , 14(12):674
2024
-
[19]
Ng, A. Y. (2004). Feature selection, l 1 vs. l 2 regularization, and rotational invariance. In Proceedings of the twenty-first international conference on Machine learning , page 78
2004
-
[20]
E., Kim, K.-Y
Scheetz, T. E., Kim, K.-Y. A., Swiderski, R. E., Philp, A. R., Braun, T. A., Knudtson, K. L., Dorrance, A. M., DiBona, G. F., Huang, J., Casavant, T. L., Sheffield, V. C., and Stone, E. M. (2006). Regulation of gene expression in the mammalian eye and its relevance to eye dise...
2006
-
[21]
and Zhang, Y
Tian, Y. and Zhang, Y. (2022). A comprehensive survey on regularization strategies in machine learning. Information Fusion , 80:146--166
2022
-
[22]
Tibshirani, R. (1996). Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series B (Methodological) , 58(1):267--288
1996
-
[23]
M., Kasper, D
Vieira, H. M., Kasper, D. P., Wang, R., Smith, L. M., Enke, C. A., Bergan, R. C., Teply, B. A., and Baine, M. J. (2023). Comparison of sequential versus concurrent chemoradiation regimens in non-metastatic muscle-invasive bladder cancer. Radiation Oncology Journal , 41(3):154
2023
-
[24]
L., Lampe, J
Wang, R., Dai, R., Huang, Y., Neuhouser, M. L., Lampe, J. W., Raftery, D., Tabung, F. K., and Zheng, C. (2023). Variable selection with fdr control for noisy data--an application to screening metabolites that are associated with breast and colorectal cancer. arXiv preprint arX...
2023 arXiv
-
[25]
and Peng, X
Wang, R. and Peng, X. (2022). Normalized iterative hard thresholding in high dimensional logistic regression with metagenomics data. International Journal of Intelligent Technologies & Applied Statistics , 15(1)
2022
-
[26]
and Wang, S
Zhang, B. and Wang, S. (2022). Stochastic functional linear models for gene-based association analysis of quantitative traits in longitudinal studies. Statistics and its interface , 15
2022
-
[27]
Zhang, C.-H. et al. (2010). Nearly unbiased variable selection under minimax concave penalty. The Annals of statistics , 38(2):894--942
2010
-
[28]
and Yu, B
Zhao, P. and Yu, B. (2006). On model selection consistency of lasso. Journal of Machine learning research , 7(Nov):2541--2563
2006
-
[29]
and Hastie, T
Zou, H. and Hastie, T. (2005). Regularization and variable selection via the elastic net. Journal of the royal statistical society: series B (statistical methodology) , 67(2):301--320
2005
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.