1 Introduction
Multitask learning is widely used learning framework in which similar tasks are considered jointly for the purpose of improving performance compared to learning the tasks separately [13]. By transferring information between related tasks it is hoped that samples will be better utilized, leading to improved generalization performance. Multitask learning has been successfully applied, for example, in natural language understanding [15], speech recognition [32], remote sensing [44], image classification [25], spam filtering [43], web search [14], disease prediction [49], and eQTL mapping [23] among other applications.
Here, we study multitask learning in a distributed setting, where each task is handled by a different machine and communication between machines is expensive. That is, each machine has access to data for a different task and needs to learn a predictor for that task, where machines communicate with each other in order to leverage the relationship between the tasks. This situation lies between a homogeneous distributed learning setting [e.g. 36], where all machines have data from the same source distribution, and inhomogeneous consensus problems [e.g. 30, 11, 6] where the goal is to reach a single consensus predictor or iterate which is the same on all machines. The main argument for this setting is that if each machine indeed has access to different data (e.g. from a different geographical region or different types of users), as in the consensus problems studied by Balcan et al. [6], then we should allow a different predictor for each distribution, instead of insisting on a single consensus predictor, while still trying to leverage the relationship and similarity between data distributions, as in classical multitask learning. As was recently pointed out by Wang et al. [41], allowing separate predictors for each task instead of insisting on a consensus predictor changes the fundamental nature of the distributed learning problem, allows for different optimization methods, and necessitates a different analysis approach, more similar to homogeneous distributed learning as studied by Shamir and Srebro [36].
The success of multitask learning relies on the relatedness between tasks. While Wang et al. [41] studied tasks related through shared sparsity, here we turn to a more general, powerful and empirically more successful model of relatedness, where the predictors for different tasks lie in some (apriori unknown) shared lowdimensional subspace and so the matrix of predictors is of low rank [3, 2, 45, 4]
. In a shared sparsity model, information from all tasks is used to learn a subset of the input features which are then used by all tasks. In contrast, in a shared subspace model, novel features, which are linear functions of the input features, are learned. The model can thus be viewed as a twolayer neural network, with the bottom layer learned jointly across tasks and the top layer taskspecific. Being arguably the most complex multilayer network that we can fully analyze, studying such models can also serve as a gateway to using deeper networks for learning shared representations.
Multitask learning with a shared subspace is wellstudied in a centralized setting, where data for all tasks are on the same machine, and some global centralized procedure is used to find a good predictor for each task. In such a situation, nuclear norm regularization is often used to leverage the low rank structure [e.g. 4, 2] and learning guarantees are known ([28] and see also Section 2). With the growth of modern massive data sets, where tasks and data often too big to handle on a single machine, it is important to develop methods also for the distributed setting. Unfortunately, the distributed multitask setting is largely unexplored and we are not aware of any prior for on distributed multitask learning with shared subspaces.
In this paper we focus on methods with efficient communication complexity (i.e. with as small as possible communication between machines), that can still leverage most of the statistical benefit of sharedsubspace multitask learning. Although all our methods are also computationally tractable and can be implemented efficiently, we are less concerned here with minimizing the runtime on each machine separately, considering communication, instead, as the main bottleneck and the main resource to be minimized [8]. This is similar to the focus in distributed optimization approaches such as ADMM [11] and DANE [37] where optimization within each machine is taken as an atomic step.
Approach  Samples  Rounds  Communication  Worker Comp.  Master Comp. 

Local  ERM  
Centralize  Nuclear Norm Minimization  
ProxGD  Gradient Comp.  SV Shrinkage  
AccProxGD  Gradient Comp.  SV Shrinkage  
ADMM  ERM  SV Shrinkage  
DFW  Gradient Comp.  Leading SV Comp.  
DGSP  ERM  Leading SV Comp.  
DNSP  ERM  Leading SV Comp. 
Summary of resources required by different approaches to distributed multitask learning with shared representations (for squared loss), in units of vector operations/communications, ignoring logfactors.
Contribution
The main contributions of this article are:

Present and formalize the sharedsubspace multitask learning [4] in the novel distributed multitask setting, identifying the relevant problems and possible approaches. We analyze two baselines, several representative firstorder distributed optimization methods, with careful sample and communication complexity analysis.

We proposed and analyzed two subspace pursuit approaches which learns the shared representation in a greedy fashion, which leverage the lowdimensional predictive structure in a communication efficient way.

We conducted comprehensive experimental comparisons of the discussed approaches on both simulated and real datasets, where we demonstrated that the proposed approaches are more communication efficient than firstorder convex optimization methods.
Table 1 summarized the approaches studied in this paper, which will be discussed in detail in the following sections.
Homogeneous, Inhomogeneous and MultiTask Distributed Learning.
We briefly review the relationship between homogeneous, inhomogeneous and multitask learning, as recently presented by Wang et al. [41].
A typical situation considered in the literature is one in which data on different machines are all drawn i.i.d from the same source distribution. In this setting, tasks on different machines are all the same, which should be taken advantage of in optimization [37]. Furthermore, as each machine has access to samples from the source distribution it can perform computations locally, without ever communicating with other machines. While having zero communication cost, this approach does not compare favorably with the centralized approach, in which all data are communicated to the central machine and used to obtain one predictor, when measured in terms of statistical efficiency. The goal in this setting is to obtain performance close to that of the centralized approach, using the same number of samples, but with low communication and computation costs [36, 20, 48, 47, 26]. Another setting considered in the distributed optimization literature is that of consensus optimization. Here each machine has data from a different distribution and the goal is to find one vector of coefficients that is good for all the separate learning or optimization problems [11, 30, 6]. The difficulty of consensus problems is that the local objectives might be rather different, and, as a result, one can obtain lower bounds on the amount of communication that must be exchanged in order to reach a joint optimum.
In this paper we suggest a novel setting that combines aspects of the above two settings. On one hand, we assume that each machine has a different source distributions , corresponding to a different task, as in consensus problems. For example, each machine serves a different geographical location, or each is at a different hospital or school with different characteristics. But if indeed there are differences between the source distributions, it is natural to learn different predictors for each machine, so that is good for the distribution typical to that machine. In this regard, our distributed multitask learning problem is more similar to singlesource problems, in that machines could potentially learn on their own given enough samples and enough time. Furthermore, availability of other machines just makes the problem easier by allowing transfer between the machine, thus reducing the sample complexity and potentially runtime. The goal, then, is to leverage as much transfer as possible, while limiting communication and runtime. As with singlesource problems, we compare our method to the two baselines, where we would like to be much better than the local approach, achieving performance nearly as good as the centralized approach, but with minimal communication and efficient runtime.
2 Setting, Formulation and Baselines
We consider a setting with tasks, each characterized by a source distribution over feature vectors and associated labels , and out goal is to find linear predictors minimizing the overall expected loss (risk) across tasks:
(2.1) 
where for convenience we denote for the matrix with columns , and
is some specified instantaneous loss function.
In the learning setting, we cannot observe directly and only have access to i.i.d. sample from each distribution , . For simplicity of presentation, we will assume that , , throughout the paper. We will denote the empirical loss where
is the local (pertask) empirical loss.
We consider a distributed setting, where each task is handled on one of separate machines, and each machine has access only to the samples drawn from . Communication between the machines is by sending realvalued vectors. Our methods work either in a broadcast communication setting, where at each iteration each machine sends a vector which is received by all other machines, or in a masteratthecenter topology where each machine sends a vector to the master node, whom in turn performs some computation and broadcasts some other vectors to all machines. Either way, we count to total number of vectors communicated.
As in standard agnosticPAC type analysis, our goal will be to obtain expected loss which is not much larger then the expected loss of some (unknown) reference predictor^{1}^{1}1Despite the notation, need not be the minimizer of the expected loss. We can think of it as the minimizer inside some restricted hypothesis class, though all analysis and statements hold for any chosen reference predictor . , and we will measure the excess error over this goal. To allow obtaining such guarantees we will assume: The loss function is Lipschitz and bounded^{2}^{2}2
This is only required for the high probability bounds.
by , be twice differentiable and smooth, that isAll the data points are bounded by unit length, i.e.
and the reference predictors have bounded norm:
for some .
The simplest approach, which we refer to as Local, is to learn a linear predictor on each machine independently of other machines. This single task learning approach ignores the fact that the tasks are related and that sharing information between them could improve statistical performance. However, the communication cost for this procedure is zero, and with enough samples it can still drive the excess error to zero. However, compared to procedures discussed later, sample complexity (number of samples required to achieve small excess error) is larger. A standard Rademacher complexity argument [7] gives the following generalization guarantee, which is an extension of Theorem 26.12 in ShalevShwartz and BenDavid [33]. Suppose Assumption 1 holds. Then with probability at least ,
where with . That is, in order to ensure excess error, we need
samples from each task.
At the other extreme, if we ignore all communication costs, and, e.g. communicate all data to a single machine, we can significantly leverage the shared subspace. To understand this, we will first need to introduce two assumptions: one about the existence of a shared subspace (i.e. that the reference predictor is indeed lowrank), and the other about the spread of the data:
There is a constant , such that
Since the data is bounded, we always have , with being a measure of how spread out the data is in different direction. A value of indicates the data is entirely contained in a onedimensional line. In this case, the predictor matrix will also always be rankone, imposing a lowrank structure is meaningless and we can’t expect to gain from it. However, when is close to , or at least high, the data is spread in many directions and the lowrank assumption is meaningful. We can think of as the “effective dimensionality” of the data, and hope to gain when .
With these two assumptions in hand, we can think of minimizing the empirical error subject to a rank constraint on . This is a hard and nonconvex optimization task, but we can instead use the nuclear norm (aka tracenorm) as a convex surrogate for the rank. This is because if Assumptions 1 and 2 hold, then we also have:
(2.2) 
With this in mind, we can define the following centralized predictor:
(2.3) 
which achieves the improved excess error guarantee: (Theorem 1 in Maurer and Pontil [28]) Suppose Assumptions 1, 2 and 2 hold. Then with probability at least ,
The sample complexity per task, up to logarithmic factors, is thus only:
When , this is a reduction by a factor of . That is, it is as if we needed to only learn linear predictors instead of .
The problem is that a naive computation of requires collecting all data on a single machine, i.e. communicating samples per machine. In the next Sections, we aim at developing methods of approximating using communication efficient methods, or computing an alternate predictor with similar statistical properties but using much less communication.
3 Distributed Convex Optimization
In this section, we study how to obtain the sharing benefit of the centralized approach using distributed convex optimization techniques, while keeping the communication requirements at low.
To enjoy the benefit of nuclearnorm regularization while avoid heavy communication cost of Centralize, a flexible strategy is to solve the convex objective (2.3) via distributed optimization techniques. Let be the solution at iteration for some iterative distributed optimization algorithm for the following constrained objective:
(3.1) 
By the generalization error decompsition [10],
Suppose satisfying with . Then will have the generalization error of order . Therefore in order to study the generalization performance, we will study how the optimization error decreases as the function of the number of iterations .
Constrained vs Regularized Objective
Note that the constrained objective (3.1) is equivalent to the following regularized objective with a proper choice of :
(3.2) 
Though they are equivalent, specific optimization algorithms might sometimes be more suitable for one particular type of objectives ^{3}^{3}3e.g. ADMM for regularized objective and FrankWolfe for constrained objective. Gradient descent methods can be adopted for both, leads to proximal and projected methods, respectively.. For convenience in the following discussion we didn’t distinguish between these two formulations.
3.1 Distributed Proximal Gradient
Maybe the simplest distributed optimization algorithm for (3.2) is the proximal gradient descent. It is not hard to see that computation of the gradient can be easily done in a distributed way as the losses are decomposable across machines:
where
Thus each machine needs to compute the gradient on the local dataset and send it to the master. The master concatenates the gradient vectors to form the gradient matrix . Finally, the master computes the proximal step
(3.3) 
which has the following closed form solution [12]: let be the SVD of , then with applied elementwise.
3.2 Distributed Accelerated Gradient
It is also possible to use Nesterov’s acceleration idea [29] to improve the convergence of the proximal gradient algorithm from to [22]. Using the distributed accelerated proximal gradient descent, one needs rounds of communication with a total bits communicated per machine to achieve generalization error. The algorithm is summarized in Algorithm 5 (in Appendix), where the master maintains two sequences: and . First, a proximal gradient update of is done based on
(3.4) 
and then is updated based on a combination of the current and the difference with previous
(3.5) 
ADMM and DFW
4 Greedy Representation Learning
In this section we propose two distributed algorithms which select the subspaces in a greedy fashion, instead of solving the nuclear norm regularized convex program.
4.1 Distributed Greedy Subspace Pursuit
Our greedy approach is inspired by the methods used for sparse signal reconstruction [39, 34]. Under the assumption that the optimal model is lowrank, say rank , we can write as a sum of rank matrices:
where , and . In the proposed approach, the projection matrix is learned in a greedy fashion. At every iteration, a new onedimensional subspace is identified that leads to an improvement in the objective. This subspace is then included into the existing projection matrix. Using the new expanded projection matrix as the current feature representation, we refit the model to obtain the coefficient vectors . In the distributed setting, there is a master that gathers local gradient information from each task. Based on this information, it then computes the subspace to be added to the projection matrix and sends it to each machine. The key step in the distributed greedy subspace pursuit algorithm is the addition of the subspace. One possible choice is the principle component of the gradient direction; after the master collected the gradient matrix , it computes the top left and right singular vectors of . Let be the largest singular vectors of . The left singular vector is used as a new subspace to be added to the projection matrix . This vector is sent to each machine, which then concatenate it to the projection matrix and refit the model with the new representation. Algorithm 1 details the steps.
Distributed gradient subspace pursuit (DGSP), detailed in Algorithm 1, creates subspaces that are orthogonal to each other, as shown in the following proposition which is proved in Appendix D: At every iteration of Algorithm 1, the columns of are orthonormal.
Both the distributed gradient subspace pursuit and the distributed FrankWolfe use the leading singular vector of the gradient matrix iteratively. Moreover, leading singular vectors of the gradient matrix have been used in greedy selection procedures for solving lowrank matrix learning problems [35, 42]. However, DGSP utilize the learned subspace in a very different way: GECO [35] refit the lowrank matrix under a larger subspace which is spanned by all left and right singular vectors; while OR1MP [42] only adjust the linear combination parameters of the rank1 matrices. The DGSP algorithm do not restrict on the joint subspaces , but focused on the lowdimensional subspace induced the projection matrix
, and estimate the task specific predictors
based on the learned representation.Next, we present convergence guarantees for the distributed gradient subspace pursuit. First, note that the smoothness of implies the smoothness property for any rank1 update. Suppose Assumption 1 holds. Then for any and unit length vectors and , we have
We defer the proof in Appendix E. The following theorem states the number of iterations needed for the distributed gradient subspace pursuit to find an suboptimal solution. Suppose Assumption 1 holds. Then the distributed gradient subspace pursuit finds such that when
We defer the proof in Appendix F. Theorem 4.1 tells us that for the distributed gradient subspace pursuit requires iterations to reach accuracy. Since each iteration requires communicating number, the communication cost per machine is . In some applications this communication cost might be still too high and in order to improve it we will try to reduce the number of rounds of communication. To that end, we develop a procedure that utilizes the secondorder information to improve the convergence. Algorithm 6 describes the Distributed Newton Subspace Pursuit algorithm (DNSP). Note that distributed optimization with secondorder information have been studied recently to achieve communication efficiency [37, 46].
Compared to the gradient based methods, the DNSP algorithm uses secondorder information to find subspaces to work with. At each iteration, each machine computes the Newton direction
based on the current solution and sends it to the master. The master computes the overall Newton direction by concatenating the Newton direction for each task
and computes the top singular vectors of . The top left singular vector is is sent back to every machine, which is then concatenated to the current projection matrix. Each machine refits the predictors using the new representation. Note that at every iteration a GramSchmidt step is performed to ensure that the learned basis are orthonormal.
DNSP is a Newtonlike method which uses secondorder information, thus its generic analysis is not immediately apparent. However empirical results in the next section illustrate good performance of the proposed DNSP.
5 Experiments
We first illustrate performance of different procedures on simulated data. We generate data according to
for regression problems and
for classification problems. We generate the lowrank as follows. We first generate two matrices
with entries sampled independently from a standard normal distribution. Then we extract the left and right singular vectors of
, denoted as . Finally, we set , where is a diagonal matrix with exponentially decaying entries: . The feature vectors are sampled from a mean zero multivariate normal with the covariance matrix , . The regularization parameters for all approaches were optimized to give the best prediction performance over a heldout validation dataset. For ProxGD and AccProxGD, we initialized the solution from Local. Our simulation results are averaged over 10 independent runs.We investigate how the performance of various procedures changes as a function of problem parameters . We compare the following procedures: i) Local
, where each machine solves an empirical risk minimization problem (ordinary least squares or logistic regression) . ii) Nuclearnorm regularization: which is a popular
Centralize approach: all machines send their data to the master, the master solves a nuclearnorm regularized loss minimization problem. iii) Learning with the best representation (BestRep): which assumes the true projection matrix is known, and just fit ordinal least squares or logistic regression model in the projected lowdimensional subspace . Note that this is not a practical approach since in practice we do not know the best lowdimensional representations of the data. iv) Convex optimization approach which runs distributed optimization algorithms over the nuclear normregularized objective: here we implemented and compared the following algorithms: distributed proximal gradient (ProxGD); distributed accelerated proximal gradient, (AccProxGD); distributed alternating direction method of multipliers (ADMM); distributed FrankWolfe (DFW) . v) The proposed DGSP and DNSP approaches. The simulation results for regression and classification problems are shown in Figure 1 and 2^{4}^{4}4For better visualization, here we omit the plot for DFW as its performance is significantly worse than others., respectively. We plot how the excess prediction error decreases as the number of rounds of communications increases (Local, Centralize and BestRep are one shot approaches thus the lines are horizontal). From the plots, we have the following observations:
Nuclear norm regularization boosts the prediction performance over plain single task learning significantly, which shows clear advantage of leveraging the shared representation in multitask learning.

ADMM and AccProxGD perform reasonably well , especially ADMM. One reason for the effectiveness of ADMM is that for the problem of nuclear norm regularized multitask learning considered here, the ADMM update solves regularized ERM problems at every iteration. ADMM and AccProxGD clearly outperform ProxGD.

ProxGD and DGSP perform similarly. DGSP usually becomes worse as the iterations increases , while ProxGD converges to a global optimum of the nuclear norm regularized objective.

DNSP is the most communicationefficient method, and usually converges to a solution that is slightly better compared to the optimum of the nuclear regularization. This shows that secondorder information helps a lot in reducing the communication cost.

The DFW performs the worst in most cases, even though DFW shares some similarity with DGSP in learning the subspace. The empirical results suggest the refitting step in DGSP is very important.
Oneshot SVD truncation
A natural question to ask is whether there exists a oneshot communication method for the shared representation problem considered here, that still matches the performance of centralized methods. One reasonable solution is to consider the following SVD truncation approach, which is based on the following derivation: consider the following well specified linear regression model:
where is drawn from meanzero Gaussian noise. It is easy to verify the following equation for OLS estimation:
Since is just plus some meanzero Gaussian noise, it is natural to consider the following lowrank matrix denoising estimator:
where the solution is a simple SVD truncation, and can be implemented in a oneshot way: each worker send its Local solution to the master, which then performs an SVD truncation step to maintain the top components
and send the resulting estimation back to each worker, where are top components of . Though this approach might work well for some simple scenarios, but will generally fail when the features are highly correlated: although the Local solution can output normal estimation of , the estimation noise might be highly correlated (depend on the correlation between features), which makes the SVD truncation estimation not reliable. To illustrate this, consider a more complex simulation which follows the same setup as above setting, except that now the feature vectors are sampled from a higher correlation matrix , . The regression simulation results are shown in Figure 3, where we see that the oneshot SVD truncation approach does not significantly outperforms Local, sometimes even slightly worse.
Besides simulation, we also conducted extensive experiments on real world datasets, which are presented in Appendix H due to space limitation.
6 Conclusion
We studied the problem of distributed representation learning for multiple tasks, discussed the implementation and guarantees for distributed convex optimization methods, and presented two novel algorithms to learn lowdimensional projection in a greedy way, which can be communication more efficient than distributed convex optimization approaches. All approaches are extensively evaluated on simulation and real world datasets.
References

Agarwal et al. (2012)
A. Agarwal, S. Negahban, and M. J. Wainwright.
Fast global convergence of gradient methods for highdimensional statistical recovery.
Ann. Stat., 40(5):2452–2482, 2012.  Amit et al. (2007) Y. Amit, M. Fink, N. Srebro, and S. Ullman. Uncovering shared structures in multiclass classification. In ICML, pages 17–24. ACM, 2007.
 Ando and Zhang (2005) R. K. Ando and T. Zhang. A framework for learning predictive structures from multiple tasks and unlabeled data. J. Mach. Learn. Res., 6:1817–1853, 2005.
 Argyriou et al. (2008) A. Argyriou, T. Evgeniou, and M. Pontil. Convex multitask feature learning. Mach. Learn., 73(3):243–272, 2008.
 Bach et al. (2011) F. Bach, R. Jenatton, J. Mairal, and G. Obozinski. Optimization with sparsityinducing penalties. Found. Trends Mach. Learn., 4(1):1–106, 2011.
 Balcan et al. (2012) M.F. Balcan, A. Blum, S. Fine, and Y. Mansour. Distributed learning, communication complexity and privacy. In JMLR W&CP 23: COLT 2012, volume 23, pages 26.1–26.22, 2012.
 Bartlett and Mendelson (2002) P. L. Bartlett and S. Mendelson. Rademacher and gaussian complexities: Risk bounds and structural results. J. Mach. Learn. Res., 3:463–482, 2002.

Bekkerman et al. (2011)
R. Bekkerman, M. Bilenko, and J. Langford.
Scaling up machine learning: Parallel and distributed approaches
. Cambridge University Press, 2011.  Bellet et al. (2015) A. Bellet, Y. Liang, A. B. Garakani, M.F. Balcan, and F. Sha. A distributed frankwolfe algorithm for communicationefficient sparse learning. In SDM, pages 478–486. 2015.
 Bousquet and Bottou (2008) O. Bousquet and L. Bottou. The tradeoffs of large scale learning. In NIPS, pages 161–168, 2008.
 Boyd et al. (2011) S. P. Boyd, N. Parikh, E. Chu, B. Peleato, and J. Eckstein. Distributed optimization and statistical learning via the alternating direction method of multipliers. Found. Trends Mach. Learn., 3(1):1–122, 2011.

Cai et al. (2010)
J.F. Cai, E. J. Candès, and Z. Shen.
A singular value thresholding algorithm for matrix completion.
SIAM Journal on Optimization, 20(4):1956–1982, 2010.  Caruana (1997) R. Caruana. Multitask learning. Mach. Learn., 28(1):41–75, 1997.
 Chapelle et al. (2010) O. Chapelle, P. Shivaswamy, S. Vadrevu, K. Weinberger, Y. Zhang, and B. Tseng. Multitask learning for boosting with application to web search ranking. In KDD, pages 1189–1198. ACM, 2010.
 Collobert et al. (2011) R. Collobert, J. Weston, L. Bottou, M. Karlen, K. Kavukcuoglu, and P. Kuksa. Natural language processing (almost) from scratch. J. Mach. Learn. Res., 12:2493–2537, 2011.
 Frank and Wolfe (1956) M. Frank and P. Wolfe. An algorithm for quadratic programming. Naval research logistics quarterly, 3(12):95–110, 1956.
 He and Yuan (2012) B. He and X. Yuan. On the convergence rate of the douglasrachford alternating direction method. SIAM Journal on Numerical Analysis, 50(2):700–709, 2012.
 Hong and Luo (2012) M. Hong and Z.Q. Luo. On the linear convergence of the alternating direction method of multipliers. ArXiv eprints, arXiv:1208.3922, 2012.
 Jaggi (2013) M. Jaggi. Revisiting frankwolfe: Projectionfree sparse convex optimization. In ICML, pages 427–435, 2013.
 Jaggi et al. (2014) M. Jaggi, V. Smith, M. Takác, J. Terhorst, S. Krishnan, T. Hofmann, and M. I. Jordan. Communicationefficient distributed dual coordinate ascent. In NIPS, pages 3068–3076, 2014.
 Jain et al. (2013) P. Jain, P. Netrapalli, and S. Sanghavi. Lowrank matrix completion using alternating minimization. In STOC, pages 665–674. ACM, 2013.
 Ji and Ye (2009) S. Ji and J. Ye. An accelerated gradient method for trace norm minimization. In ICML, pages 457–464. ACM, 2009.
 Kim and Xing (2010) S. Kim and E. P. Xing. Treeguided group lasso for multitask regression with structured sparsity. In ICML, pages 543–550, 2010.
 LacosteJulien and Jaggi (2015) S. LacosteJulien and M. Jaggi. On the global linear convergence of frankwolfe optimization variants. In NIPS, pages 496–504, 2015.

Lapin et al. (2014)
M. Lapin, B. Schiele, and M. Hein.
Scalable multitask representation learning for scene classification.
In CVPR, pages 1434–1441, 2014.  Lee et al. (2015) J. D. Lee, Y. Sun, Q. Liu, and J. E. Taylor. Communicationefficient sparse regression: a oneshot approach. ArXiv eprints, arXiv:1503.04337, 2015.
 Lenk et al. (1996) P. J. Lenk, W. S. DeSarbo, P. E. Green, and M. R. Young. Hierarchical bayes conjoint analysis: Recovery of partworth heterogeneity from reduced experimental designs. Marketing Science, 15(2):173–191, 1996.
 Maurer and Pontil (2013) A. Maurer and M. Pontil. Excess risk bounds for multitask learning with trace norm regularization. pages 55–76, 2013.
 Nesterov (1983) Y. Nesterov. A method of solving a convex programming problem with convergence rate . In Soviet Mathematics Doklady, volume 27, pages 372–376, 1983.
 Ram et al. (2010) S. S. Ram, A. Nedić, and V. V. Veeravalli. Distributed stochastic subgradient projection algorithms for convex optimization. Journal of optimization theory and applications, 147(3):516–545, 2010.
 Sander and Schneider (1991) C. Sander and R. Schneider. Database of homologyderived protein structures and the structural meaning of sequence alignment. Proteins: Structure, Function, and Bioinformatics, pages 56–68, 1991.
 Seltzer and Droppo (2013) M. L. Seltzer and J. Droppo. Multitask learning in deep neural networks for improved phoneme recognition. In ICASSP, pages 6965–6969. IEEE, 2013.
 ShalevShwartz and BenDavid (2014) S. ShalevShwartz and S. BenDavid. Understanding machine learning: From theory to algorithms. Cambridge University Press, 2014.
 ShalevShwartz et al. (2010) S. ShalevShwartz, N. Srebro, and T. Zhang. Trading accuracy for sparsity in optimization problems with sparsity constraints. SIAM Journal on Optimization, 20(6):2807–2832, 2010.
 ShalevShwartz et al. (2011) S. ShalevShwartz, A. Gonen, and O. Shamir. Largescale convex minimization with a lowrank constraint. In ICML, 2011.
 Shamir and Srebro (2014) O. Shamir and N. Srebro. Distributed stochastic optimization and learning. In Allerton, pages 850–857. IEEE, 2014.
 Shamir et al. (2014) O. Shamir, N. Srebro, and T. Zhang. Communication efficient distributed optimization using an approximate newtontype method. In ICML, pages 1000–1008, 2014.
 SpyromitrosXioufis et al. (2012) E. SpyromitrosXioufis, G. Tsoumakas, W. Groves, and I. Vlahavas. Multitarget regression via input space expansion: Treating targets as inputs. ArXiv eprints, arXiv:1211.6581, 2012.
 Tropp (2004) J. A. Tropp. Greed is good: Algorithmic results for sparse approximation. IEEEit, 50(10):2231–2242, 2004.
 Turnbull et al. (2008) D. Turnbull, L. Barrington, D. Torres, and G. Lanckriet. Semantic annotation and retrieval of music and sound effects. IEEE Transactions on Acoustics, Speech and Signal Processing, 16(2):467–476, 2008.
 Wang et al. (2015a) J. Wang, M. Kolar, and N. Srebro. Distributed multitask learning. ArXiv eprints, arXiv:1510.00633, 2015a.
 Wang et al. (2015b) Z. Wang, M.J. Lai, Z. Lu, W. Fan, H. Davulcu, and J. Ye. Orthogonal rankone matrix pursuit for low rank matrix completion. SIAM Journal on Scientific Computing, 37(1):A488–A514, 2015b.
 Weinberger et al. (2009) K. Weinberger, A. Dasgupta, J. Langford, A. Smola, and J. Attenberg. Feature hashing for large scale multitask learning. In ICML, pages 1113–1120. ACM, 2009.
 Xue et al. (2007) Y. Xue, X. Liao, L. Carin, and B. Krishnapuram. Multitask learning for classification with dirichlet process priors. J. Mach. Learn. Res., 8:35–63, 2007.
 Yuan et al. (2007) M. Yuan, A. Ekici, Z. Lu, and R. Monteiro. Dimension reduction and coefficient estimation in multivariate linear regression. J. R. Stat. Soc. B, 69(3):329–346, 2007.
 Zhang and Xiao (2015) Y. Zhang and L. Xiao. Communicationefficient distributed optimization of selfconcordant empirical loss. ArXiv eprints, arXiv:1501.00263, 2015.
 Zhang et al. (2012) Y. Zhang, M. J. Wainwright, and J. C. Duchi. Communicationefficient algorithms for statistical optimization. In NIPS, pages 1502–1510, 2012.
 Zhang et al. (2013) Y. Zhang, J. C. Duchi, M. I. Jordan, and M. J. Wainwright. Informationtheoretic lower bounds for distributed statistical estimation with communication constraints. In NIPS, pages 2328–2336, 2013.
 Zhou et al. (2013) J. Zhou, J. Liu, V. A. Narayan, and J. Ye. Modeling disease progression via multitask learning. NeuroImage, 78:233 – 248, 2013.
Appendix
Appendix A Distributed Alternating Direction Methods of Multipliers
The Alternating Direction Methods of Multipliers (ADMM) is also a popular method for distributed optimization [11] and can be used to solve the distributed lowrank multitask learning problem. We first write the objective (2.3) as
By introducing the Lagrangian and augmented terms, we get the following unconstrained problem:
where is a parameter controlling the augmentation level. Note that except for , the augmented Lagrangian objective are decomposable across tasks. To implement the distributed ADMM algorithm, we let the workers maintain the data and , while the master maintains and . At round , each machine separately solves
(A.1) 
which is minimizing the local loss plus a regularization term. Next, each worker sends their solution to the master, which performs the following updates for and
(A.2)  
(A.3) 
which have closedform solutions.
The algorithm ADMM is summarized in Algorithm 2. Note that compared to methods discussed before, ADMM needs to communicate three dimensional vectors between each worker and the master at each round, while the proximal gradient approaches only communicate two dimensional vectors per round. Based on convergence results of ADMM [17], rounds of communication are needed to obtain generalization error.
Appendix B Distributed FrankWolfe Method
Another approach we consider is the distributed FrankWolfe method [16, 19, 9]. This methods does not require performing SVD, which might bring additional computational advantages. Instead of directly minimizing the nuclear norm regularized objective, the FrankWolfe algorithm considers the equivalent constrained minimization problem
At each step, FrankWolfe algorithm considers the following direction to update
where is the leading singular vectors of . The next iterate is obtained as
where is a step size parameter. To implement this algorithm in a distributed way, the master first collects the gradient matrix and computes and . The vector is sent to th machine, which performs the following update:
(B.1) 
The algorithm is summarized in Algorithm 3. Similar to the distributed (accelerated) proximal gradient descent, the distributed FrankWolfe only requires communication of two dimensional vectors per round. Though computationally cheaper compared to other methods considered in this section, the distributed FrankWolfe algorithm enjoys similar convergence guarantees to the distributed proximal gradient descent [19], that is, after iterations, the solution will be suboptimal.
Appendix C Pseudocode of the algorithms
Appendix D Proof of Proposition 4.1
Proof.
It is sufficient to prove that at every iteration, the current projection matrix and the subspace to be added are orthogonal to each other. Note that by the optimality condition:
Since is the leading left singular vector of , we have . Each column of has unit length, since it is a left singular vector of some matrix. ∎
Appendix E Proof of Proposition 4.1
Proof.
It is sufficient to prove that the largest eigenvalue of
does not exceed . Since is a block diagonal matrix, it is sufficient to show that for every block , the largest eigenvalue of the block is not larger than .This is true by the smoothness of and the fact that the data points have bounded length:
∎
Appendix F Proof of Theorem 4.1
Proof.
By the smoothness of , we know
(F.1) 
Let . Since is a minimizer of with respect to , we have and therefore . From convexity of , we have
Combining with the display above
By choosing
we have
Using Lemma G in Appendix we know that after
iterations, we have . ∎
Appendix G An auxiliary lemma
(Lemma B.2 of [34]) Let and let be a sequence such that for all . Let be a positive scalar and be a positive integer such that . Then .
Appendix H Evaluation on Real World Datasets
We also evaluate discussed algorithms on several real world data sets, with of the whole dataset as training set, as heldout validation, then report the testing performance on the remaining . For the real data, we have observed that adding regularization usually helps improving the generalization performance. For the Local procedure we added an
regularization term (leads to ridge regression or
regularized logistic regression). For DGSP and DNSP, we also add an regularization in finding the subspaces and refitting . We have worked on the following multitask learning datasets:School.^{5}^{5}5http://cvn.ecp.fr/personnel/andreas/code/mtl/index.html The dataset consists of examination scores of students from London’s secondary schools during the years 1985, 1986, 1987. There are 27 schoolspecific and studentspecific features to describe each student. The instances are divided by different schools, and the task is to predict the students’ performance. We only considered schools with at least records, which results in 72 tasks in total. The maximum number of records for each individual school is 260.
Computer Survey. The data is taken from a conjoint analysis experiment [27] which surveyed 180 persons about the probability of purchasing 20 kinds of personal computers. There are 14 variables for each computer, the response is an integer rating with scale .
ATP.^{6}^{6}6http://mulan.sourceforge.net/datasets.html The task here is to predict the airline ticket price [38]. We are interested in the minimum prices next day for some specific observation date and departure date pairs. Each case is described by 411 features, and there are 6 target minimum prices for different airlines to predict. The sample size is 337.
Protein. Given the amino acid sequence, we are interested predicting the protein secondary structure [31]. We tackle the problem by considering the following three binary classification tasks: coil vs helix, helix vs strand, strand vs coil. Each sequence is described by 357 features. There are 24,387 instances in total.
Landmine. The data is collected from 19 landmine detection tasks [44]
. Each landmine field is represented by a 9dimensional vector extracted from radar images, containing momentbased, correlationbased, energy ratio, and spatial variance features. The sample size for each task varies from 445 to 690.
Cal500.^{7}^{7}7http://eceweb.ucsd.edu/~gert/calab/ This music dataset [40] consists of 502 songs, where for each song 68 features are extracted. Each task is to predict whether a particular musically relevant semantic keyword should be an annotation for the song. We only consider tags with at least 50 times apperance, which results in 78 prediction tasks.
We compared various approaches as in the simulation study, except the BestRep as the best lowdimensional representation is unknown. We also compared with AltMin, which learns lowrank prediction matrix using the alternating minimization [21]. The results are shown in Figure 4. Since the labels for the real world classification datasets are often unbalanced, we report averaged area under the curve (AUC) instead of classification accuracy. We have the following observations:

The distributed firstorder approaches converge much slower than in simulations, especially on ATP and Cal500. We suspect this is because in the simulation study, the generated data are usually well conditioned, which makes faster convergence possible for such methods [1, 18]. On real data, the condition number can be much worse.

In most case, DNSP is the best in terms of communicationefficiency. DGSP also has reasonable performance with fewer round of communications compared to distributed firstorder approaches.

Among the firstorder distributed convex optimization methods, AccProxGD is overall the most communicationefficient, while DFW is the worst, though it might have some advantages in terms of computation. Also, we observed significant zigzag behavior of the DFW algorithm, as discussed in [24].
Comments
There are no comments yet.