Correlated Feature Selection with Extended Exclusive Group Lasso

2020·Arxiv

Abstract

Abstract

In many high dimensional classification or regression problems set in a biological context, the complete identification of the set of informative features is often as important as predictive accuracy, since this can provide mechanistic insight and conceptual understanding. Lasso and related algorithms have been widely used since their sparse solutions naturally identify a set of informative features. However, Lasso performs erratically when features are correlated. This limits the use of such algorithms in biological problems, where features such as genes often work together in pathways, leading to sets of highly correlated features. In this paper, we examine the performance of a Lasso derivative, the exclusive group Lasso, in this setting. We propose fast algorithms to solve the exclusive group Lasso, and introduce a solution to the case when the underlying group structure is unknown. The solution combines stability selection with random group allocation and introduction of artificial features. Experiments with both synthetic and real-world data highlight the advantages of this proposed methodology over Lasso in comprehensive selection of informative features.

1. Introduction

New technologies, especially the introduction of massively parallel Next Generation Sequencing of DNA have generated an explosion of very high dimensional biological and biomedical data sets. Intensively studied examples include global transcriptomics, typically involving the simultaneous measurement of 20–100 thousand RNA transcripts, and adaptive immune repertoires, typically involving millions of different RNA or DNA sequences per sample. Generally, these data sets are very sparse, in the sense that only a small proportion of the measurements (features) are relevant to the biological process being investigated, while the rest contribute only noise. Furthermore, features often show high levels of correlation, since they represent the activity of biological components acting together in groups to perform a task. Neither the correlation structure, nor the classification of features as informative or not, is usually known a priori.

In this context, identifying the complete set of informative features connected to a particular biological process or outcome is an important task in itself, since it can provide fundamental mechanistic insights into how the process works, and is regulated, or misregulated in disease. Lasso (Tibshirani, 1996) offers a promising approach since the sparse solution given by the -norm regularization corresponds naturally to the set of informative features. However, Lasso has some limitations in optimal feature selection. First, Lasso finds at most m features when sample size m is much smaller than feature size n. This is often the case in many real-world problems. Secondly, choosing the regularization parameter is non-trivial, and optimizing classification does not always give the optimal solution in terms of accurate feature selection (Leng et al., 2006; Zhang et al., 2008). Thirdly, Lasso performs erratically when features are correlated, selecting only a few features from each group of correlated features (Zou & Hastie, 2005). Theoretical analysis shows that for Lasso to select correct features, informative features (features with nonzero weights) cannot be overly correlated with irrelevant features (features with zero weights), or other informative features (Zhao & Yu, 2006; Wainwright, 2009; Raskutti et al., 2011). Therefore, Lasso will not discover all informative features when there is a strong correlation structure in the data.

A number of Lasso-type algorithms have therefore been developed (Figure 1), which introduce a group structure into the global feature set. Group Lasso (Yuan & Lin, 2006) forces inter-group sparsity through an -norm regularization so that features from the same group, such as genes in the same biological pathways, are either selected or deselected simultaneously. However, group Lasso requires the group structure to be predefined through prior knowledge, while in many biological problems, such group structure is unknown. In contrast, the exclusive group Lasso (Kong et al., 2014) adopts a complementary approach, introducing an -norm regularization term, which forces intra-group sparsity but relaxes inter-group sparsity. Thus only a few features within each group are selected. If co-correlated informative features are allocated to different groups, the exclusive group Lasso offers a potential solution to the problem of selecting correlated features.

Figure 1. A comparison of the regularization terms of Lasso, group Lasso (GL), and exclusive group Lasso (EGL).

Although the properties of exclusive group Lasso have been investigated in several studies, these have focused either on examples where group structure was known, or where the informative features were uncorrelated, or where an proper group number was required (Kong et al., 2014; Huai et al., 2018; Manabe et al., 2018; Huang & Liu, 2018; Ming et al., 2019; Zhao et al., 2019). In this paper we examine the performance of exclusive group Lasso in the accurate selection of synthetic and real-world correlated features, and introduce new methods for relaxing its limitations.

In Section 2, we examine the limitations of Lasso on synthetic data and real-world immunological problems, which motivates the evolution of exclusive group Lasso in Section 3 to select biologically meaningful features. In this section, we develop efficient algorithms for solving the exclusive group Lasso. Lasso can be solved by efficient coordinate descent algorithms (Wu & Lange, 2008; Friedman et al., 2007; 2010). In an analogous manner, group Lasso can be solved by fast implementation of block coordinate descent (Qin et al., 2013). Other strategies such as active set have also been used in Lasso (Friedman et al., 2010) and group Lasso (Meier et al., 2008; Vogt & Roth, 2012). In this paper we extend the algorithm proposed in (Kong et al., 2014) in a different re-weighting approach and develop it to carry out coordinate descent with active set to achieve a fast solution for exclusive group Lasso. Section 4 proposes random group allocation with stability selection (Meinshausen & B¨uhlmann, 2010) and artificial features for solving the problem of unknown group structure. Synthetic experiments that compare Lasso and our methodology are also included in Section 4. Performance of the proposed methodology on real-world problems is validated in Section 5. We summarize our main findings and contributions in Section 6.

Notation Throughout this paper, we represent matrices as upper case, non-bold letters. Vectors are represented by lower-case, bold letters. We denote by the input labels and by the weight vector. We let X represent an matrix of input data, denote the j-th feature, and denote the i-th sample. The term “relevant features” or “informative features” represents features with nonzero weights such that , and “irrelevant features” represents features with zero weights where We assume the X are standardized by default.

2. Lasso and correlation

In Section 1, we have reviewed that Lasso does not perform well on correlated features. Let us consider a simple experiment to verify this. We generate a dataset such that

where , is Gaussian noise, informative features , and irrelevant features . We split the first features into two groups of 25. Features in the same group are pairwise correlated with correlation . We set the weight vector , where the first 50 entries equal to with random signs. Therefore, the first 50 features are relevant and informative, while the remaining features are irrelevant. Figure 2 shows the false discovery rate (FDR) of Lasso with the largest possible so that all informative features become discovered based on 200 random datasets. As increases, Lasso can only select all informative features at the cost of including irrelevant features.

Figure 2. FDR in a feature selection task as a function of the strength of the correlation among the informative features. Selection threshold was set such that all informative features become chosen. Selection method: Lasso, on 200 random datasets.

To study the performance of Lasso on real-world biological problems, we consider a large dataset (Emerson et al.,

Table 1. Description of the CMV dataset.

Table 2. Classification accuracy on CMV dataset.

2017). The study reported the T-cell receptors (TCRs) from a cohort of healthy individuals, a proportion of which are known to be infected by human cytomegalovirus (CMV), a common virus which does not cause any symptoms in otherwise healthy individuals. A very small proportion of T cells are known to recognize CMV, and control its growth. A fundamental dogma of the current model of adaptive immunity is the clonal selection theory (Burnet, 1957), which proposes that T cells recognizing a virus proliferate, and that CMV-related TCRs are therefore more likely to be detected in individuals infected with CMV. The task which the study therefore addressed was to identify a small set of TCRs which would be over-represented in the infected group (reflecting exposure to the virus) and could be used to predict which individuals were carrying the virus. The dataset contained two independent US cohorts, Cohort 1 and 2, for training and test. To further test the performance, we also use another geographically independent dataset (Cohort 3) that was collected by Belgium researchers (Neuter et al., 2019). Table 1 summarizes the three cohorts.

We used Cohort 1 for training, and Cohorts 2 and 3 for test. To improve the efficiency, we preselected TCRs that occurred in at least 10 samples in Cohort 1. This reduced the number of unique TCRs to 241,008. We applied a 10-fold cross-validation for parameter tuning. As 35 samples were published for 33 patients in Cohort 3, we also included the two additional measurements without affecting the training model. We implemented a two stage strategy, identifying informative TCRs using -regularized SVM (the classification version of Lasso) followed by RBF or linear -regularized SVM for classification. Classification results are shown in Table 2. Among the 1950 TCRs selected by the -SVM as informative features, only 44 were found to be over-represented in individuals infected with CMV.

Clearly, the -regularized algorithm neither selected the biologically meaningful features, nor achieved high accuracy with the selected features. In the meantime, it is not possible to apply group Lasso to this problem as underlying group structure is unknown. This motivates us to apply and adjust exclusive group Lasso to select correlated CMV-related TCRs without knowing the explicit group structure.

3. From Lasso to exclusive group Lasso

The exclusive group Lasso (Kong et al., 2014) is defined as

where f(w) is a convex loss function, G represents group allocation and g represents indices of group . Kong et al. (2014) presented an iterative algorithm to solve the optimization problem. Let us denote by group indicator of the i-th feature. The algorithm proceeds by alternating the following steps.

3.1. A re-weighting scheme

In this section, we provide an alternative solution to the exclusive group Lasso particularly by solving w in (3) via a re-weighting approach. For convenience we rewrite (2) as

which can be further rewritten as

Let and substitute w with; the loss function then becomes . By denoting , (6) can be rewritten as

Therefore, the exclusive group Lasso converts into an optimization problem with re-weighted features and weights , which can make use of the widely-studied -regularized optimization techniques. An additional benefit of the re-weighting step is that it allows us to introduce additional constraints on the values of the weights. For example, biological knowledge would suggest that CMV infection would lead to over-representation of some TCRs, but not necessarily to under-representation of TCRs in the infected group. To capture this biological prior knowledge, we imposed the constraint that weights should be all be of the same sign. We summarize the re-weighting approach with adjustments in Algorithm 1.

3.2. An iterative approach

Exclusive group Lasso with square loss is given by

Luo & Tseng (1992); Tseng (2001) proved that coordinate descent algorithms converge for the problem

under these conditions: 1) g(x) is convex and differentiable; 2) each is convex; and 3) when is a vector, it cannot overlap with other vectors. Let represent the complement of g in the set of features. Then (8) can be solved by iteratively solving subproblems of individual blocks,

where

Solving problem (10) is equivalent to solving

which can be further rewritten as

Now (10) has become converted to Lasso with a new regularization parameter . Theoretically we can solve (13) with every possible using Lasso, and the optimal solution is recovered when the solution of (13) satisfies

, under the current choice of . In practice, solving the problem is fairly fast with bisection approaches, given the piecewise linear property of Lasso (Osborne et al., 2000).

Algorithm 2 shows a bisection algorithm to find the optimal solution to (12) and Algorithm 3 shows solution to (8) with block coordinate descent. It is possible to extend Algorithm 3 to other types of loss functions, as long as the loss function is convex and differentiable and the corresponding -norm solver is piecewise linear.

3.3. A faster solution with the active set of features

For large-scale datasets coordinate descent in -regularized methods such as Lasso is often faster than with -regularized methods like ridge regression. We can therefore reduce the computational time of Algorithm 1 by taking advantage of the sparsity of -regularized algorithms such as Algorithm 3. We first run a single iteration over all groups in Algorithm 3, then apply Algorithm 1 with a reduced feature space of features receiving nonzero weights from the other algorithm until convergence. We usually need another cycle of Algorithm 3 on the complete feature space or to check KKT conditions to ensure convergence—if not converged, we simply run Algorithm 1 again on nonzero features. The double-iteration procedure is repeated until convergence. We summarize the combined algorithm in Algorithm 4.

Table 3. Computational time (s) of exclusive group Lasso.

Optimization in (Kong et al., 2014; Ming et al., 2019) computed for (3), which resulted in 10,221s even when . Therefore we compared the computational time of exclusive group lasso with active set (EGL-AS) and re-weighting scheme (EGL-RW). We randomly generated 1,000 or 10,000 samples with various numbers of features (n) and pairwise correlated informative features () at . We set with random signs. Exclusive group Lasso was run with fixed groups on 10 s using Scikit-learn (Pedregosa et al., 2011) in Python 2.7. The computational times are reported in Table 3. There is a pronounced speed-up in convergence time once n exceeds 5000.

4. Random group allocation and beyond

4.1. Fixed and random group allocation

As group Lasso, exclusive group Lasso also requires reliable group allocation. Figure 3 shows the empirical probabilities that a feature is selected when other informative or irrelevant features coexist in the same group. The probabilities were averaged over 100 random datasets, each containing 50 relevant features and 150 irrelevant features. As expected, informative features are more likely to be discovered when no other informative features coexist in the same group,

Figure 3. Empirical probabilities as a function of strength of regularization (regularization parameter : Prob[a relevant feature is selected | no other relevant features exist in the same group]; Prob[an irrelevant feature is selected | no relevant features exist in the same group]; : Prob[a relevant feature is selected | other relevant features exist in the same group]; : Prob[an irrelevant feature is selected | at least one relevant feature exists in the same group].

while the probability of selecting an irrelevant feature is lower if a group contains at least one informative feature. The ideal group allocation is therefore to allocate each informative (and correlated) feature into a separate group, and set the number of groups to the number of informative features. We refer to this ideal group allocation as “fixed groups”. In most real-world cases informative features are unknown, so groups are allocated randomly (“random groups”).

We compared the performance of Lasso and exclusive group Lasso with fixed and random groups on synthetic datasets that were generated such that , with being the covariance matrix. We inherit the weight vector in Section 2 with the first entries being some scalar with random signs, and the last entries being 0. Noise followed a Gaussian distribution . We computed . To validate the performance on different covariances, four types of were considered.

1. Example 1 (pairwise correlation): relevant features are pairwise correlated at correlation . This simulates the case when active biological patterns such as genes and TCRs are correlated because of functional similarity.

2. Example 2 (random pairwise correlation): relevant features are pairwise correlated at , while some irrelevant features are correlated with randomly selected relevant features at correlation . Similar real-world examples

include discovery of virus-related patterns, which may share lower correlation with other virus mutations, e.g. a certain type of flu virus and its mutations.

3. Example 3 (blocked diagonal correlation): relevant features are grouped into equal-sized groups and features in the same group are pairwise correlated at . If features from the same group are adjacent to each other, all entries in are zero, except for the blocks along the diagonal. This resembles the situation when desired features have underlying subgroup structure such as disease subtypes and active areas in brain scan images.

4. Example 4 (Erdos-Renyi correlation): relevant features are correlated at based on the connectivity of an Erdos-Renyi graph. This is a more complexed case of Example 3, perhaps closer to real-world scenarios.

We set . For Example 2, set to 0.6 because it was the maximum correlation under this setting. We used the same for Examples 1 and 3. In addition, 20 irrelevant features correlated with two randomly selected relevant features at correlation 0.3. For Example 3, we used 5 blocks. For Example 4, because 0.3 was the maximum value to ensure to be positive definite. The connection probability was set to to mimic the structure of Example 3. We set and . The experiments were based on 50 randomly generated datasets. For random exclusive group Lasso, 50 random groups were used. We evaluated the feature selection results with F measure, a common metric in information retrieval. F measure is computed by

where , is the number of relevant features, is the number of selected relevant features, and is the number of all selected features. Table 4 compares the results of Lasso and exclusive Lasso with fixed and random groups.

In most experiments, exclusive group Lasso with fixed groups outperformed Lasso. Exclusive group Lasso with random group allocation also outperformed Lasso, but performed worse than using fixed groups. However, fixed group allocation is usually not possible in real-world scenarios where group structure is unknown for most problems. We therefore introduce stability selection to improve the stability of feature selection with random group allocation.

4.2. Stability selection on random groups

Stability selection combines subsampling with high-dimensional selection algorithms. Meinshausen & B¨uhlmann (2010) showed that stability selection performs

Table 4. F measure of Lasso, exclusive group Lasso with fixed groups (EGL-F) and random groups (EGL-R). EGL outperformed Lasso in almost all experiments.

well even when Lasso fails to select correct features. We adapt stability selection with random group allocation in exclusive group Lasso in Algorithm 5. In practice, a threshold for selection probability is required to distinguish informative and irrelevant features. This can be done with cross-validation with an extra parameter to tune. In our experiments, we rely on k-means to partition probabilities for informative and irrelevant features. We applied stability selection with 50 iterations to the same datasets used in Table 4 and compared Lasso and exclusive group Lasso with fixed and random groups in Table 5. The introduction of stability selection provided moderate improvement of exclusive group Lasso with random groups in almost all cases. So far, we have tested the performance when m = n under various correlation types in the above experiments. We further examine the results when m < n where , and . Table 6 shows the average F measure using stability selection. Exclusive group Lasso with random groups substantially outperformed Lasso in all types of data under these high-dimensional conditions.

Table 5. F measure of Lasso, EGL-F, and EGL-R with stability selection (S). (S) showed improvements and EGL-R(S) performed better than EGL-R particularly when EGL did not perform well.

Table 6. F measure for . Compared to n = 100, Lasso performed worse while EGL had similar F measure.

4.3. Artificial features

Recalling Figure 3, it is clear that when a group contains only irrelevant features, some irrelevant features from this group are selected. However, to ensure that two informative features are rarely found in the same group, especially taking into consideration that the real number of informative features is usually unknown, the number of groups will be chosen as larger than the estimated number of informative features. A significant proportion of the groups may therefore contain no informative feature. Hence, we propose a novel method to further reduce the risk of including too many irrelevant features by inserting some artificially generated informative features into feature groups.

For a regression problem, this approach involves:

1. Randomly generate labels

2. Insert a feature from into each group;

3. Find solution with

4. Remove selected artificial features after convergence.

The procedure for classification is slightly different since the labels need to remain unchanged. We let the artificial features corresponding to each class follow different distributions, so the first step of the above method becomes

1. Randomly generate , where rows of positive samples , and rows of negative samples are diagonal co- variance matrices.

In practice, we find Gaussian distributed and , such as , generally achieve good performance. of fixed values also works well in some experiments.

To evaluate the performance of artificial features, we focused on Examples 3 and 4, where Lasso performed similarly or better, with a slightly more challenging task. We set and . We increased the correlation level for Example 3 to 0.99, which is very challenging to Lasso, and for Example 4 to ensure positive definite covariance. We set and furthermore, we restricted the number of positive and negative weights for informative features to be equal (25 positively weighted and 25 negatively weighted informative features). This is challenging for Lasso since the condition for discovering correct features, the irrepresentable condition (Zhao & Yu, 2006), is easily violated. For Example 3, we used 10 blocks of 5 features. The connection probability in Example 4 was set to to simulate the connection in Example 3. We used 70 groups in random group allocation. Artificial features and weights were generated such that (Example 3) or (Example 4). All experiments were performed with stability selection of 50 iterations on 100 randomly generated datasets. Apart from F measure, we also report the average selection probabilities for informative and irrelevant features in Table 7. Selection probabilities are shown in Figure 4. Lasso performed well on Example 4, perhaps because of the low correlation (). However, Lasso failed to select a large proportion of informative features in Example 3. Exclusive group Lasso performed well on both examples. The introduction of artificial features led to lower selection probabilities (more “compressed” points as shown in Figure 4) for irrelevant features but higher probabilities for informative features. Therefore, EGL-R(SA) is promising in more difficult cases.

5. Real-world application

For real-world application, let us consider the datasets introduced in Section 2. We generated artificial features for CMV+ and CMV- samples separately such that

the artificial features of CMV+ and CMV- samples, respec-

Table 7. F measure (F) and mean selection probabilities of informative (I) and irrelevant (IR) features using Lasso, EGL-F(S), EGL-R(S), EGL-R(SA) (random groups with stability selection and artificial features).

Figure 4. “Stacked” selection probabilities from 100 random datasets. EGL clearly distinguished informative and irrelevant features in both examples while Lasso failed Example 3. Probabilities using EGL-R(SA) were lower for irrelevant features, but higher for informative features than EGL-R(S).

tively. The artificial features were added to 200 random groups, each of which contained a single artificial feature. Since CMV-specific TCRs are likely to be more frequent in the CMV-infected group, we restricted the selected features to those that have positive weights with respect to CMV

Table 8. Feature selection and classification accuracy on CMV data. #TCR: number of selected TCRs, #CMV+TCR: number of selected TCRs that are over-represented in CMV-infected samples.

class. Again, exclusive group Lasso was only used for TCR selection, while the classification task was performed by SVM. Stability selection with 30 iterations was also applied to feature selection. Similar to Section 2, the proposed method was applied to Cohort 1 for training and Cohort 2 and 3 for test. Parameters were selected by a 10-fold cross-validation. Feature selection results and classifica-tion accuracy are reported in Table 8. For feature selection, EGL-R(SA) selects 47 TCRs, all of which are over-represented in CMV+ samples. Compared to -regularized SVM, EGL-R(SA) achieved much higher accuracy with a sparse set of biologically meaningful features. The finding that EGL-R(SA) can select a very small set of features (47 from > 200, 000) which predict CMV status accurately in two independently-generated validation sets suggests the algorithm finds a stable set of informative features. Further biological experiments will be required to validate that these TCRs are indeed specific for CMV.

6. Conclusion

In this paper, we developed fast algorithms to solve exclusive group Lasso, which also allows the possibility to introduce a priori constraints on the values of the weights which may be derived from prior knowledge. Motivated by the challenges of high-dimensional biological datasets, we addressed the problems that informative features are often highly correlated, and neither the number of informative features, nor the underlying correlation structure (i.e. group structure) is usually known. We proposed solutions with stability selection and artificial features for unknown underlying group structure. The exclusive group Lasso, incorporating random group allocation, with stability selection and artificial features, outperformed Lasso and achieved high performance and robust feature selection in a wide range of challenging problems. Future work will be required to further improve model selection and parameter tuning. However, we believe these new feature selection algorithms will prove useful in extracting mechanistic and conceptual understanding from the increasingly complex “omic” datasets being generated across the biomedical domain.

References

Burnet, F. M. A modification of Jerne’s theory of antibody production using the concept of clonal selection. Australian Journal of Science, 1957.

Emerson, R., DeWitt, W., Vignali, M., Gravley, J., Hu, J., Osborne, E., Desmarais, C., Klinger, M., Carlson, C., Hansen, J., Rieder, M., and Robins, H. Immunosequencing identifies signatures of cytomegalovirus exposure history and HLA-mediated effects on the T cell repertoire. Nature Genetics, 49(5):659–665, 2017.

Friedman, J., Hastie, T., Hofling, H., and Tibshirani, R. Pathwise coordinate optimization. The Annals of Applied Statistics, 1(2):302–332, 2007.

Friedman, J., Hastie, T., and Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. Journal of Statistical Software, 33(1):1–21, 2010.

Huai, M., Miao, C., Suo, Q., Li, Y., Gao, J., and Zhang, A. Uncorrelated patient similarity learning. In Proceedings of the 2018 SIAM International Conference on Data Mining, pp. 270–278, 2018.

Huang, Y. and Liu, J. Exclusive sparsity norm minimization with random groups via cone projection. IEEE Transactions on Neural Networks and Learning Systems, 29(12):6145–6153, 2018. ISSN 2162-2388. doi: 10.1109/TNNLS.2018.2819958.

Kong, D., Fujimaki, R., Liu, J., Nie, F., and Ding, C. Exclusive feature learning on arbitrary structures via In Ghahramani, Z., Welling, M., Cortes, C., Lawrence, N. D., and Weinberger, K. Q. (eds.), Advances in Neural Information Processing Systems 27, pp. 1655–1663. Curran Associates, Inc., 2014.

Leng, C., Lin, Y., and Wahba, G. A note on the LASSO and related procedures in model selection. Statistica Sinica, 16:1273–1284, 2006.

Luo, Z. Q. and Tseng, P. On the convergence of the coordinate descent method for convex differentiable minimization. Journal of Optimization Theory and Applications, 72(1):7–35, 1992.

Manabe, H., Hayashi, K., and Shimbo, M. Data-dependent learning of symmetric/antisymmetric relations for knowledge base completion. In Proceedings of the 32th AAAI Conference on Artificial Intelligence, volume 32, 2018.

Meier, L., van de Geer, S., and B¨uhlmann, P. The group lasso for logistic regression. Journal of the Royal Statistical Society. Series B (Statistical Methodology), 70(1):53–71, 2008.

Meinshausen, N. and B¨uhlmann, P. Stability selection. Journal of the Royal Statistical Society. Series B (Statistical Methodology), 72(4):417–473, 2010.

Ming, D., Ding, C., and Nie, F. A probabilistic derivation of lasso and l12-norm feature selections. In Proceedings of the 33rd AAAI Conference on Artificial Intelligence, volume 33, 2019.

Neuter, N. D., Bartholomeus, E., Elias, G., Keersmaekers, N., Suls, A., Jansens, H., Smits, E., Hens, N., Beutels, P., Damme, P. V., Mortier, G., Tendeloo, V. V., Laukens, K., Meysman, P., and Ogunjimi, B. Memory CD4+ T cell receptor repertoire data mining as a tool for identifying cytomegalovirus serostatus. Genes & Immunity, 20(3): 255–260, 2019.

Osborne, M., Presnell, B., and Turlach, B. A new approach to variable selection in least squares problems. IMA Journal of Numerical Analysis, 20(3):389–404, 2000.

Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., and Duchesnay, E. Scikit-learn: Machine Learning in Python . Journal of Machine Learning Research, 12:2825–2830, 2011.

Qin, Z., Scheinberg, K., and Goldfarb, D. Efficient blockcoordinate descent algorithms for the group Lasso. Mathematical Programming Computation, 5(2):143–169, 2013.

Raskutti, G., Wainwright, M., and Yu, B. Minimax rates of estimation for high-dimensional linear regression over balls. IEEE Transactions on Information Theory, 57(10): 6976–6994, 2011.

Tibshirani, R. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series B (Methodological), 58(1):267–288, 1996.

Tseng, P. Convergence of a block coordinate descent method for nondifferentiable minimization. Journal of Optimization Theory and Applications, 109(3):475–494, 2001.

Vogt, J. E. and Roth, V. A complete analysis of the group-Lasso. In Proceedings of the 29th International Conference on Machine Learning, Madison, WI, USA, 2012. Omnipress.

Wainwright, M. Sharp thresholds for high-dimensional and noisy sparsity recovery using -constrained quadratic programming (lasso). IEEE Transactions on Information Theory, 55(5):2183–2202, 2009.

Wu, T. T. and Lange, K. Coordinate descent algorithms for Lasso penalized regression. The Annals of Applied Statistics, 2(1):224–244, 2008.

Yuan, M. and Lin, Y. Model selection and estimation in regression with grouped variables. Journal of the Royal Statistical Society. Series B (Statistical Methodology), 68: 49–67, 2006.

Zhang, J., Jeng, X. J., and Liu, H. Some two-step procedures for variable selection in high-dimensional linear regression. arXiv:0810.1644 [math.ST], 2008.

Zhao, H., Ding, S., Li, X., and Zhao, L. norm independently interpretable regularization based sparse coding for highly correlated data. IEEE Access, 7:53542–53554, 2019. doi: 10.1109/ACCESS.2019.2911004.

Zhao, P. and Yu, B. On model selection consistency of Lasso. Journal of Machine Learning Research, 7:2541– 2563, 2006.

Zou, H. and Hastie, T. Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society. Series B (Statistical Methodology), 67(2):301–320, 2005.