# On The Projection Operator to A Three-view Cardinality Constrained Set

The cardinality constraint is an intrinsic way to restrict the solution structure in many domains, for example, sparse learning, feature selection, and compressed sensing. To solve a cardinality constrained problem, the key challenge is to solve the projection onto the cardinality constraint set, which is NP-hard in general when there exist multiple overlapped cardinality constraints. In this paper, we consider the scenario where the overlapped cardinality constraints satisfy a Three-view Cardinality Structure (TVCS), which reflects the natural restriction in many applications, such as identification of gene regulatory networks and task-worker assignment problem. We cast the projection into a linear programming, and show that for TVCS, the vertex solution of this linear programming is the solution for the original projection problem. We further prove that such solution can be found with the complexity proportional to the number of variables and constraints. We finally use synthetic experiments and two interesting applications in bioinformatics and crowdsourcing to validate the proposed TVCS model and method.

## Authors

• 10 publications
• 9 publications
• 7 publications
• 17 publications
• 9 publications
• 101 publications
• ### On Minimal Constraint Networks

In a minimal binary constraint network, every tuple of a constraint rela...
03/08/2011 ∙ by Georg Gottlob, et al. ∙ 0

• ### Synthesizing Linked Data Under Cardinality and Integrity Constraints

The generation of synthetic data is useful in multiple aspects, from tes...
03/26/2021 ∙ by Amir Gilad, et al. ∙ 0

• ### A Comparison of Lex Bounds for Multiset Variables in Constraint Programming

Set and multiset variables in constraint programming have typically been...
06/29/2011 ∙ by Yat-Chiu Law, et al. ∙ 0

• ### Integrating Cardinality Constraints into Constraint Logic Programming with Sets

Formal reasoning about finite sets and cardinality is an important tool ...
02/10/2021 ∙ by Maximiliano Cristiá, et al. ∙ 0

• ### Reflections on "Incremental Cardinality Constraints for MaxSAT"

To celebrate the first 25 years of the International Conference on Princ...
10/10/2019 ∙ by Ruben Martins, et al. ∙ 0

• ### Faster Maximum Feasible Subsystem Solutions for Dense Constraint Matrices

Finding the largest cardinality feasible subset of an infeasible set of ...
02/10/2021 ∙ by Fereshteh Fakhar Firouzeh, et al. ∙ 0

• ### Fair Robust Assignment using Redundancy

We study the consideration of fairness in redundant assignment for multi...
02/11/2021 ∙ by Matthew Malencia, et al. ∙ 0

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## 1 Introduction

The cardinality constraint is an intrinsic way to restrict the solution structure in many real problems, for example, sparse learning (Olshausen and Field, 1997), feature selection (Zhang, 2009), and compressed sensing (Candes et al., 2006). The generic cardinality constrained optimization can be expressed as

 minw∈Rp f(w) (1a) subject to ∥wg∥0≤sg∀g∈G (1b)

where is the optimization variable, is an index subset of ,

is the sub-vector of

indexed by . denotes the cardinality of the sub-vector, i.e., the number of nonzeros in , is the hyper set of all predefined groups, and is the upper bound vector - refers to the upper bound of the sparsity over group . Objective

is the loss function which could be defined with different form according to the specific application. The problem (

1) refers to a nonconvex optimization (NP-hard) due to the cardinality constraint. Some efficient iterative methods such as IHT (Yuan et al., 2014), CoSaMP (Needell and Tropp, 2009), GradMP (Nguyen et al., 2012), and their variants can guarantee to solve the original problem under some mild conditions. A key component in all of these methods is the projection operator

 PΩ(G,s)(v):=argminw∈Ω(G,s) ∥w−v∥2 (2)

where denotes the feasible set to the constraint (1b). While in some special case, for example, , the projection is trivial, it is quite challenging, especially when includes multiple overlapped index sets (even NP-hard in some cases).

In this paper, we consider the scenario where the overlapped cardinality constraints (1b) satisfy a Three-view Cardinality Structure (TVCS):

###### Definition 1.

(Three-view Cardinality Structure (TVCS)) For , the hyper set consisting of subsets of admits the TVCS structure if the following conditions are satisfied:

• There exists a partition , and such that ;

• ;

• All element sets in have no overlap;

• All element sets in have no overlap.

This definition basically requires that can be partitioned into three hyper sets , , and , and overlaps can only happen between element sets in different hyper sets. is usually used to restrict the overall sparsity. Figure 1 provides two examples of for TVCS.

The TVCS model is motivated from some important applications, for example, in recommendation, task-worker assignment, and bioinformatics.

• Online recommendation. Suppose we want to recommend a certain number of books (among books) to a customer - corresponding to the based sparsity constraint. Among the selected books, we want to maintain some diversities - the recommended books by the same author should not exceed a certain number ( based sparsity constraint) and about the same topic should not exceed a certain number either ( based sparsity constraint). One can refer to the top graph in Figure 1: is grouped by authors and is grouped by topics.

• Task-worker assignment. Suppose we have a bunch of tasks and workers, and we want to assign the tasks to workers. For example, in crowdsourcing, we usually assign several different workers to each task since we want to use the answers from multiple workers to improve the accuracy. On the other hand, each worker is usually assigned to multiple tasks so there is a “many to many” relationship in this assignment. The goal is to pursue the optimal assignment under a certain criteria in crowdsourcing, while satisfying some restrictions. For example, the total assignments should be bounded by the total budget (corresponding to ), the total cost of assignments to a single worker cannot exceed a certain threshold (corresponding to ), and the total cost of assignments on a single task cannot exceed a certain threshold (corresponding to ). Let be the assignment matrix, and its rows are indexed by workers and the columns are indexed by tasks. These constraints can be illustrated by the bottom graph in Figure 1.

• Identification of gene regulatory networks. The essential goal of identifying gene regulatory network is to identify a weighted directed graph, which can be represented by a square matrix with elements in total where is the number of vertices. A sparse network constraint is to restrict the in-degree and out-degree for each vertex, which corresponds to the sparsity in each row and column of .

To solve the TVCS constrained projection (2), we show an interesting connection between the projection and a linear programming (LP) that the vertex solution to this linear programming is an integer solution which solves the original problem.

To find an integer solution to such LP efficiently, we formulate it into a feasibility problem, and further an equivalent quadratic convex optimization. By using the rounding technique, we can avoid solving the exact solution of this LP problem. We propose an iterative algorithm to solve it and each iteration can be completed in linear time. We also show that the iterate linearly converges to the optimal point. Finally, the proposed TVCS model is validated by the synthetic experiment and two important and novel applications in identification of gene regulatory networks and task assignment problem of crowdsourcing.

## 2 Related Works

Recent years have witnessed many research works in the field of structured sparsity and group-based sparsity. Yuan and Lin (2006) introduced the group LASSO, which pursues group-wise sparsity that restricts the number of groups for the selected variables. Jenatton et al. (2011) construct a hierarchical structure over the variables and use group LASSO with overlapped groups to solve it. Exclusive LASSO (Zhou et al., 2010; Kong et al., 2014) was proposed for the exclusive group sparsity which can be treated as relaxing our cardinality constraints to convex regularizations. In (Kong et al., 2014), the authors discussed the overlapping situation and tried to solve the problem using convex relaxation, which is different from our approach. Besides the aforementioned works, some proposed more general models to cover various sparsity structures. Bach et al. (2012) extended the usage of -norm relaxation to several different categories of structures. And recently, another generalization work (El Halabi and Cevher, 2015) proposed convex envelopes for various sparsity structures. They built the framework by defining a totally unimodular penalty, and showed how to formulate different sparsity structures using the penalty. The work above concentrated on using convex relaxation to control the sparsity.

Besides using convex relaxation, there are also some works focusing on projection-based methods. When the exact projection operator was provided, Baraniuk et al. (2010) extended the traditional IHT and CoSaMP methods to general sparsity structures. In this work, the authors also introduced the projection operator for block sparsity and tree sparsity. Cevher et al. (2009) investigated cluster sparsity and they applied dynamic programming to solve the projection operator for their sparsity model. Hegde et al. (2009) introduced a “spike trains” signal model, which is also related to exclusive group sparsity. Its groups always have consecutive coordinates, and each group cannot contain more than one nonzero element. To solve the projection problem of their model, they showed the basic feasible solutions of the relaxed linear programming (LP) are always integer points. In our work, we also use LP to solve the projection problem. But our model defines the group structure differently and aims at different applications.

In addition, there are some works for the cases without an efficient exact projection operator (Hegde et al., 2015a, b; Nguyen et al., 2014)

. This is meaningful since the projection operator for complex structured sparsity often involves solving complicated combinatorial optimization problems.

Hegde et al. (2015a) discussed how to guarantee convergence if using approximate projection in IHT and CoSaMP for compressive sensing. They proved that the convergence needs a “head approximation” to project the update (gradient) before applying it. Hegde et al. (2015b) proposed a general framework to formulate a series of models as a weighted graph, and designed an efficient approximate projection operator for the models. Nguyen et al. (2014) applied the approximate projection-based IHT and CoSaMP to general convex functions and stochastic gradients.

## 3 Preliminary: GradMP and IHT Frameworks

This section briefly reviews two commonly used algorithm frameworks to solve the cardinality constrained optimization (1): iterative hard thresholding (IHT) (Yuan et al., 2014; Nguyen et al., 2014) and gradient matching pursuit (GradMP) (Nguyen et al., 2012, 2014) (the general version of CoSaMP (Needell and Tropp, 2009)) for solving cardinality constrained problem. Other methods like hard thresholding pursuit (HTP) also follows similar steps and has been shown to be effective both empirically and theoretically (Yuan et al., 2016). The procedures of IHT and GradMP for our model are shown in Algorithms LABEL:alg:stoiht and LABEL:alg:gradmp, where is the support set of the argument vector.

Therefore, one can see that the efficiency of both algorithms relies on the computation of the gradient and the projection. To avoid the expensive computation of the gradient, GradMP and IHT can be extended to the stochastic versions (Nguyen et al., 2014) by assigning the stochastic gradient at the gradient computation step.

Both Algorithms LABEL:alg:stoiht and  LABEL:alg:gradmp (and their stochastic variants) guarantee some nice properties: the iterate converges to a small ball surrounding the true solution at a linear rate under certain RIP-type conditions (Nguyen et al., 2014) and the radius of such ball converges to zero when the number of samples goes to infinity.

algocf[htbp]

algocf[htbp]

A common component in Algorithms LABEL:alg:stoiht and LABEL:alg:gradmp is the projection operator. If all the groups except in do not overlap each other, the projection problem can be easily solved by sequential projections (Yang et al., 2016). But for those cases involving overlapped groups, it is generally challenging to solve them efficiently.

## 4 Projection Operator

This section introduces how to solve the essential projection step. Note that the projection onto a nonconvex set is NP-hard in general. By utilizing the special structure of TVCS, we show that the projection can be solved efficiently. Due to the page limitation, all proofs are provided in the supplementary material.

### 4.1 LP Relaxation

Firstly, we can cast the projection problem (2) to an equivalent integer linear programming problem (ILP) according to Lemma 1.

###### Lemma 1.

The projection problem (2) is equivalent to the following integer linear programming problem (ILP):

 maxx ⟨v2,x⟩ (3) subject to Ax≤s x∈{0,1}p

where is applying element-wise square operation on vector . is a matrix which is defined as:

 A=[1⊤C] (4)

where , whose rows represent the indicator vector of each group and .

Each row in corresponds to one group from . For example, if the -th coordinate is in the -th group, otherwise . The first row corresponds to the overall sparsity i.e. .

It is NP-hard to solve an ILP in general. One common way to handle such ILP is making a linear programming (LP) relaxation. In our case, we can use a box constraint to replace the integer constraint :

 maxx ⟨v2,x⟩ (5) subject to Ax≤s x∈[0,1]p

However, there is no guarantee that a general ILP can be solved via its LP relaxation, because the solution of the relaxed LP is not always integer. Although one can make a rounding to the LP solution and acquire a integer solution, such solution is not guaranteed to be optimal (or even feasible) to the original ILP.

Fortunately, due to the special structure of our TVCS model, we find that its relaxed LP has some nice properties which make it possible to get the optimal solution of the ILP efficiently. The following theorem reveals the relationship between the ILP problem and the relaxed LP problem.

###### Theorem 2.

Given satisfying , all the vertices of the feasible set to (5) are integer points. Furthermore, there is an optimal solution on the vertex that solves the ILP (3).

This theorem suggests that finding a vertex solution of the relaxed LP can solve the original projection problem onto a TVCS . The proof basically shows that matrix (for TVCS) is a totally unimodular matrix (Papadimitriou and Steiglitz, 1982). We provide the detailed proof in the supplementary material.

### 4.2 Linearly Convergent Algorithm for Projection Operator onto TVCS

To find a solution on the vertex, one can use the Simplex method. Although Simplex method guarantees to find an optimal solution on the vertex and could be very efficient in practice, it does not have a deterministic complexity bound. In the IHT and GradMP algorithms, projection operator is only a sub-procedure in one iteration. Hence, we are usually supposed to solve lots of instances of problem (3). Simplex might be efficient practically, but its worst case may lead to exponential time complexity (Papadimitriou and Steiglitz, 1982). In this section, the integer solution to the linear programming can be found within the complexity proportional to the number of variables and constraints.

#### Equivalent Feasibility Problem Formulation.

The dual of LP problem (5) can be written as:

 miny ⟨[s⊤ 1⊤]⊤,y⟩ (6) subject to [A⊤ I]y≥v2,y≥0

Since the duality gap of LP is zero, combining the primal LP (5) and dual LP (6), we can formulate an equivalent problem, i.e. the feasibility problem over the following constraints:

 find x,y subject to ⟨[s⊤ 1⊤]⊤,y⟩=⟨v2,x⟩ [A⊤ I]y≥v2 [AI]x≤[s1] y≥0,x≥0

#### Iterative Algorithm.

The feasibility problem with linear constraints above is equivalent to the following optimization problem:

 minx,y 12(⟨[s⊤ 1⊤]⊤,y⟩−⟨v2,x⟩)2 +12∥[v2−[A⊤ I]y]+∥2+12∥[Ax−s]+∥2 (7) subject to 0≤x≤1,y≥0

where is the element-wise hinge operator, i.e. it transforms each element to .

This is a convex optimization problem with a quadratic objective and box constraints. We adopt the projected gradient descent to solve this problem, and show it converges linearly.

###### Theorem 3.

For the optimization problem with the form

 minz f(z):=∥[Az−a]+∥2+∥Bz−b∥2 subject to z∈Ω

where , the projected gradient descent algorithm has a linear convergence rate with some (depending on and ):

 ∥zt+1−Pz∗(zt+1)∥≤α∥zt−Pz∗(zt)∥,

where is the projection onto the optimal solution set.

Notice that the objective function in Theorem 3 is not necessarily strongly convex, which means the well recognized linear convergence conclusion from the strong convexity is not applicable here.

Theorem 3 mainly applies Hoffman’s Theorem (Hoffman, 2003) to show that is an optimal strongly convex function (Liu and Wright, 2015). This leads to a linear convergence rate.

The convergence rate , where is the Hoffman constant (Hoffman, 2003) that depends on and is always positive. is the Lipschitz continuous gradient constant. More details are included in the supplementary materials.

To show the complexity of this algorithm, we firstly count how many iterations we need. Since we know that we can just make a rounding111

Acute readers may notice that the convergent point may be on the face of the polytope in some cases instead of vertex. However, we can add a small random perturbation to ensure the optimal point to be vertices with probability 1.

to the result when we attain . Let represent all the variables in (7). Because , we can do the rounding safely when , where are the optimal points of this problem. According to Theorem 3, we have the linear convergence rate , so the number of iterations we need is

 t>logα12∥z0−z∗∥

Therefore, we claim that we can obtain the solution by rounding after iterations.

Secondly, we show that the computation complexity in each iteration is linear with dimensionality and the amount of groups . Since each column of contains at most 3 nonzero elements, the complexity of the matrix multiplications in computing the gradient of (7) is . Together with other computation, the complexity for each iteration is .

## 5 Empirical Study

This section will validate the proposed method on both synthetic data and two practical applications: crowdsourcing and identification of gene regulatory networks.

### 5.1 Linear Regression and Classification on Synthetic Data

In this section, we validate the proposed

method with linear regression objective and squared hinge objective (classification) on

synthetic data. Let be a matrix, and are defined as groups with all rows and all columns respectively. The linear regression loss is defined as and the squared hinge loss is defined as , where is the total number of training samples. and are the features and label of the -th sample respectively.

In the linear regression experiment, the true model

is generated from the following procedure: generate a random vector and apply the projection operator to get a support set which satisfy our sparsity constraints; the elements of positions in support set are drawn from standard normal distribution.

is fixed as and is gradually increased. The group sparsity upper bounds for and are uniformly generated from the integers in the range. The overall sparsity upper bound is set by . Each ’s is an

i.i.d. Gaussian random matrix.

is generated from , where is the i.i.d. Gaussian random noise drawn from . We compare the proposed methods to bilevel exclusive sparsity with non-overlapped groups (row-wise or column-wise) (Yang et al., 2016), overall sparsity (Needell and Tropp, 2009), and exclusive LASSO (Kong et al., 2014). For fairness we project the final result of all the compared methods to satisfy all constraints. All the experiments are repeated 30 times and we use the averaged result. We use selection recall and successful recovery rate to evaluate the performance. Selection recall is defined as , where is the optimization result. Successful recovery rate is the ratio of the successful feature selection i.e. to the total number of repeated experiments. In Figure 2 we can observe that our model with all sparsity constraints always have the best performance. While the performance of exclusive LASSO and our method is comparable when the number of samples are very limited, our method outperforms exclusive LASSO when the number of samples increases.

For classification experiments, we use the same settings of sparsity with linear regression. Here we set , and change from to . The true model and feature matrices are generated by the same way as the linear regression experiment. The class label is got by . Besides the selection recall, we also compare the classification error. In Figure 3, we can see that the superiority of our method is even more significant in the classification experiment. Although the overall sparsity has the lowest selection recall, it still has a similar classification error with the methods that consider row or column groups.

### 5.2 Application in Crowdsourcing

This section applies the proposed method to the worker-task assignment problem in crowdsourcing. Take the image labeling task as an example. Given workers and images, each image can be assigned to multiple workers and each worker can label multiple images. The predicted label for each image is decided by all the labels provided by the assigned workers and the quality of each worker on the image. The goal is to maximize the expected prediction accuracy based on the assignment. Let be the assignment matrix, i.e. if assign the -th worker to -th task, otherwise .

is the corresponding quality matrix, which is usually estimated from the golden standard test

(Ho et al., 2013). The whole formulation is defined to maximize the average of the expected prediction accuracy over tasks over a TVCS constraint:

 maxX 1mm∑j=1Eacc(Q⋅,j,X⋅,j) (8) subject to n∑i=1Xij≤sworker,∀j; m∑j=1Xij≤stask,∀i; ∑i,jXij≤stotal; X∈{0,1}n×m

where is the expected prediction accuracy, is the “worker sparsity”, i.e. the largest number of assigned workers for each task, and is the “task sparsity”, i.e. each worker can be assigned with at most tasks, and is the total sparsity to control the budget, i.e., the maximal number of assignment. In image labeling task, we assume that each image can only have two possible classes and the percentage of images in each class is one half. We use the Bayesian rule to infer the predicted labels given the workers’ answer. Here we consider the binary classification task. Let be the true label of the -th task and be the prediction given labels by selected workers, i.e.,

 ^yj={0, if P(yj=1|^YΩj,j)

where is the -th worker’s predication on -th task. Set contains the indices of the selected workers for -th task, i.e. , and

Then can be defined in the following:

 Eacc(Q⋅,j,X⋅,j)=P(^yj=1,yj=1)+P(^yj=0,yj=0)

By this way, the expected accuracy will not be continuous, so we use smooth function to approximate the expected accuracy and adopt the stochastic gradient with the proposed projection operator to optimize it. Due to the space limitation, the detailed derivation of the objective formulation can be found in the supplemental material.

We conduct experiment for crowdsourcing task assignment on synthetic data. Specifically, we generate the quality matrix from uniformly random distribution with interval

. The prior probability

and are set as for all the tasks.

To avoid evaluating the expectation term, we apply the stochastic iterative hard thresholding framework (Nguyen et al., 2014). Each iteration we get and by sampling based on i.e. , . Then we can get a stochastic gradient based on the sampled .

We evaluate the experiments on different value of by setting them as different ratios of the total number of tasks and workers. The overall sparsity is set by the same way as in Section 5.1. To measure the performance, we compare the sampled expected accuracy. The samples (i.e., ) are independent to the samples used in training. Figure 4 shows the comparison of the expected accuracy of the three approaches. We can observe that the accuracy increases with larger ratio (i.e. more assignments). The random assignment strategy needs more assignments to get the same accuracy compared with the other two methods.

### 5.3 Application in Identification of Gene Regulatory Networks

In this section, we apply the projection operator to the identification of gene regulatory networks (GRN).

#### Background.

Gene regulatory network represents the relations between different genes, which plays important roles in biological processes and activities by controlling the expression level of RNAs. There is a well-known biological competition named DREAM challenge about identifying

GRN. Based on the time series gene expression data which are RNAs’ level along time sequence, contestants are required to recover the whole gene network of given size. One popular way to infer GRN is to utilize the sparse property of GRN: e.g., one gene in the network is only related to a small number of genes and we already know that there exists no relationship between some genes. Therefore, the amount of edges connecting to one vertex is far less than the dimension of the graph. It is a practical case of row-wise and column-wise sparsity for matrix. We could apply the projection operator to constrain the number of edges related to each vertex to identify the whole network. Recently, the dynamic Bayesian network (DBN)

(Zou and Conzen, 2005) is supposed to be an effective model to recover GRNs. The RNAs’ level of all genes in GRN at time is stored in gene expression vector , where each entry corresponds to one gene respectively and is the number of genes in GRN. Hence, We define the total amount of time points in the experiment as . Gene activity model is usually assumed to be

 xt+1=Pxt+et,t=1…T−1,

where is the covariance matrix of GRN and

is Gaussian white noise. Then the difference of RNA levels between time points

and , i.e. is defined as follows:

 yt+1,t:=xt+1−xt=¯Wxt+et,t=1…T−1,

where is the true sparse -by- matrix. Therefore, the GRN is only considered between different genes and we eliminate edges whose start and end vertex are the same. We define that and . The objective function is

 f(W)=12∥Y−¯WX∥2F=12T−1∑t=1∥(xt+1−xt)−¯Wxt∥2.

#### Time-course Gene Expression Data.

To evaluate our method, we employ GeneNetWeaver (Marbach et al., 2009; Schaffter et al., 2011)

, the official DREAM Challenge tool for time-series expression data generation. With typical gene network structure and ordinary differential equation (ODE) models, GeneNetWeaver will produce the time-course gene expression data at pre-specified time points. In the simulation studies, we control the size of gene network to be

vertexes and the gene expression data are generated under 10% Gaussian white noise.

The network is shown in Figure 5. In this Figure, it is clear that one gene only has a few connections to other genes. Therefore, the GRN is sparse and we are able to restrict the in-degree and out-degree of every vertex by representing the network as a matrix and controlling the sparsity within each row and column.

#### Performance evaluation.

Six commonly-used criteria are considered to measure the performance, i.e., sensitivity (SN), specificity (SP), accuracy (ACC), F-measure, Matthews correlation coefficient (MCC), and the Area Under ROC Curve (AUC):

 SN=TPTP+FN,SP=TNTN+FP, ACC=TP+TNTP+FP+% TN+FN, F-measure=2×SN×SPSN+SP, MCC=TP×TN−FP×FN% √(TP+FP)(TP+FN)(TN+FP)(%TN+FN),

where TP and TN denote the true positive and true negative, and FP and FN denote the false positive and false negative, respectively. With these criteria, we compare the performance of our method with six representative algorithms, including PCC, ARACNE (Margolin et al., 2006), CLR (Faith et al., 2007), MINET (Meyer et al., 2008), GENIE3 (Huynh-Thu et al., 2010), TIGRESS (Haury et al., 2012). The results are summarized in Table 1. Our method outperforms other six state-of-art methods: the AUC of our method achieve 0.7 higher which is far more than other methods; 5 out of 6 different measure show that our method has significant advantage compared to other algorithms.

## 6 Conclusion

This paper considers the TVCS constrained optimization, motivated by the intrinsic restrictions for many important applications, for example, in bioinformatics, recommendation system, and crowdsourcing. To solve the cardinality constrained problem, the key step is the projection onto the cardinality constraints. Although the projection onto the overlapped cardinality constraints is NP-hard in general, we prove that if the TVCS condition is satisfied the projection can be reduced to a linear programming. We further prove that there is an iterative algorithm which finds an integer solution to the linear programming within time complexity , where is the distance from the initial point to the optimization solution and is the convergence rate. We finally use synthetic experiments and two interesting applications in bioinformatics and crowdsourcing to validate the proposed TVCS model.

## Acknowledgements

This project is supported in part by the NSF grant CNS-1548078 and the NEC fellowship.

## References

• Bach et al. [2012] F. Bach, R. Jenatton, J. Mairal, G. Obozinski, et al. Structured sparsity through convex optimization. Statistical Science, 27(4):450–468, 2012.
• Baraniuk et al. [2010] R. G. Baraniuk, V. Cevher, M. F. Duarte, and C. Hegde. Model-based compressive sensing. Information Theory, IEEE Transactions on, 56(4):1982–2001, 2010.
• Candes et al. [2006] E. J. Candes, J. K. Romberg, and T. Tao. Stable signal recovery from incomplete and inaccurate measurements. Communications on pure and applied mathematics, 59(8):1207–1223, 2006.
• Cevher et al. [2009] V. Cevher, P. Indyk, C. Hegde, and R. G. Baraniuk. Recovery of clustered sparse signals from compressive measurements. Technical report, DTIC Document, 2009.
• El Halabi and Cevher [2015] M. El Halabi and V. Cevher. A totally unimodular view of structured sparsity. In

Proceedings of the Eighteenth International Conference on Artificial Intelligence and Statistics

, pages 223–231, 2015.
• Faith et al. [2007] J. J. Faith, B. Hayete, J. T. Thaden, I. Mogno, J. Wierzbowski, G. Cottarel, S. Kasif, J. J. Collins, and T. S. Gardner. Large-scale mapping and validation of escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS biol, 5(1):e8, 2007.
• Haury et al. [2012] A.-C. Haury, F. Mordelet, P. Vera-Licona, and J.-P. Vert. Tigress: trustful inference of gene regulation using stability selection. BMC systems biology, 6(1):145, 2012.
• Hegde et al. [2009] C. Hegde, M. F. Duarte, and V. Cevher. Compressive sensing recovery of spike trains using a structured sparsity model. In SPARS’09-Signal Processing with Adaptive Sparse Structured Representations, 2009.
• Hegde et al. [2015a] C. Hegde, P. Indyk, and L. Schmidt. Approximation algorithms for model-based compressive sensing. Information Theory, IEEE Transactions on, 61(9):5129–5147, 2015a.
• Hegde et al. [2015b] C. Hegde, P. Indyk, and L. Schmidt. A nearly-linear time framework for graph-structured sparsity. In

Proceedings of the 32nd International Conference on Machine Learning (ICML-15)

, pages 928–937, 2015b.
• Ho et al. [2013] C.-J. Ho, S. Jabbari, and J. W. Vaughan. Adaptive task assignment for crowdsourced classification. In Proceedings of The 30th International Conference on Machine Learning, pages 534–542, 2013.
• Hoffman [2003] A. J. Hoffman. On approximate solutions of systems of linear inequalities. In Selected Papers Of Alan J Hoffman: With Commentary, pages 174–176. 2003.
• Huynh-Thu et al. [2010] V. A. Huynh-Thu, A. Irrthum, L. Wehenkel, and P. Geurts. Inferring regulatory networks from expression data using tree-based methods. PloS one, 5(9):e12776, 2010.
• Jenatton et al. [2011] R. Jenatton, J. Mairal, G. Obozinski, and F. Bach. Proximal methods for hierarchical sparse coding. Journal of Machine Learning Research, 12(Jul):2297–2334, 2011.
• Kong et al. [2014] D. Kong, R. Fujimaki, J. Liu, F. Nie, and C. Ding. Exclusive feature learning on arbitrary structures via -norm. In Advances in Neural Information Processing Systems, pages 1655–1663, 2014.
• Liu and Wright [2015] J. Liu and S. J. Wright. Asynchronous stochastic coordinate descent: Parallelism and convergence properties. SIAM Journal on Optimization, 25(1):351–376, 2015.
• Marbach et al. [2009] D. Marbach, T. Schaffter, C. Mattiussi, and D. Floreano. Generating realistic in silico gene networks for performance assessment of reverse engineering methods. Journal of computational biology, 16(2):229–239, 2009.
• Margolin et al. [2006] A. A. Margolin, I. Nemenman, K. Basso, C. Wiggins, G. Stolovitzky, R. D. Favera, and A. Califano. Aracne: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC bioinformatics, 7(Suppl 1):S7, 2006.
• Meyer et al. [2008] P. E. Meyer, F. Lafitte, and G. Bontempi. minet: Ar/bioconductor package for inferring large transcriptional networks using mutual information. BMC bioinformatics, 9(1):461, 2008.
• Needell and Tropp [2009] D. Needell and J. A. Tropp. Cosamp: Iterative signal recovery from incomplete and inaccurate samples. Applied and Computational Harmonic Analysis, 26(3):301–321, 2009.
• Nguyen et al. [2012] N. Nguyen, S. Chin, and T. D. Tran. A unified iterative greedy algorithm for sparsity constrained optimization. 2012.
• Nguyen et al. [2014] N. Nguyen, D. Needell, and T. Woolf. Linear convergence of stochastic iterative greedy algorithms with sparse constraints. arXiv preprint arXiv:1407.0088, 2014.
• Olshausen and Field [1997] B. A. Olshausen and D. J. Field. Sparse coding with an overcomplete basis set: A strategy employed by v1? Vision research, 37(23):3311–3325, 1997.
• Papadimitriou and Steiglitz [1982] C. H. Papadimitriou and K. Steiglitz. Combinatorial optimization: algorithms and complexity. Courier Corporation, 1982.
• Schaffter et al. [2011] T. Schaffter, D. Marbach, and D. Floreano. Genenetweaver: in silico benchmark generation and performance profiling of network inference methods. Bioinformatics, 27(16):2263–2270, 2011.
• Yang et al. [2016] H. Yang, Y. Huang, L. Tran, J. Liu, and S. Huang. On benefits of selection diversity via bilevel exclusive sparsity. In Computer Vision and Pattern Recognition (CVPR), 2016 IEEE Conference on. IEEE, 2016.
• Yuan and Lin [2006] M. Yuan and Y. Lin. Model selection and estimation in regression with grouped variables. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 68(1):49–67, 2006.
• Yuan et al. [2014] X. Yuan, P. Li, and T. Zhang. Gradient hard thresholding pursuit for sparsity-constrained optimization. In Proceedings of The 31st International Conference on Machine Learning, pages 127–135, 2014.
• Yuan et al. [2016] X. Yuan, P. Li, and T. Zhang. Exact recovery of hard thresholding pursuit. In Advances in Neural Information Processing Systems, pages 3558–3566, 2016.
• Zhang [2009] T. Zhang. On the consistency of feature selection using greedy least squares regression. Journal of Machine Learning Research, 10(Mar):555–568, 2009.
• Zhou et al. [2010] Y. Zhou, R. Jin, and S. C. Hoi. Exclusive lasso for multi-task feature selection. In AISTATS, volume 9, pages 988–995, 2010.
• Zou and Conzen [2005] M. Zou and S. D. Conzen. A new dynamic bayesian network (dbn) approach for identifying gene regulatory networks from time course microarray data. Bioinformatics, 21(1):71–79, 2005.

## 1 Proof of Lemma 1

Firstly we show how to convert the projection problem (2) to a support set selection problem. For any vector , let vector indicate the nonzero positions of , then we can claim that

 ∥v−vx∥2≤∥v−w∥2

where is a vector having same dimension with , and it keeps elements at positions where has “1”, and fills zeros at positions where has “0”. In addition, vector if and only if its support set indicator vector satisfies , given is defined in (4).

So the problem (2) can be converted to integer programming:

 minx∈{0,1}p ∥v−vx∥2 (9) subject to Ax≤s

and the objective can be further simplified:

 ∥v−vx∥2 =⟨v2,1−x⟩ =∥v∥2−⟨v2,x⟩

Since is constant here, then the problem (9) is equivalent to the ILP (3), which means problem (2) is equivalent to ILP (3). Then we complete the proof.

## 2 Proof of Theorem 2

To prove Theorem 2, we use the concept of totally unimodular matrix.

###### Definition 2.

Totally Unimodular (TU) Matrix. An integer matrix is TU, if the determinant of any square submatrices333Submatrix here is a square smaller matrix obtained by removing certain rows and columns is in the set .

###### Proposition 1.

If is TU, then is TU, and their concatenations with identity matrices (i.e. ) are still TU.

###### Proof.

Since transposing matrix will not change the determinant, so it is obvious is TU.

Then we prove stacking with identity matrix

preserves TU property. We prove it by induction. Firstly, we show that submatrix with size 1 always has determinant in , because any element from is either 1 or 0. Now consider a submatrix with size having determinant in , then submatrix with size will still have determinant in . To show this, we only need to prove that adding a new row/column from will not change the determinant out of set . Since any row/column from only has one nonzero element “1”, we can eliminate other elements in the same position by subtracting a multiple of this row/column to other rows/columns. After that, we can remove this row and column, and the determinant can only change its sign. So we know that submatrix with size has determinant in if submatrix with size has determinant in . ∎

###### Lemma 4.

If is TU, is an integer vector, then all vertices of the following polytope are integer points:

 {x∣Ax≤s,x∈[0,1]p} (10)
###### Proof.

We can have an equivalent form of this polytope:

 {x∣∣∣[AI]x≤[s1],x≥0} (11)

From Proposition 1, we know that matrix

 [AI]

is TU if is TU, so this meets the case in Theorem 13.2 (see Papadimitriou and Steiglitz, 1982, chap 13). Then we complete the proof. ∎

###### Lemma 5.

If is the matrix whose each row is the indicator vector of a group in our TVCS model, then is a TU matrix.

###### Proof.

Since ’s rows are the indicator vectors of groups in , so . From the definition 1, we know that there are at most two “1”s in each column. For the column which has two “1”s, the two corresponding groups are from and respectively. (Because we know that groups within or do not overlap.)

By this way, our matrix meets the case in Theorem 13.3 (see Papadimitriou and Steiglitz, 1982, chap 13), and it is a TU matrix. ∎

###### Lemma 6.

Recall the matrix in equation (4):

 A=[1⊤C]

If it is constructed from the TVCS model, then is a TU matrix.

###### Proof.

From Lemma 5, we know is a TU matrix for any of our TVCS model. In other words, any submatrix restricted in has determinant , , or . Therefore, we only need to consider the submatrix of has overlaps with the first row . There are only three possible forms of such submatrix . We will show all of their determinants are in .

1. At least one column of has a single “”, so it must appear in the first row . By exchanging such column with the last column (which can only influence the sign of determinant), we can transform it with form:

 [1⊤1¯C0]

where is any submatrix of . From the matrix determinant property, we have . Therefore, submatrix in such form have determinants in .

2. All columns of have three “” elements (the last row has “1” for every column). For the rows which are from , we can sum all the rows to a certain row (this will not change the determinant). By this way we transform to the following form:

 ⎡⎢⎣1⊤2⊤¯C⎤⎥⎦

where is a submatrix of . In this case, is not full rank, so its determinant is 0.

3. Each column in contains at least two “” elements, and there exists one column which has exactly two “”s. By exchanging it to the last column, we can transform it to be:

 ⎡⎢ ⎢ ⎢ ⎢⎣1⊤1⋯0¯Ci1⋯0⎤⎥ ⎥ ⎥ ⎥⎦

This means that one “1” is in the first row, and another is in the row from , let us say it’s . Since subtracting one row from another row will not change the determinant, we can let the first row subtract :

 ⎡⎢ ⎢ ⎢ ⎢⎣1⊤1⋯0¯Ci1⋯0⎤⎥ ⎥ ⎥ ⎥⎦→⎡⎢ ⎢ ⎢ ⎢⎣1⊤−¯Ci0⋯0¯Ci1⋯0⎤⎥ ⎥ ⎥ ⎥⎦

Now the last column only has a single “1” in the -th row. We can generate a smaller matrix by removing the -th row and the last column, and if has determinant in , so does .

 ⎡⎢ ⎢ ⎢ ⎢⎣1⊤−¯Ci0⋯0¯Ci1⋯0⎤⎥ ⎥ ⎥ ⎥⎦→S′=⎡⎢⎣1⊤−¯Ci⋯⋯⎤⎥⎦

If , then there are some positions (including -th column) in the first row will become zeros. For any column of matrix which has “0” element in the first row, there are two cases:

1. This column only contains zeros, i.e. has zero determinant.

2. This column contains a single “1”, we can generate a smaller matrix by removing this column and the row where this “1” sits. If has determinant in , so does .

Notice that it is impossible for the case that such column has two “1”s. (Since each column can have at most three “1”s, and we already remove the “1” in the first row by subtraction, and discard another “1” by removing .) In the above case (b), we can repeat removing columns and rows until we get a degenerate matrix (has 0 determinant), or a matrix whose first row does not have zeros. For the later situation, we can process it by same procedures as the original matrix unless it only has one row and one column, i.e. a matrix having single element “1” (has determinant 1).

If , we can also process it by same procedures as the original matrix .

Therefore, we have proved that any square submatrix in has determinant in , which means is TU, and hence is TU. ∎

Applying Lemma 4 and Lemma 6, we complete the proof of Theorem 2.

## 3 Proof of Theorem 3

###### Lemma 7.

Formulate the feasibility problem as problem (7), let be the objective of the formulation in Theorem 3. If , there exists a such that

 f(z)−f∗≥λ2∥z−Pz∗(z)∥2,∀z∈Ω

where is the optimal point which is closet to .

###### Proof.

Since , there exists at least an such that

 Az∗−a≤0 Bz∗=b Cz∗≤c

From Hoffman’s Theorem (Hoffman, 2003), we know that there exists a , such that

 λ2∥z−Pz∗(z)∥2≤∥[Az−a]+∥2+∥Bz−b∥2+∥[Cz−c]+∥2

Therefore, we know for any in ,

 Cz≤c

and

 λ2∥z−Pz∗(z)∥2≤∥[Az−a]+∥2+∥Bz−b∥2

Using the lemma above, we can prove the Theorem 3 now.

###### Proof.

Denote by . We have

 Δt+1 =∥zt+1−Pz∗(zt+1)∥2 ≤∥zt+1−Pz∗(zt)∥2 ≤∥zt−Pz∗(zt)+zt+1−zt∥2 =∥zt−Pz∗(zt)∥2+∥zt+1−zt∥2+2⟨zt+1−zt,zt−Pz∗(zt)⟩ =Δt−∥zt+1−zt∥2+2⟨zt+1−zt,zt+1−Pz∗(zt)⟩ ≤Δt−∥zt+1−zt∥2−2γ⟨∇f(zt),zt+1−Pz∗(zt)⟩

Let . Then we have

 T =⟨∇f(zt),zt+1−Pz∗(zt)⟩ =⟨∇f(zt),zt−Pz∗(zt)⟩+⟨∇f(zt),zt+1−zt⟩ ≥−f∗+f(zt)+f(zt+1)−f(zt)−L2∥zt+1−zt∥2 =f(zt+1)−f∗−L2∥zt+1−zt∥2

where is the Lipschitz continuous gradient constant. Back to the original inequality, we have

 ∥zt+1−Pz∗(zt+1)∥2 ≤Δt−∥zt+1−zt∥2−2γ⟨∇f(zt),zt+1−Pz∗(zt)⟩ ≤Δt−(1−Lγ)∥zt+1−zt∥2−2γ(f(zt+1)−f∗) ≤Δt−(1−Lγ)∥zt+1−zt∥2−2γλ2∥zt+1−Pz∗(zt+1)∥2

where the last inequality comes from Lemma 7.

Let and we have

 (1+λL)Δt+1 ≤Δt Δt+1 ≤11+λLΔt

which shows the linear convergence rate , then it completes the proof. ∎

## 4 Formulation of the Expected Accuracy in Crowdsourcing Task Assignment

In crowdsourcing task assignment problem, recall the objective function of problem (8):

 1mm∑j=1Eacc(Q⋅,j,X⋅,j)

For the -th task, is defined in the following:

 Eacc(Q⋅,j,X⋅,j)= P(^yj=1|yj=1)P(yj=1)+P(^yj=0|yj=0)P(yj=0) = E^YΩj,j|yj=1[I(^yj=1)]P(yj=1)+E^YΩj,j|yj=0[I(^yj=0)]P(yj=0) (12)

where is the indicator function. We can further specify this formulation by considering the equivalent forms for and :

 ^yj=1 ⇔ P(yj=1|^YΩj,j)≥P(yj=0|^YΩj,j) ⇔ P(^YΩj,j|yj=1)P(^YΩj,j|yj=0)≥P(yj=0)P(yj=1) ⇔ ∏i∈Ωj(Qij1−Qij)2^Yi,j−1≥P(yj=0)P(yj=1) ⇔ ∑i∈Ωj(2^Yi,j−1)log(Qij1−Qij)≥log(P(yj=0)P(yj=1)) ⇔ n∑i=1Xij(2^Yi,j−1)log(Qij1−Qij)−log(P(yj=0)P(yj=1))≥0

Similar derivation can be applied to (change “” to “”). Here we substitute the indicator function

to obtain a smooth approximation. Denote by and for short. The (smooth) objective turns out to be: