Lasso Screening Rules via Dual Polytope Projection

2012·Arxiv

Abstract

Abstract

Lasso is a widely used regression technique to find sparse representations. When the dimension of the feature space and the number of samples are extremely large, solving the Lasso problem remains challenging. To improve the efficiency of solving large-scale Lasso problems, El Ghaoui and his colleagues have proposed the SAFE rules which are able to quickly identify the inactive predictors, i.e., predictors that have 0 components in the solution vector. Then, the inactive predictors or features can be removed from the optimization problem to reduce its scale. By transforming the standard Lasso to its dual form, it can be shown that the inactive predictors include the set of inactive constraints on the optimal dual solution. In this paper, we propose an efficient and effective screening rule via Dual Polytope Projections (DPP), which is mainly based on the uniqueness and nonexpansiveness of the optimal dual solution due to the fact that the feasible set in the dual space is a convex and closed polytope. Moreover, we show that our screening rule can be extended to identify inactive groups in group Lasso. To the best of our knowledge, there is currently no exact screening rule for group Lasso. We have evaluated our screening rule using synthetic and real data sets. Results show that our rule is more effective in identifying inactive predictors than existing state-of-the-art screening rules for Lasso.

1 Introduction

Data with various structures and scales comes from almost every aspect of daily life. To effectively extract patterns in the data and build interpretable models with high prediction accuracy is always desirable. One popular technique to identify important explanatory features is by sparse regularization. For instance, consider the widely used -regularized least squares regression problem known as Lasso [31]. The most appealing property of Lasso is the sparsity of the solutions, which is equivalent to feature selection. Suppose we have N observations and p features. Let y denote the N dimensional response vector and feature matrix. Let be the regularization parameter. The Lasso problem is formulated as the following optimization problem:

Lasso has achieved great success in a wide range of applications [13, 12, 38, 9, 34] and in recent years many algorithms have been developed to efficiently solve the Lasso problem [15, 20, 26, 14, 18, 4, 19]. However, when the dimension of feature space and the number of samples are very large, solving the Lasso problem remains challenging because we may not even be able to load the data matrix into main memory. The idea of screening has been shown very promising in solving Lasso for large-scale problems. Essentially, screening aims to quickly identify the inactive features that have 0 components in the solution and then remove them from the optimization. Therefore, we can work on a reduced feature matrix to solve the Lasso problem, which may lead to substantial savings in computational cost and memory usage.

Existing screening methods for Lasso can be roughly divided into two categories: the Heuristic Screening Methods and the Safe Screening Methods. As the name indicated, the heuristic screening methods can not guarantee that the discarded features have zero coefficients in the solution vector. In other words, they may mistakenly discard the active features which have nonzero coefficients in the sparse representations. Well-known heuristic screening methods for Lasso include SIS [17] and strong rules [32]. SIS is based on the associations between features and the prediction task, but not from an optimization point of view. Strong rules rely on the assumption that the absolute values of the inner products between features and the residue are nonexpansive [3] with respect to the parameter values. Notice that, in real applications, this assumption is not always true. In order to ensure the correctness of the solutions, strong rules check the KKT conditions for violations. In case of violations, they weaken the screened set and repeat this process. In contrast to the heuristic screening methods, the safe screening methods for Lasso can guarantee that the discarded features are absent from the resulting sparse models. Existing safe screening methods for Lasso includes SAFE [16] and DOME [36, 35], which are based on an estimation of the dual optimal solution. The key challenge of searching for effective safe screening rules is how to accurately estimate the dual optimal solution. The more accurate the estimation is, the more effective the resulting screening rule is in discarding the inactive features. Moreover, Xiang et al. [36] have shown that the SAFE rule for Lasso can be read as a special case of their testing rules.

In this paper, we develop novel efficient and effective screening rules for the Lasso problem; our screening rules are safe in the sense that no active features will be discarded. As the name indicated (DPP), the proposed approaches heavily rely on the geometric properties of the Lasso problem. Indeed, the dual problem of problem (1) can be formulated as a projection problem. More specifically, the dual optimal solution of the Lasso problem is the projection of the scaled response vector onto a nonempty closed and convex polytope (the feasible set of the dual problem). This nice property provides us many elegant approaches to accurately estimate the dual optimal solutions, e.g., nonexpansiveness, firmly nonexpansiveness [3]. In fact, the estimation of the dual optimal solution in DPP is a direct application of the nonexpansiveness of the projection operators. Moreover, by further exploiting the properties of the projection operators, we can significantly improve the estimation of the dual optimal solution. Based on this estimation, we develop the so called enhanced DPP (EDPP) rules which are able to detect far more inactive features than DPP. Therefore, the speedup gained by EDPP is much higher than the one by DPP.

In real applications, the optimal parameter value of is generally unknown and needs to be estimated. To determine an appropriate value of , commonly used approaches such as cross validation and stability selection involve solving the Lasso problems over a grid of tuning parameters . Thus, the process can be very time consuming. To address this challenge, we develop the sequential version of the DPP families. Briefly speaking, for the Lasso problem, suppose we are given the solution . We then apply the screening rules to identify the inactive features of problem (1) at by making use of ). The idea of the sequential screening rules is proposed by [16] and [32] and has been shown to be very effective for the aforementioned scenario. In [32], the authors demonstrate that the sequential strong rules are very effective in discarding inactive features especially for very small parameter values and achieve the state-of-the-art performance. However, in contrast to the recursive SAFE (the sequential version of SAFE rules) and the sequential version of DPP rules, it is worthwhile to mention that the sequential strong rules may mistakenly discard active features because they are heuristic methods. Moreover, it is worthwhile to mention that, for the existing screening rules including SAFE and strong rules, the basic versions are usually special cases of their sequential versions, and the same applies to our DPP and EDPP rules. For the DOME rule [36, 35], it is unclear whether its sequential version exists.

The rest of this paper is organized as follows. We present the family of DPP screening rules, i.e., DPP and EDPP, in detail for the Lasso problem in Section 2. Section 3 extends the idea of DPP screening rules to identify inactive groups in group Lasso [37]. We evaluate our screening rules on synthetic and real data sets from many different applications. In Section 4, the experimental results demonstrate that our rules are more effective in discarding inactive features than existing state-of-the-art screening rules. We show that the efficiency of the solver can be improved by several orders of magnitude with the enhanced DPP rules, especially for the high-dimensional data sets (notice that, the screening methods can be integrated with any existing solvers for the Lasso problem). Some concluding remarks are given in Section 5.

2 Screening Rules for Lasso via Dual Polytope Projections

In this section, we present the details of the proposed DPP and EDPP screening rules for the Lasso problem. We first review some basics of the dual problem of Lasso including its geometric properties in Section 2.1; we also briefly discuss some basic guidelines for developing safe screening rules for Lasso. Based on the geometric properties discussed in Section 2.1, we then develop the basic DPP screening rule in Section 2.2. As a straightforward extension in dealing with the model selection problems, we also develop the sequential version of DPP rules. In Section 2.3, by exploiting more geometric properties of the dual problem of Lasso, we further improve the DPP rules by developing the so called enhanced DPP (EDPP) rules. The EDPP screening rules significantly outperform DPP rules in identifying the inactive features for the Lasso problem.

2.1 Basics

Different from [36, 35], we do not assume have unit length.The dual problem of problem (1) takes the form of (to make the paper self-contained, we provide the detailed derivation of the dual form in the appendix):

where is the dual variable. For notational convenience, let the optimal solution of problem (2) be ) [recall that the optimal solution of problem (1) with parameter is denoted by Then, the KKT conditions are given by:

where [denotes the component. In view of the KKT condition in (4), we have

In other words, we can potentially make use of (R1) to identify the inactive features for the Lasso problem. However, since ) is generally unknown, we can not directly apply (R1) to identify the inactive features. Inspired by the SAFE rules [16], we can first estimate a region which contains ). Then, (R1) can be relaxed as follows:

Clearly, as long as we can find a region which contains ), (R1’) will lead to a screening rule to detect the inactive features for the Lasso problem. Moreover, in view of (R1) and (R1’), we can see that the smaller the region is, the more accurate the estimation of ) is. As a result, more inactive features can be identified by the resulting screening rules.

(2), we can observe that the dual optimal solution is the feasible point which is closest to notational convenience, let the feasible set of problem (2) be F. Clearly, F is the intersection of 2p closed half-spaces, and thus a closed and convex polytope. (Notice that, F is also nonempty since 0 .) In other words, ) is the projection of onto the polytope F. Mathematically, for an arbitrary vector w and a convex set C in a Hilbert space H, let us define the projection operator as

Then, the dual optimal solution ) can be expressed by

Indeed, the nice property of problem (2) illustrated by Eq. (6) leads to many interesting results. For example, it is easy to see that would be an is large enough. If this is the case, we immediately have the following assertions: 1) is an interior point of F implies that none of the constraints of problem (2) would be all ) is an interior point of F as well since by Eq. (6) and the fact . Combining the results in 1) and 2), it is easy to see that all i = 1, . . . , p. By (R1), we can conclude that ) = 0, under the assumption that enough.

The above analysis may naturally lead to a question: does there exist a specific parameter value such that the optimal solution of problem (1) is 0 whenever ? The answer is affirmative. Indeed, let us define

It is well known [32] that defined by Eq. (7) is the smallest parameter such that problem (1) has a trivial solution, i.e.,

Combining the results in (8) and Eq. (3), we immediately have

Therefore, through out the rest of this paper, we will focus on the cases with

In the subsequent sections, we will follow (R1’) to develop our screening rules. More specifically, the derivation of the proposed screening rules can be divided into the following three steps:

Step 1. We first estimate a region which contains the dual optimal solution

Step 2. We solve the maximization problem in (R1’), i.e., sup

Step 3. By plugging in the upper bound we find in Step 2, it is straightforward to develop the screening rule based on (R1’).

The geometric property of the dual problem illustrated by Eq. (6) serves as a fundamentally important role in developing our DPP and EDPP screening rules.

2.2 Fundamental Screening Rules via Dual Polytope Projections (DPP)

In this Section, we propose the so called DPP screening rules for discarding the inactive features for Lasso. As the name indicates, the idea of DPP heavily relies on the properties of projection operators, e.g., the nonexpansiveness [5]. We will follow the three steps stated in Section 2.1 to develop the DPP screening rules.

First, we need to find a region which contains the dual optimal solution ). Indeed, the result in (9) provides us an important clue. That is, we may be able to estimate a possible region for ) in terms of a known . Notice that, we can always set make use of the fact that implied by (9). Another key ingredient comes from Eq. (6), i.e., the dual optimal solution ) is the projection of onto the feasible set F, which is nonempty closed and convex. A nice property of the projection operators defined in a Hilbert space with respect to a nonempty closed and convex set is the so called nonexpansiveness. For convenience, we restate the definition of nonexpansiveness in the following theorem.

Theorem 1. Let C be a nonempty closed convex subset of a Hilbert space H. Then the projection operator defined in Eq. (5) is continuous and nonexpansive, i.e.,

In view of Eq. (6), a direct application of Theorem 1 leads to the following result:

Theorem 2. For the Lasso problem, let be two regularization parameters. Then,

2 actually implies that the dual optimal solution must be inside a ball centered at ) with radius

We thus complete the first step for developing DPP. Because it is easy to find the upper bound of a linear functional over a ball, we combine the remaining two steps as follows.

Theorem 3. For the Lasso problem, assume we are given the solution of its dual problem for a specific be a positive value different from

Proof. The dual optimal solution ) is estimated to be inside the ball given by Eq. (12). To simplify notations, let . To develop a screening rule based on (R1’), we need to solve the optimization problem: sup

Indeed, for any ), it can be expressed by:

Therefore, the optimization problem can be easily solved as follows:

By plugging the upper bound in Eq. (14) to (R1’), we obtain the statement in Theorem 3, which completes the proof.

Theorem 3 implies that we can develop applicable screening rules for Lasso as long as the dual optimal solution ) is known for a certain parameter value . By simply setting and noting that [please refer to Eq. (9)], Theorem 3 immediately leads to the following result.

Corollary 4. Basic DPP: For the Lasso problem (1), let [. Otherwise, [

Remark 1. Notice that, DPP is not the same as ST1 [36] and SAFE [16], which discards the

feature if

From the perspective of the sphere test, the radius of ST1/SAFE and DPP are the same. But the centers of ST1 and DPP are respectively, which leads to different formulas, i.e., Eq. (15) and Corollary 4.

In real applications, the optimal parameter value of is generally unknown and needs to be estimated. To determine an appropriate value of , commonly used approaches such as cross validation and stability selection involve solving the Lasso problem over a grid of tuning parameters , which is very time consuming. Motivated by the ideas of [32] and [16], we develop a sequential version of DPP rules. We first apply the DPP screening rule in Corollary 4 to discard inactive features for the Lasso problem (1) with parameter being . After solving the reduced optimization problem for , we obtain the exact solution ). Hence by Eq. (3), we can find ). According to Theorem 3, once we know the optimal dual solution construct a new screening rule by setting to identify inactive features for problem (1) with parameter being . By repeating the above process, we obtain the sequential version of the DPP rule as in the following corollary.

Corollary 5. Sequential DPP: For the Lasso problem (1), suppose we are given a sequence of parameter values . Then for any integer 0 [is known and the following holds:

Remark 2. From Corollaries 4 and 5, we can see that both of the DPP and sequential DPP rules discard the inactive features for the Lasso problem with a smaller parameter value by assuming a known dual optimal solution at a larger parameter value. This is in fact a standard way to construct screening rules for Lasso [32, 16, 36, 35].

Remark 3. For illustration purpose, we present both the basic and sequential version of the DPP screening rules. However, it is easy to see that the basic DPP rule can be easily derived from its sequential version by simply setting . Therefore, in this paper, we will focus on the development and evaluation of the sequential version of the proposed screening rules. To avoid any confusions, DPP and EDPP all refer to the corresponding sequential versions.

2.3 Enhanced DPP Rules for Lasso

In this section, we further improve the DPP rules presented in Section 2.2 by a more careful analysis of the projection operators. Indeed, from the three steps by which we develop the DPP rules, we can see that the first step is a key. In other words, the estimation of the dual optimal solution serves as a fundamentally important role in developing the DPP rules. Moreover, (R1’) implies that the more accurate the estimation is, the more effective the resulting screening rule is in discarding the inactive features. The estimation of the dual optimal solution in DPP rules is in fact a direct consequence of the nonexpansiveness of the projection operators. Therefore, in order to improve the performance of the DPP rules in discarding the inactive features, we propose two different approaches to find more accurate estimations of the dual optimal solution. These two approaches are presented in detail in Sections 2.3.1 and 2.3.2 respectively. By combining the ideas of these two approaches, we can further improve the estimation of the dual optimal solution. Based on this estimation, we develop the enhanced DPP rules (EDPP) in Section 2.3.3. Again, we will follow the three steps in Section 2.1 to develop the proposed screening rules.

2.3.1 Improving the DPP rules via Projections of Rays

In the DPP screening rules, the dual optimal solution ) is estimated to be inside the ball ) given. In this section, we show that ) lies inside a ball centered at ) with a smaller radius.

Indeed, it is well known that the projection of an arbitrary point onto a nonempty closed convex set C in a Hilbert space H always exists and is unique [3]. However, the converse is not true, i.e., there may exist ). In fact, it is known that the following result holds:

Lemma 6. [3] Let C be a nonempty closed convex subset of a Hilbert space H. For a point let . Then, the projection of the point i.e.,

Clearly, when 0 is the ray starting from pointing in the same direction as ). By Lemma 6, we know that the projection of the ray 0 onto the set C is a single point with 0 becomes a single point and the statement in Lemma 6 is trivial.]

By making use of Lemma 6 and the nonexpansiveness of the projection operators, we can improve the estimation of the dual optimal solution in DPP [please refer to Theorem 2 and Eq. (12)]. More specifically, we have the following result:

Theorem 7. For the Lasso problem, suppose the dual optimal solution known. For any , let us define

Then, the dual optimal solution can be estimated as follows:

Proof. By making use of Lemma 6, we present the proof of the statement for the cases with ). We postpone the proof of the statement for the case with introduce more general technical results.

In view of the assumption ), it is easy to see that

For each ), let us define

By the result in (21), we can see that ) defined by Eq. (22) is a ray which starts at points in the same direction as ). In view of Eq. (6), a direct application of Lemma 6 leads to that:

By applying Theorem 1 again, we have

Because the inequality in (24) holds for all 0, it is easy to see that

The inequality in (25) implies that, to prove the first half of the statement, i.e.,

Therefore, we can see that:

The last inequality in (28) is due to the result in (27).

Clearly, in view of (25) and (28), we can see that the first half of the statement holds, i.e., ). The second half of the statement, i.e., ), can be easily obtained by noting that the inequality in (24) reduces to the one in (12) when t = 1. This completes the proof of the statement with

Before we present the proof of Theorem 7 for the case with , let us briefly review some technical results from convex analysis first.

Definition 8. [28] Let C be a nonempty closed convex subset of a Hilbert space The set

is called the normal cone to C at w.

In terms of the normal cones, the following theorem provides an elegant and useful characterization of the projections onto nonempty closed convex subsets of a Hilbert space.

Theorem 9. [3] Let C be a nonempty closed convex subset of a Hilbert space H. Then, for every is the projection of w onto C if and only if

In view of the proof of Theorem 7, we can see that Eq. (23) is a key step. When similar to Eq. (22), let us define

By Theorem 9, the following lemma shows that Eq. (23) also holds for

Lemma 10. For the Lasso problem, let be given by Eq. (17) and Eq. (31), then the following result holds:

Proof. To prove the statement, Theorem 9 implies that we only need to show:

Recall that [Eq. (17)], and [Eq. (9)]. It is easy to see that

Moreover, assume is an arbitrary point of F. Then, we have 1, and thus

Therefore, the inequality in (33) easily follows by combing the results in (34) and (35), which completes the proof.

We are now ready to give the proof of Theorem 7 for the case with

Proof. In view of Theorem 1 and Lemma 10, we have

Because the inequality in (36) holds for all 0, we can see that

Clearly, we only need to show that

By noting that 0

Moreover, because ), it is easy to see that

Therefore, in view of (37) and (39), we can see that the first half of the statement holds, i.e., ). The second half of the statement, i.e.,

can be easily obtained by noting that the inequality in (37) reduces to the one in (12) when t = 0. This completes the proof of the statement with Thus, the proof of Theorem 7 is completed.

Theorem 7 in fact provides a more accurate estimation of the dual optimal solution than the one in DPP, i.e., ) lies inside a ball centered at ) with a radius . Based on this improved estimation and (R1’), we can develop the following screening rule to discard the inactive features for Lasso.

Theorem 11. For the Lasso problem, assume the dual optimal solution known. Then, for each

We omit the proof of Theorem 11 since it is very similar to the one of Theorem 3. By Theorem 11, we can easily develop the following sequential screening rule.

values . Then for any integer 0 is known and the following holds:

The screening rule in Improvement 1 is developed based on (R1’) and the estimation of the dual optimal solution in Theorem 7, which is more accurate than the one in DPP. Therefore, in view of (R1’), the screening rule in Improvement 1 are more effective in discarding the inactive features than the DPP rule.

2.3.2 Improving the DPP rules via Firmly Nonexpansiveness

In Section 2.3.1, we improve the estimation of the dual optimal solution in DPP by making use of the projections of properly chosen rays. (R1’) implies that the resulting screening rule stated in Improvement 1 is more effective in discarding the inactive features than DPP. In this Section, we present another approach to improve the estimation of the dual optimal solution in DPP by making use of the so called firmly nonexpansiveness of the projections onto nonempty closed convex subset of a Hilbert space.

Theorem 12. [3] Let C be a nonempty closed convex subset of a Hilbert space H. Then the projection operator defined in Eq. (5) is continuous and firmly nonexpansive. In other words, for any

where Id is the identity operator.

In view of the inequalities in (40) and (10), it is easy to see that firmly nonexpansiveness implies nonexpansiveness. But the converse is not true. Therefore, firmly nonexpansiveness of the projection operators is a stronger property than the nonexpansiveness. A direct application of Theorem 12 leads to the following result.

Theorem 13. For the Lasso problem, let be two parameter values. Then

Proof. In view of Eq. (6) and the firmly nonexpansiveness in (40), we have

which completes the proof of the first half of the statement. The second half of the statement is trivial by noting that the first inequality in (42) (firmly nonexpansiveness) implies the inequality in (11) (nonexpansiveness) but not vice versa. Indeed, it is easy to see that the ball in the middle of (41) is inside the right one and has only a half radius.

Clearly, Theorem 13 provides a more accurate estimation of the dual optimal solution than the one in DPP, i.e., the dual optimal solution must be inside a ball which is a subset of the one in DPP and has only a half radius. Again, based on the estimation in Theorem 13 and (R1’), we have the following result.

Theorem 14. For the Lasso problem, assume the dual optimal solution known. Then, for each

We omit the proof of Theorem 14 since it is very similar to the proof of Theorem 3. A direct application of Theorem 14 leads to the following sequential screening rule.

values . Then for any integer 0 is known and the following holds:

Because the screening rule in Improvement 2 is developed based on (R1’) and the estimation in Theorem 13, it is easy to see that Improvement 2 is more effective in discarding the inactive features than DPP.

2.3.3 The Proposed Enhanced DPP Rules

In Sections 2.3.1 and 2.3.2, we present two different approaches to improve the estimation of the dual optimal solution in DPP. In view of (R1’), we can see that the resulting screening rules, i.e., Improvements 1 and 2, are more effective in discarding the inactive features than DPP. In this section, we give a more accurate estimation of the dual optimal solution than the ones in Theorems 7 and 13 by combining the aforementioned two approaches together. The resulting screening rule for Lasso is the so called enhanced DPP rule (EDPP). Again, (R1’) implies that EDPP is more effective in discarding the inactive features than the screening rules in Improvements 1 and 2. We also present several experiments to demonstrate that EDPP is able to identify more inactive features than the screening rules in Improvements 1 and 2. Therefore, in the subsequent sections, we will focus on the generalizations and evaluations of EDPP.

To develop the EDPP rules, we still follow the three steps in Section 2.1. Indeed, by combining the two approaches proposed in Sections 2.3.1 and 2.3.2, we can further improve the estimation of the dual optimal solution in the following theorem.

Theorem 15. For the Lasso problem, suppose the dual optimal solution known, and be given by Eq. (19). Then, we have

Proof. Recall that ) is defined by Eq. (22) and Eq. (31). In view of (40), we have

By expanding the second term on the left hand side of (44) and rearranging the terms, we obtain the following equivalent form:

In view of Eq. (6), Eq. (23) and Eq. (32), the inequality in (45) can be rewritten as

[Recall that ) are defined by Eq. (17) and Eq. (18) respectively.] Clearly, the inequality in (46) is equivalent to

The statement follows easily by minimizing the right hand side of the inequality in (47), which has been done in the proof of Theorem 7.

Indeed, Theorem 15 is equivalent to bounding ) in a ball as follows:

Based on this estimation and (R1’), we immediately have the following result.

Theorem 16. For the Lasso problem, assume the dual optimal problem known, and if the following holds:

We omit the proof of Theorem 16 since it is very similar to the one of Theorem 3. Based on Theorem 16, we can develop the EDPP rules as follows.

Corollary 17. EDPP: For the Lasso problem, suppose we are given a sequence of parameter values . Then for any integer 0 is known and the following holds:

It is easy to see that the ball in (48) has the smallest radius compared to the ones in Theorems 7 and 13, and thus it provides the most accurate estimation of the dual optimal solution. According to (R1’), EDPP is more effective in discarding the inactive features than DPP, Improvements 1 and 2.

Figure 1: Comparison of the family of DPP rules on three real data sets: Prostate Cancer digit data set (left), PIE data set (middle) and MNIST image data set (right). The first row shows the rejection ratios of DPP, Improvement 1, Improvement 2 and EDPP. The second row presents the speedup gained by these four methods.

Comparisons of the Family of DPP rules We evaluate the performance of the family of DPP screening rules, i.e., DPP, Improvement 1, Improvement 2 and EDPP, on three real data sets: a) the Prostate Cancer [27]; b) the PIE face image data set [30]; c) the MNIST handwritten digit data set [21]. To measure the performance of the screening rules, we compute the following two quantities:

1. the rejection ratio, i.e., the ratio of the number of features discarded by screening rules to the actual number of zero features in the ground truth;

2. the speedup, i.e., the ratio of the running time of the solver with screening rules to the running time of the solver without screening.

Table 1: Running time (in seconds) for solving the Lasso problems along a sequence of 100 tuning parameter values equally spaced on the scale of 05 to 1 by (a): the solver [22] (reported in the second column) without screening; (b): the solver combined with different screening methods (reported in the 3columns). The last four columns report the total running time (in seconds) for the screening methods.

For each data set, we run the solver with or without the screening rules to solve the Lasso problem along a sequence of 100 parameter values equally spaced on the scale from 0.05 to 1.0. Fig. 1 presents the rejection ratios and speedup by the family of DPP screening rules. Table 1 reports the running time of the solver with or without the screening rules for solving the 100 Lasso problems, as well as the time for running the screening rules.

The Prostate Cancer Data Set The Prostate Cancer data set [27] is obtained by protein mass spectrometry. The features are indexed by time-of-flight values, which are related to the mass over charge ratios of the constituent proteins in the blood. The data set has 15154 measurements of 132 patients. 69 of the patients have prostate cancer and the rest are healthy. Therefore, the data matrix X is of size 132 15154, and the response vector contains the binary labels of the patients.

The PIE Face Image Data Set The PIE face image data set used in this experimentcontains 11554 gray face images of 68 people, taken under different poses, illumination conditions and expressions. Each of the images has 32 32 pixels. Therefore, in each trial, we first randomly pick an image as the response , and then use the remaining images to form the data matrix . We run 100 trials and report the average performance of the screening rules.

The MNIST Handwritten Digit Data Set This data set contains grey images of scanned handwritten digits, including 60, 000 for training and 10, 000 for testing. The dimension of each image is 28 28. We first randomly select 5000 images for each digit from the training set (and in total we have 50000 images) and get a data matrix . Then in each trial, we randomly select an image from the testing set as the response . We run 100 trials and report the average performance of the screening rules.

From Fig. 1, we can see that both Improvements 1 and 2 are able to discard more inactive features than DPP, and thus lead to a higher speedup. Compared to Improvement 2, we can also observe that Improvement 1 is more effective in discarding the inactive features. For the three data sets, the second row of Fig. 1 shows that Improvement 1 leads to about 20, 60, 70 times speedup respectively, which are much higher than the ones gained by Improvement 1 (roughly 10 times for all the three cases).

Moreover, the EDPP rule, which combines the ideas of both Improvements 1 and 2, is even more effective in discarding the inactive features than Improvement 1. We can see that, for all of the three data sets and most of the 100 parameter values, the rejection ratios of EDPP are very close to 100%. In other words, EDPP is able to discard almost all of the inactive features. Thus, the resulting speedup of EDPP is significantly better than the ones gained by the other three DPP rules. For the PIE and MNIST data sets, we can see that the speedup gained EDPP is about 150 and 230 times, which are two orders of magnitude. In view of Table 1, for the MNIST data set, the solver without screening needs about 2566.26 seconds to solve the 100 Lasso problems. In contrast, the solver with EDPP only needs 11.12 seconds, leading to substantial savings in the computational cost. Moreover, from the last four columns of Table 1, we can also observe that the computational cost of the family of DPP rules are very low. Compared to that of the solver without screening, the computational cost of the family of DPP rules is negligible.

In Section 4, we will only compare the performance of EDPP against several other state-of-the-art screening rules.

3 Extensions to Group Lasso

To demonstrate the flexibility of the family of DPP rules, we extend the idea of EDPP to the group Lasso problem [37] in this section. Although the Lasso and group Lasso problems are very different from each other, we will see that their dual problems share a lot of similarities. For example, both of the dual problems can be formulated as looking for projections onto nonempty closed convex subsets of a Hilbert space. Recall that, the EDPP rule for the Lasso problem is entirely based on the properties of the projection operators. Therefore, the framework of the EDPP screening rule we developed for Lasso is also applicable for the group Lasso problem. In Section 3.1, we briefly review some basics of the group Lasso problem and explore the geometric properties of its dual problem. In Section 3.2, we develop the EDPP rule for the group Lasso problem.

3.1 Basics

With the group information available, the group Lasso problem takes the form of:

where is the data matrix for the . The dual problem of (50) is (see detailed derivation in the appendix):

The KKT conditions are given by

for g = 1, 2, . . . , G. Clearly, in view of Eq. (53), we can see that

However, since ) is generally unknown, (R2) is not applicable to identify the inactive groups, i.e., the groups which have 0 coefficients in the solution vector, for the group Lasso problem. Therefore, similar to the Lasso problem, we can first find a region which contains ), and then (R2) can be relaxed as follows:

Therefore, to develop screening rules for the group Lasso problem, we only need to estimate the region which contains ), solve the maximization problem in (R2), and plug it into (R2other words, the three steps proposed in Section 2.1 can also be applied to develop screening rules for the group Lasso problem. Moreover, (R2) also implies that the smaller the region more accurate the estimation of the dual optimal solution is. As a result, the more effective the resulting screening rule is in discarding the inactive features.

Geometric Interpretations For notational convenience, let F be the feasible set of problem (51). Similar to the case of Lasso, problem (51)implies that the dual optimal ) is the projection of onto the feasible set F, i.e.,

Compared to Eq. (6), the only difference in Eq. (54) is that the feasible set F is the intersection of a set of ellipsoids, and thus not a polytope. However, similar to F, F is also a nonempty closed and convex (notice that 0 is a feasible point). Therefore, we can make use of all the aforementioned properties of the projection operators, e.g., Lemmas 6 and 10, Theorems 9 and 12, to develop screening rules for the group Lasso problem.

Moreover, similar to the case of Lasso, we also have a specific parameter value [32] for the group Lasso problem, i.e.,

Indeed, is the smallest parameter value such that the optimal solution of problem (50) is 0. More specifically, we have:

Combining the result in (56) and Eq. (52), we immediately have

Therefore, all through the subsequent sections, we will focus on the cases with

3.2 Enhanced DPP rule for Group Lasso

In view of (R2), we can see that the estimation of the dual optimal solution is the key step to develop a screening rule for the group Lasso problem. Because ) is the projection of the nonempty closed convex set F [please refer to Eq. (54)], we can make use of all the properties of projection operators, e.g., Lemmas 6 and 10, Theorems 9 and 12, to estimate the dual optimal solution. First, let us develop a useful technical result as follows.

Lemma 18. For the group Lasso problem, let be given by Eq. (55) and

Suppose the dual optimal solution is known at , let us define

Then, we have the following result holds

Proof. Let us first consider the cases with ). In view of the definition of easy to see that . Therefore, in view of Eq. (54) and Lemma 6, the statement in Eq. (61) follows immediately.

We next consider the case with . By Theorem 9, we only need to check if

Indeed, in view of Eq. (55) and Eq. (57), we can see that

On the other hand, by Eq. (55) and Eq. (58), we can see that

where is the number of columns of . By plugging Eq. (64) into Eq. (63), we have

Moreover, for any feasible point , we can see that

In view of the result in (66) and Eq. (64), it is easy to see that

Combining the result in Eq. (63) and (67), it is easy to see that the inequality in (62) holds for all , which completes the proof.

By Lemma 18, we can accurately estimate the dual optimal solution of the group Lasso problem in the following theorem. It is easy to see that the result in Theorem 19 is very similar to the one in Theorem 15 for the Lasso problem.

Theorem 19. For the group Lasso problem, suppose the dual optimal solution is known, and is given by Eq. (59). For any , let us define

Then, the dual optimal solution can be estimated as follows:

We omit the proof of Theorem 19 since it is exactly the same as the one of Theorem 15. Indeed, Theorem 19 is equivalent to estimating ) in a ball as follows:

Based on this estimation and (R2), we immediately have the following result.

Theorem 20. For the group Lasso problem, assume the dual optimal solution is known at if the following holds

To simplify notations, let

It is easy to see that

which completes the proof. The second and third inequalities in (73) are due to (72) and Theorem 19, respectively.

In view of Eq. (52) and Theorem 20, we can derive the EDPP rule to discard the inactive groups for the group Lasso problem as follows.

Corollary 21. EDPP: For the group Lasso problem (50), suppose we are given a sequence of parameter values . For any integer 0 is known and the following holds:

4 Experiments

In this section, we evaluate the proposed EDPP rules for Lasso and group Lasso on both synthetic and real data sets. To measure the performance of our screening rules, we compute the rejection ratio and speedup (please refer to Section 2.3.3 for details). Because the EDPP rule is safe, i.e., no active features/groups will be mistakenly discarded, the rejection ratio will be less than one.

In Section 4.1, we conduct two sets of experiments to compare the performance of EDPP against several state-of-the-art screening methods. We first compare the performance of the basic versions of EDPP, DOME, SAFE, and strong rule. Then, we focus on the sequential versions of EDPP, SAFE, and strong rule. Notice that, SAFE and EDPP are safe. However, strong rule may mistakenly discard features with nonzero coefficients in the solution. Although DOME is also safe for the Lasso problem, it is unclear if there exists a sequential version of DOME. Recall that, real applications usually favor the sequential screening rules because we need to solve a sequence of of Lasso problems to determine an appropriate parameter value [32]. Moreover, DOME assumes special structure on the data, i.e., each feature and the response vector should be normalized to have unit length.

In Section 4.2, we compare EDPP with strong rule for the group Lasso problem on synthetic data sets. We are not aware of any safe screening rules for the group Lasso problem at this point. For SAFE and Dome, it is not straightforward to extend them to the group Lasso problem.

4.1 EDPP for the Lasso Problem

For the Lasso problem, we first compare the performance of the basic versions of EDPP, DOME, SAFE and strong rule in Section 4.1.1. Then, we compare the performance of the sequential versions of EDPP, SAFE and strong rule in Section 4.1.2.

4.1.1 Evaluation of the Basic EDPP Rule

In this section, we perform experiments on six real data sets to compare the performance of the basic versions of SAFE, DOME, strong rule and EDPP. Briefly speaking, suppose that we are given a parameter value . Basic versions of the aforementioned screening rules always make use of to identify the zero components of ). Take EDPP for example. The basic version of EDPP can be obtained by replacing ), respectively, in (49) for all k = 1, . . . , K.

In this experiment, we report the rejection ratios of the basic SAFE, DOME, strong rule and EDPP along a sequence of 100 parameter values equally spaced on the scale from 0.05 to 1.0. We note that DOME requires that all features of the data sets have unit length. Therefore, to compare the performance of DOME with SAFE, strong rule and EDPP, we normalize the features of all the data sets used in this section. However, it is worthwhile to mention that SAFE, strong rule and EDPP do not assume any specific structures on the data set. The data sets used in this section are listed as follows:

Figure 2: Comparison of basic versions of SAFE, DOME, Strong Rule and EDPP on six real data sets.

The Colon Cancer Data Set This data set contains gene expression information of 22 normal tissues and 40 colon cancer tissues, and each has 2000 gene expression values.

The Lung Cancer Data Set This data set contains gene expression information of 186 lung tumors and 17 normal lung specimens. Each specimen has 12600 expression values.

The COIL-100 Image Data Set The data set consists of images of 100 objects. The images of each object are taken every 5 degree by rotating the object, yielding 72 images per object. The dimension of each image is 32 32. In each trial, we randomly select one image as the response vector and use the remaining ones as the data matrix. We run 100 trials and report the average performance of the screening rules.

The description and the experimental settings for the Prostate Cancer data set, the PIE face image data set and the MNIST handwritten digit data set are given in Section 2.3.3.

Fig. 2 reports the rejection ratios of the basic versions of SAFE, DOME, strong rule and EDPP. We can see that EDPP significantly outperforms the other three screening methods on five of the six data sets, i.e., the Colon Cancer, Lung Cancer, Prostate Cancer, MNIST, and COIL-100 data sets. On the PIE face image data set, EDPP and DOME provide similar performance and both significantly outperform SAFE and strong rule.

However, as pointed out by Tibshirani et al. [32], the real strength of screening methods stems from their sequential versions. The reason is because the optimal parameter value is unknown in real applications. Typical approaches for model selection usually involve solving the Lasso problems many times along a sequence of parameter values. Thus, the sequential screening methods are more suitable in facilitating the aforementioned scenario and more useful than their basic-version counterparts in practice [32].

4.1.2 Evaluation of the Sequential EDPP Rule

In this section, we compare the performance of the sequential versions of SAFE, strong rule and EDPP by the rejection ratio and speedup. We first perform experiments on two synthetic data sets. We then apply the three screening rules to six real data sets.

in the sparse learning literature [7, 39, 31]. We simulate data from the true model

We generate two data sets with 250 10000 entries: Synthetic 1 and Synthetic 2. For Synthetic 1, the entries of the data matrix X are i.i.d. standard Gaussian with pairwise correlation zero, i.e., corr() = 0. For Synthetic 2, the entries of the data matrix X are drawn from i.i.d. standard Gaussian with pairwise correlation 0, i.e., corr(. To generate the response vector by the model in (74), we need to set the parameter and construct the ground truth . Throughout this section, is set to be 0.1. To construct , we randomly select p components which are populated from a uniform [1] distribution, and set the remaining ones as 0. After we generate the data matrix X and the response vector y, we run the solver with or without screening rules to solve the Lasso problems along a sequence of 100 parameter values equally spaced on the scale from 0.05 to 1.0. We then run 100 trials and report the average performance.

We first apply the screening rules, i.e., SAFE, strong rule and EDPP to Synthetic 1 with p = 100, 1000, 5000 respectively. Fig. 3(a), Fig. 3(b) and Fig. 3(c) present the corresponding rejection ratios and speedup of SAFE, strong rule and EDPP. We can see that the rejection ratios of strong rule and EDPP are comparable to each other, and both of them are more effective in discarding inactive features than SAFE. In terms of the speedup, EDPP provides better performance than strong rule. The reason is because strong rule is a heuristic screening method, i.e., it may mistakenly discard active features which have nonzero components in the solution. Thus, strong rule needs to check the KKT conditions to ensure the correctness of the screening result. In contrast, the EDPP rule does not need to check the KKT conditions since the discarded features are guaranteed to be absent from the resulting sparse representation. From the last two columns of Table 2, we can observe that the running time of strong rule is about twice of that of EDPP.

Fig. 3(d), Fig. 3(e) and Fig. 3(f) present the rejection ratios and speedup of SAFE, strong rule and EDPP on Synthetic 2 with p = 100, 1000, 5000 respectively. We can observe patterns similar to

Figure 3: Comparison of SAFE, Strong Rule and EDPP on two synthetic datasets with different numbers of nonzero components of the groud truth.

Synthetic 1. Clearly, our method, EDPP, is very robust to the variations of the intrinsic structures of the data sets and the sparsity of the ground truth.

Real Data Sets

In this section, we compare the performance of the EDPP rule with SAFE and strong rule on six real data sets along a sequence of 100 parameter values equally spaced on the scale from 0.05 to 1.0. The data sets are listed as follows:

Table 2: Running time (in seconds) for solving the Lasso problems along a sequence of 100 tuning parameter values equally spaced on the scale of 05 to 1 by (a): the solver [22] (reported in the third column) without screening; (b): the solver combined with different screening methods (reported in the 4columns). The last four columns report the total running time (in seconds) for the screening methods.

We present the rejection ratios and speedup of EDPP, SAFE and strong rule in Fig. 4. Table 3 reports the running time of the solver with or without screening for solving the 100 Lasso problems, and that of the screening rules.

The Breast Cancer Data Set This data set contains 44 tumor samples, each of which is represented by 7129 genes. Therefore, the data matrix 7129. The response vector contains the binary label of each sample.

The Leukemia Data Set This data set is a DNA microarray data set, containing 52 samples and 11225 genes. Therefore, the data matrix 11225. The response vector contains the binary label of each sample.

The SVHN Data set The SVHN data set contains color images of street view house numbers, including 73257 images for training and 26032 for testing. The dimension of each image is 32 In each trial, we first randomly select an image as the response , and then use the remaining ones to form the data matrix . We run 100 trials and report the average performance.

The description and the experiment settings for the Prostate Cancer data set, the PIE face image data set and the MNIST handwritten digit data set are given in Section 2.3.3.

From Fig. 4, we can see that the rejection ratios of strong rule and EDPP are comparable to each other. Compared to SAFE, both of strong rule and EDPP are able to identify far more inactive features, leading to a much higher speedup. However, because strong rule needs to check the KKT conditions to ensure the correctness of the screening results, the speedup gained by EDPP is higher than that by strong rule. When the size of the data matrix is not very large, e.g., the

Figure 4: Comparison of SAFE, Strong Rule, and EDPP on six real data sets.

Breast Cancer and Leukemia data sets, the speedup gained by EDPP are slightly higher than that by strong rule. However, when the size of the data matrix is large, e.g., the MNIST and SVHN data sets, the speedup gained by EDPP are significantly higher than that by strong rule. Moreover, we can also observe from Fig. 4 that, the larger the data matrix is, the higher the speedup can be gained by EDPP. More specifically, for the small data sets, e.g., the Breast Cancer, Leukemia and Prostate Cancer data sets, the speedup gained by EDPP is about 10, 17 and 30 times. In contrast, for the large data sets, e.g., the PIE, MNIST and SVHN data sets, the speedup gained by EDPP is two orders of magnitude. Take the SVHN data set for example. The solver without screening needs about 3 hours to solve the 100 Lasso problems. Combined with the EDPP rule, the solver

Table 3: Running time (in seconds) for solving the Lasso problems along a sequence of 100 tuning parameter values equally spaced on the scale of 05 to 1 by (a): the solver [22] (reported in the second column) without screening; (b): the solver combined with different screening methods (reported in the 3columns). The last three columns report the total running time (in seconds) for the screening methods.

only needs less than 1 minute to complete the task.

Clearly, the proposed EDPP screening rule is very effective in accelerating the computation of Lasso especially for large-scale problems, and outperforms the state-of-the-art approaches like SAFE and strong rule. Notice that, the EDPP method is safe in the sense that the discarded features are guaranteed to have zero coefficients in the solution.

EDPP with Least-Angle Regression (LARS)

As we mentioned in the introduction, we can combine EDPP with any existing solver. In this experiment, we integrate EDPP and strong rule with another state-of-the-art solver for Lasso, i.e., Least-Angle Regression (LARS) [15]. We perform experiments on the same real data sets used in the last section with the same experiment settings. Because the rejection ratios of screening methods are irrelevant to the solvers, we only report the speedup. Table 4 reports the running time of LARS with or without screening for solving the 100 Lasso problems, and that of the screening methods. Fig. 5 shows the speedup of these two methods. We can still observe a substantial speedup gained by EDPP. The reason is that EDPP has a very low computational cost (see Table 4) and it is very effective in discarding inactive features (see Fig. 4).

Figure 5: The speedup gained by Strong Rule and EDPP with LARS on six real data sets.

Table 4: Running time (in seconds) for solving the Lasso problems along a sequence of 100 tuning parameter values equally spaced on the scale of 05 to 1 by (a): the solver [15, 23] (reported in the second column) without screening; (b): the solver combined with different screening methods (reported in the 3columns). The last two columns report the total running time (in seconds) for the screening methods.

4.2 EDPP for the Group Lasso Problem

In this experiment, we evaluate the performance of EDPP and strong rule with different numbers of groups. The data matrix X is fixed to be 250 200000. The entries of the response vector y and the data matrix X are generated i.i.d. from a standard Gaussian distribution. For each experiment, we repeat the computation 20 times and report the average results. Moreover, let denote the number of groups and be the average group size. For example, if is 10000, then

Figure 6: Comparison of EDPP and strong rules with different numbers of groups.

From Figure 6, we can see that EDPP and strong rule are able to discard more inactive groups when the number of groups increases. The intuition behind this observation is that the estimation of the dual optimal solution is more accurate with a smaller group size. Notice that, a large implies a small average group size. Figure 6 also implies that compared to strong rule, EDPP is

Table 5: Running time (in seconds) for solving the group Lasso problems along a sequence of 100 tuning parameter values equally spaced on the scale of 0 by (a): the solver from SLEP (reported in the second column) without screening; (b): the solver combined with different screening methods (reported in the 3columns). The last two columns report the total running time (in seconds) for the screening methods. The data matrix X is of size 250

able to discard more inactive groups and is more robust with respect to different values of

Table 5 further demonstrates the effectiveness of EDPP in improving the efficiency of the solver. When = 10000, the efficiency of the solver is improved by about 80 times. When = 20000 and 40000, the efficiency of the solver is boosted by about 120 and 160 times with EDPP respectively.

5 Conclusion

In this paper, we develop new screening rules for the Lasso problem by making use of the properties of the projection operators with respect to a closed convex set. Our proposed methods, i.e., DPP screening rules, are able to effectively identify inactive predictors of the Lasso problem, thus greatly reducing the size of the optimization problem. Moreover, we further improve DPP rule and propose the enhanced DPP rule, which is more effective in discarding inactive features than DPP rule. The idea of the family of DPP rules can be easily generalized to identify the inactive groups of the group Lasso problem. Extensive numerical experiments on both synthetic and real data demonstrate the effectiveness of the proposed rules. It is worthwhile to mention that the family of DPP rules can be combined with any Lasso solver as a speedup tool. In the future, we plan to generalize our ideas to other sparse formulations consisting of more general structured sparse penalties, e.g., tree/graph Lasso, fused Lasso.

Appendix A.

In this appendix, we give the detailed derivation of the dual problem of Lasso.

A1. Dual Formulation

Assuming the data matrix is , the standard Lasso problem is given by:

For completeness, we give a detailed deviation of the dual formulation of (75) in this section. Note that problem (75) has no constraints. Therefore the dual problem is trivial and useless. A common trick [8] is to introduce a new set of variables such that problem (75) becomes:

By introducing the dual variables , we get the Lagrangian of problem (76):

For the Lagrangian, the primal variables are . And the dual function

In order to get ), we need to solve the following two optimization problems.

and

Let us first consider problem (79). Denote the objective function of problem (79) as

) is convex but not smooth. Therefore let us consider its subgradient

where . In other words, should satisfy

which is equivalent to

Then we plug into Eq. (81):

Therefore, the optimum value of problem (79) is 0.

us rewrite

Combining everything above, we get the dual problem:

which is equivalent to

By a simple re-scaling of the dual variables , i.e., let , problem (86) transforms to:

A2. The KKT Conditions

Problem (76) is clearly convex and its constraints are all affine. By Slater’s condition, as long as problem (76) is feasible we will have strong duality. Denote as optimal primal and dual variables. The Lagrangian is

From the KKT condition, we have

From Eq. (90) and (91), we have:

From Eq. (89), we know there exists

which is equivalent to

From Eq. (93), it is easy to conclude:

Appendix B.

In this appendix, we present the detailed derivation of the dual problem of group Lasso.

B1. Dual Formulation

By introducing the dual variables , the Lagrangian of problem (96) is:

and the dual function

In order to get ), let us solve the following two optimization problems.

and

Let us first consider problem (99). Denote the objective function of problem (99) as

then we can split problem (99) into a set of subproblems. Clearly ˆ) is convex but not smooth because it has a singular point at 0. Consider the subgradient of ˆ

where is the subgradient of

All together, we can conclude the

and thus

The second equality is due to the fact that ’s are independent.

should satisfy:

Next, let us consider problem (100). Since problem (100) is exactly the same as problem (80), we conclude:

Therefore the dual function

Combining everything above, we get the dual formulation of the group Lasso:

which is equivalent to

By a simple re-scaling of the dual variables , i.e., let , problem (104) transforms to:

B2. The KKT Conditions

Clearly, problem (96) is convex and its constraints are all affine. By Slater’s condition, as long as problem (96) is feasible we will have strong duality. Denote as optimal primal and dual variables. The Lagrangian is

From the KKT condition, we have

From Eq. (108) and (109), we have:

From Eq. (107), we know there exists

and

Then the following holds:

for g = 1, 2, . . . , G. Clearly, if , we can conclude

References

[1] U. Alon, N. Barkai, D. Notterman, K. Gish, S. Ybarra, D. Mack, and A. Levine. Broad patterns of gene expression revealed by clustering analysis of tumor and normal colon tissues probed by oligonucleotide arrays. Cell Biology, 96:6745–6750, 1999.

[2] S. Armstrong, J. Staunton, L. Silverman, R. Pieters, M. den Boer, M. Minden, S. Sallan, E. Lander, T. Golub, and S. Korsmeyer. MLL translocations specify a distinct gene expression profile that distinguishes a unique leukemia. Nature Genetics, 30:41–47, 2002.

[3] H. H. Bauschke and P. L. Combettes. Convex Analysis and Monotone Operator Theory in Hilbert Spaces. Springer, 2011.

[4] S. R. Becker, E. Cand`es, and M. Grant. Templates for convex cone problems with applications to sparse signal recovery. Technical report, Standford University, 2010.

[5] D. P. Bertsekas. Convex Analysis and Optimization. Athena Scientific, 2003.

[6] A. Bhattacharjee, W. Richards, J. Staunton, C. Li, S. Monti, P. Vasa, C. Ladd, J. Beheshti, R. Bueno, M. Gillette, M. Loda, G. Weber, E. Mark, E. Lander, W. Wong, B. Johnson, T. Golub, D. Sugarbaker, and M. Meyerson. Classification of human lung carcinomas by mrna expression profiling reveals distinct adenocarcinoma subclasses. Proceedings of the National Academy of Sciences, 98:13790–13795, 2001.

[7] H. Bondell and B. Reich. Simultaneous regression shrinkage, variable selection and clustering of predictors with OSCAR. Biometrics, 64:115–123, 2008.

[8] S. Boyd and L. Vandenberghe. Convex Optimization. Cambridge University Press, 2004.

[9] A. Bruckstein, D. Donoho, and M. Elad. From sparse solutions of systems of equations to sparse modeling of signals and images. SIAM Review, 51:34–81, 2009.

[10] D. Cai, X. He, and J. Han. Efficient kernel discriminant analysis via spectral regression. In ICDM, 2007.

[11] D. Cai, X. He, J. Han, and T. Huang. Graph regularized non-negative matrix factorization for data representation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 33: 1548–1560, 2011.

[12] E. Cand`es. Compressive sampling. In Proceedings of the International Congress of Mathematics, 2006.

[13] S. S. Chen, D. L. Donoho, and M. A. Saunders. Atomic decomposition by basis pursuit. SIAM Review, 43:129–159, 2001.

[14] D. L. Donoho and Y. Tsaig. Fast solution of l-1 norm minimization problems when the solution may be sparse. IEEE Transactions on Information Theory, 54:4789–4812, 2008.

[15] B. Efron, T. Hastie, I. Johnstone, and R. Tibshirani. Least angle regression. Annals of Statistics, 32:407–499, 2004.

[16] L. El Ghaoui, V. Viallon, and T. Rabbani. Safe feature elimination in sparse supervised learning. Pacific Journal of Optimization, 8:667–698, 2012.

[17] J. Fan and J. Lv. Sure independence screening for ultrahigh dimensional feature spaces. Journal of the Royal Statistical Society Series B, 70:849–911, 2008.

[18] J. Friedman, T. Hastie, H. H¨efling, and R. Tibshirani. Pathwise coordinate optimization. Annals of Applied Statistics, 1:302–332, 2007.

[19] J. Friedman, T. Hastie, and R. Tibshirani. Regularization paths for generalized linear models via coordinate descent. Journal of Statistical Software, 33:1–22, 2010.

[20] S. J. Kim, K. Koh, M. Lustig, S. Boyd, and D. Gorinevsky. An interior-point method for large scale l1-regularized least squares. IEEE Journal on Selected Topics in Signal Processing, 1: 606–617, 2007.

[21] Y. Lecun, L. Bottou, Y. Bengio, and P. Haffner. Gradient-based learning applied to document recognition. In Proceedings of the IEEE, 1998.

[22] J. Liu, S. Ji, and J. Ye. SLEP: Sparse Learning with Efficient Projections. Arizona State University, 2009.

[23] J. Mairal, F. Bach, J. Ponce, and G. Sapiro. Online learning for matrix factorization and sparse coding. Journal of Machine Learning Research, 11:19–60, 2010.

[24] S. Nene, S. Nayar, and H. Murase. Columbia object image library (coil-100). Technical report, CUCS-006-96, Columbia University, 1996.

[25] Y. Netzer, T. Wang, A. Coates, A. Bissacco, B. Wu, and A. Ng. Reading digits in nature images with unsupervised feature learning. In NIPS Workshop on Deep Learning and Unsupervised Feature Learning, 2001.

[26] M. Y. Park and T. Hastie. L1-regularized path algorithm for generalized linear models. Journal of the Royal Statistical Society Series B, 69:659–677, 2007.

[27] E. Petricoin, D. Ornstein, C. Paweletz, A. Ardekani, P. Hackett, B. Hitt, A. Velassco, C. Trucco, L. Wiegand, K. Wood, C. Simone, P. Levine, W. Linehan, M. Emmert-Buck, S. Steinberg, E. Kohn, and L. Liotta. Serum proteomic patterns for detection of prostate cancer. Journal of National Cancer Institute, 94:1576–1578, 2002.

[28] A. Ruszczy´nski. Nonlinear Optimization. Princeton University Press, 2006.

[29] S. Shevade and S. Keerthi. A simple and efficient algorithm for gene selection using sparse logistic regression. Bioinformatics, 19:2246–2253, 2003.

[30] T. Sim, B. Baker, and M. Bsat. The CMU pose, illumination, and expression database. IEEE Transactions on Pattern Analysis and Machine Intelligence, 25:1615–1618, 2003.

[31] R. Tibshirani. Regression shringkage and selection via the lasso. Journal of the Royal Statistical Society Series B, 58:267–288, 1996.

[32] R. Tibshirani, J. Bien, J. Friedman, T. Hastie, N. Simon, J. Taylor, and R. Tibshirani. Strong rules for discarding predictors in lasso-type problems. Journal of the Royal Statistical Society Series B, 74:245–266, 2012.

[33] M. West, C. Blanchette, H. Dressman, E. Huang, S. Ishida, R. Spang, H. Zuzan, J. Olson, J. Marks, and J. Nevins. Predicting the clinical status of human breast cancer by using gene expression profiles. Proceedings of the National Academy of Sciences, 98:11462–11467, 2001.

[34] J. Wright, Y. Ma, J. Mairal, G. Sapiro, T. Huang, and S. Yan. Sparse representation for computer vision and pattern recognition. In Proceedings of IEEE, 2010.

[35] Z. J. Xiang and P. J. Ramadge. Fast lasso screening tests based on correlations. In IEEE ICASSP, 2012.

[36] Z. J. Xiang, H. Xu, and P. J. Ramadge. Learning sparse representation of high dimensional data on large scale dictionaries. In NIPS, 2011.

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

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

[39] H. Zou and T. Hastie. Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society Series B, 67:301–320, 2005.

designed for accessibility and to further open science