We consider the networked multi-agent reinforcement learning (MARL) problemin a fully decentralized setting, where agents learn to coordinate to achievethe joint success. This problem is widely encountered in many areas includingtraffic control, distributed control, and smart grids. We assume that thereward function for each agent can be different and observed only locally bythe agent itself. Furthermore, each agent is located at a node of acommunication network and can exchanges information only with its neighbors.Using softmax temporal consistency and a decentralized optimization method, weobtain a principled and data-efficient iterative algorithm. In the first stepof each iteration, an agent computes its local policy and value gradients andthen updates only policy parameters. In the second step, the agent propagatesto its neighbors the messages based on its value function and then updates itsown value function. Hence we name the algorithm value propagation. We prove anon-asymptotic convergence rate 1/T with the nonlinear function approximation.To the best of our knowledge, it is the first MARL algorithm with convergenceguarantee in the control, off-policy and non-linear function approximationsetting. We empirically demonstrate the effectiveness of our approach inexperiments.
Quick Read (beta)
Value Propagation for Decentralized Networked Deep Multi-agent Reinforcement Learning
We consider the networked multi-agent reinforcement learning (MARL) problem in a fully decentralized setting, where agents learn to coordinate to achieve joint success. This problem is widely encountered in many areas including traffic control, distributed control, and smart grids. We assume each agent is located at a node of a communication network and can exchange information only with its neighbors. Using softmax temporal consistency, we derive a primal-dual decentralized optimization method and obtain a principled and data-efficient iterative algorithm named value propagation. We prove a non-asymptotic convergence rate of with nonlinear function approximation. To the best of our knowledge, it is the first MARL algorithm with a convergence guarantee in the control, off-policy, non-linear function approximation, fully decentralized setting.
Value Propagation for Decentralized Networked Deep Multi-agent Reinforcement Learning
noticebox[b]33rd Conference on Neural Information Processing Systems (NeurIPS 2019), Vancouver, Canada.\[email protected]
Multi-agent systems have applications in a wide range of areas such as robotics, traffic control, distributed control, telecommunications, and economics. For these areas, it is often difficult or simply impossible to predefine agents’ behaviour to achieve satisfactory results, and multi-agent reinforcement learning (MARL) naturally arises (Bu et al., 2008; Tan, 1993). For example, El-Tantawy et al. (2013) model a traffic signal control problem as a multi-player stochastic game and solve it with MARL. MARL generalizes reinforcement learning by considering a set of agents (decision makers) sharing a common environment. However, multi-agent reinforcement learning is a challenging problem since the agents interact with both the environment and each other. For instance, independent Q-learning—treating other agents as a part of the environment—often fails as the multi-agent setting breaks the theoretical convergence guarantee of Q-learning and makes the learning process unstable (Tan, 1993). Rashid et al. (2018); Foerster et al. (2018); Lowe et al. (2017) alleviate such a problem using a centralized network (i.e., being centralized for training, but decentralized during execution.). Its communication pattern is illustrated in the left panel of Figure 1.
Despite the great success of (partially) centralized MARL approaches, there are various scenarios, such as sensor networks (Rabbat and Nowak, 2004) and intelligent transportation systems (Adler and Blue, 2002) , where a central agent does not exist or may be too expensive to use. In addition, privacy and security are requirements of many real world problems in multi-agent system (also in many modern machine learning problems) (Abadi et al., 2016; Kurakin et al., 2016) . For instance, in Federated learning (McMahan et al., 2016), the learning task is solved by a lose federation of participating devices (agents) without the need to centrally store the data, which significantly reduces privacy and security risk by limiting the attack surface to only the device. In the agreement problem (DeGroot, 1974; Mo and Murray, 2017), a group of agents may want to reach consensus on a subject without leaking their individual goal or opinion to others. Obviously, centralized MARL violates privacy and security requirements. To this end, we and others have advocated the fully decentralized approaches, which are useful for many applications including unmanned vehicles (Fax and Murray, 2002), power grid (Callaway and Hiskens, 2011), and sensor networks (Cortes et al., 2004). For these approaches, we can use a network to model the interactions between agents (see the right panel of Figure 1). Particularly, We consider a fully cooperative setting where each agent makes its own decision based on its local reward and messages received from their neighbors. Thus each agent preserves the privacy of its own goal and policy. At the same time, through the message-passing all agents achieve consensus to maximize the averaged cumulative rewards over all agents; see Equation (3).
In this paper, we propose a new fully decentralized networked multi-agent deep reinforcement learning algorithm. Using softmax temporal consistency (Nachum et al., 2017; Dai et al., 2018) to connect value and policy updates, we derive a new two-step primal-dual decentralized reinforcement learning algorithm inspired by a primal decentralized optimization method (Hong et al., 2017) 11 1 The objective in Hong et al. (2017) is a primal optimization problem with constraint. Thus they introduce a Lagrange multiplier like method to solve it (so they call it primal-dual method ). Our objective function is a primal-dual optimization problem with constraint. . In the first step of each iteration, each agent computes its local policy, value gradients and dual gradients and then updates only policy parameters. In the second step, each agent propagates to its neighbors the messages based on its value function (and dual function) and then updates its own value function. Hence we name the algorithm value propagation. It preserves the privacy in the sense that no individual reward function is required for the network-wide collaboration. We approximate the local policy, value function and dual function of each agent by deep neural networks, which enables automatic feature generation and end-to-end learning.
Contributions:  We propose the value propagation algorithm and prove that it converges with the rate even with the non-linear deep neural network function approximation. To the best of our knowledge, it is the first deep MARL algorithm with non-asymptotic convergence guarantee. At the same time, value propagation can use off-policy updates, making it data efficient. When it reduces to the single agent case, it provides a proof of (Dai et al., 2018) in the realistic setting; see remarks of algorithm 3.3 in Section 3.3.  The objective function in our problem is a primal-dual decentralized optimization form (see (3.2)), while the objective function in (Hong et al., 2017) is a primal problem. When our method reduces to pure primal analysis, we extend (Hong et al., 2017) to the stochastic and biased gradient setting which may be of independent interest to the optimization community. In the practical implementation, we extend ADAM into the decentralized setting to accelerate training.
MDP Markov Decision Process (MDP) can be described by a 5-tuple (): is the finite state space, is the finite action space, are the transition probabilities, are the real-valued immediate rewards and is the discount factor. A policy is used to select actions in the MDP. In general, the policy is stochastic and denoted by , where is the conditional probability density at associated with the policy. Define to be the optimal value function. It is known that is the unique fixed point of the Bellman optimality operator, The optimal policy is related to by the following equation:
Softmax Temporal Consistency Nachum et al. (2017) establish a connection between value and policy based reinforcement learning based on a relationship between softmax temporal value consistency and policy optimality under entropy regularization. Particularly, the soft Bellman optimality is as follows,
where and controls the degree of regularization. When , above equation reduces to the standard Bellman optimality condition. An important property of soft Bellman optimality is the called temporal consistency, which leads to the Path Consistency Learning.
A straightforward way to apply temporal consistency is to optimize the following problem, Dai et al. (2018) get around the double sampling problem of above formulation by introduce a primal-dual form
where , controls the trade-off between bias and variance.
In the following discussion, we use to denote the Euclidean norm over the vector, stands for the transpose of , and denotes the entry-wise product between two vectors.
3 Value Propagation
In this section, we present our multi-agent reinforcement learning algorithm, i.e., value propagation. To begin with, we extend the MDP model to the Networked Multi-agent MDP model following the definition in (Zhang et al., 2018). Let be an undirected graph with agents (node). represents the set of edges. means agent and can communicate with each other through this edge. A networked multi-agent MDP is characterized by a tuple : is the global state space shared by all agents (It could be partially observed, i.e., each agent observes its own state , see our experiment). is the action space of agent , is the joint action space, is the transition probability, denotes the local reward function of agent . We assume rewards are observed only locally to preserve the privacy of the each agent’s goal. At each time step, agents observe and make the decision . Then each agent just receives its own reward , and the environment switches to the new state according to the transition probability. Furthermore, since each agent make the decisions independently, it is reasonable to assume that the policy can be factorized, i.e., (Zhang et al., 2018). We call our method fully-decentralized method, since reward is received locally, the action is executed locally by agent, critic (value function) are trained locally.
3.1 Multi-Agent Softmax Temporal Consistency
The goal of the agents is to learn a policy that maximizes the long-term reward averaged over the agent, i.e.,
In the following, we adapt the temporal consistency into the multi-agent version. Let be the optimal value function and be the corresponding policy. Apply the soft temporal consistency, we obtain that for all , is the unique pair that satisfies
A optimization problem inspired by (4) would be
There are two potential issues of above formulation: First, due to the inner conditional expectation, it would require two independent samples to obtain the unbiased estimation of gradient of (Dann et al., 2014). Second, is a global variable over the network, thus can not be updated in a decentralized way.
Change the variable , the objective function becomes
3.2 Decentralized Formulation
So far the problem is still in a centralized form, and we now turn to reformulating it in a decentralized way. We assume that policy, value function, dual variable are all in the parametric function class. Particularly, each agent’s policy is and The value function is characterized by the parameter , while represents the parameter of . Similar to (Dai et al., 2018), we optimize a slightly different version from (6).
where controls the bias and variance trade-off. When , it reduces to the pure primal form.
We now consider the second issue that is a global variable. To address this problem, we introduce the local copy of the value function, i.e., for each agent . In the algorithm, we have a consensus update step, such that these local copies are the same, i.e., , or equivalently , where are parameter of respectively. Notice now in (7), there is a global dual variable in the primal-dual form. Therefore, we also introduce the local copy of the dual variable, i.e., to formulate it into the decentralized optimization problem. Now the final objective function we need to optimize is
where We are now ready to present the value propagation algorithm. In the following, for notational simplicity, we assume the parameter of each agent is a scalar, i.e., . We pack the parameter together and slightly abuse the notation by writing , , Similarly, we also pack the stochastic gradient , .
3.3 Value propagation algorithm
Solving (3.2) even without constraints is not an easy problem when both primal and dual parts are approximated by the deep neural networks. An ideal way is to optimize the inner dual problem and find the solution , such that . Then we can do the (decentralized) stochastic gradient decent to solve the primal problem.
However in practice, one tricky issue is that we can not get the exact solution of the dual problem. Thus, we do the (decentralized) stochastic gradient for steps in the dual problem and get an approximated solution in the Algorithm 3.3. In our analysis, we take the error generated from this inexact solution into the consideration and analyze its effect on the convergence. Particularly, since , the primal gradient is biased and the results in (Dai et al., 2018; Hong et al., 2017) do not fit this problem.
In the dual update we do a consensus update using the stochastic gradient of each agent, where is some auxiliary variable to incorporate the communication, is the degree matrix, is the node-edge incidence matrix, is sign-less graph Laplacian. We defer the detail definition and the derivation of this algorithm to Appendix A.1 and Appendix A.5 due to space limitation. After updating the dual parameters, we optimize the primal parameters , . Similarly, we use a mini-batch data from the replay buffer and then do a consensus update on . The same remarks on also hold for the primal parameter . Notice here we do not need the consensus update on , since each agent’s policy is different than each other. This update rule is adapted from a primal decentralized optimization algorithm (Hong et al., 2017). Notice even in the pure primal case, Hong et al. (2017) only consider the batch gradient case while our algorithm and analysis include the stochastic and biased gradient case. In practicals implementation, we consider the decentralized momentum method and multi-step temporal consistency to accelerate the training; see details in Appendix A.2 and Appendix A.3.
Remarks on Algorithm 3.3. (1) In the single agent case, Dai et al. (2018) assume the dual problem can be exactly solved and thus they analyze a simple pure primal problem. However such assumption is unrealistic especially when the dual variable is represented by the deep neural network. Our multi-agent analysis considers the inexact solution. This is much harder than that in (Dai et al., 2018), since now the primal gradient is biased. (2) The update of each agent just needs the information of the agent itself and its neighbors. See this from the definition of , , in the appendix. (3) The topology of the Graph affects the convergence speed. In particular, the rate depends on and , which are related to spectral gap of the network.
\STATE Input: Environment ENV, learning rate , , , discount factor , number of step to train dual parameter , replay buffer capacity , node-edge incidence matrix , degree matrix , signless graph Laplacian .
\STATEInitialization of , , .
\STATE sample trajectory and add it into the replay buffer.
1. Update the dual parameter
Do following dual update times:
Random sample a mini-batch of transition from the replay buffer.
\STATECalculate the stochastic gradient of w.r.t. .
// Do consensus update on
2. Update primal parameters \STATERandom sample a mini-batch of transition from the replay buffer. \FOR agent to \STATECalculate the stochastic gradient , of , w.r.t. , \ENDFOR
// Do gradient decent on : for each agent .
// Do consensus update on :
4 Theoretical Result
In this section, we give the convergence result on Algorithm 3.3. We first make two mild assumptions on the function approximators of , , .
i) The function approximator is differentiable and has Lipschitz continuous gradient, i.e., This is commonly assumed in the non-convex optimization. ii) The function approximator is lower bounded. This can be easily satisfied when the parameter is bounded, i.e., for some positive constant .
In the following, we give the theoretical analysis for Algorithm 3.3 in the same setting of (Antos et al., 2008; Dai et al., 2018) where samples are prefixed and from one single -mixing off-policy sample path. We denote
Let the function approximators of , and satisfy Assumption 1, snd denote the total training step be . We solve the inner dual problem with a approximated solution , such that , and . Assume the variance of the stochastic gradient , and (estimated by a single sample) are bounded by , the size of the mini-batch is , the step size . Then value propagation in Algorithm 3.3 converges to the stationary solution of with rate
Remarks: (1) The convergence criteria and its dependence on the network structure are involved. We defer the definition of them to the proof section in the appendix (Equation (44)). (2) We require that the approximated dual solution are not far from such that the estimation of the primal gradient of and are not far from the true one (the distance is less than ). Once the inner dual problem is concave, we can get this approximated solution easily using vanilla decentralized stochastic gradient method after at most steps. If the dual problem is non-convex, we still can show the dual problem converges to some stationary solution with rate by our proof. (3) In the theoretical analysis, the stochastic gradient estimated from the mini-batch (rather than the estimation from a single sample ) is common in non-convex analysis, see the work (Ghadimi and Lan, 2016). In practice, a mini-batch of samples is commonly used in training deep neural network.
5 Related work
Among related work on MARL, the setting of (Zhang et al., 2018) is close to ours, where the authors proposed a fully decentralized multi-agent Actor-Critic algorithm to maximize the expected time-average reward . They provide the asymptotic convergence analysis on the on-policy and linear function approximation setting. In our work, we consider the discounted reward setup, i.e., Equation (3). Our algorithm includes both on-policy and off-policy setting thus can exploit data more efficiently. Furthermore, we provide a convergence rate in the non-linear function approximation setting which is much stronger than the result in (Zhang et al., 2018). Littman (1994) proposed the framework of Markov games which can be applied to collaborative and competitive setting (Lauer and Riedmiller, 2000; Hu and Wellman, 2003). These early works considered the tabular case thus can not apply to real problems with large state space. Recent works (Foerster et al., 2016, 2018; Rashid et al., 2018; Raileanu et al., 2018; Jiang et al., 2018; Lowe et al., 2017) have exploited powerful deep learning and obtained some promising empirical results. However most of them lacks theoretical guarantees while our work provides convergence analysis. We emphasize that most of the research on MARL is in the fashion of centralized training and decentralized execution. In the training, they do not have the constraint on the communication, while our work has a network decentralized structure.
6 Experimental result
The goal of our experiment is two-fold: To better understand the effect of each component in the proposed algorithm; and to evaluate efficiency of value propagation in the off-policy setting. To this end, we first do an ablation study on a simple random MDP problem, we then evaluate the performance on the cooperative navigation task (Lowe et al., 2017). The settings of the experiment are similar to those in (Zhang et al., 2018). Some implementation details are deferred to Appendix A.4 due to space constraints.
6.1 Ablation Study
In this experiment, we test effect of several components of our algorithm such as the consensus update, dual formulation in a random MDP problem. Particularly we answer following three questions: (1) Whether an agent can get consensus through message-passing in value propagation even when each agent just knows its local reward. (2) How much performance does the decentralized approach sacrifice comparing with centralized one? (3) What is the effect of the dual part in our formulation ( and corresponds to the pure primal form)?
We compare value propagation with the centralized PCL. The centralized PCL means that there is a central node to collect rewards of all agent, thus it can optimize the objective function (5) using the single agent PCL algorithm (Nachum et al., 2017; Dai et al., 2018). Ideally, value propagation should converges to the same long term reward with the one achieved by the centralized PCL. In the experiment, we consider a multi-agent RL problem with and agents, where each agent has two actions. A discrete MDP is randomly generated with states. The transition probabilities are distributed uniformly with a small additive constant to ensure ergodicity of the MDP, which is . For each agent and each state-action pair , the reward is uniformly sampled from .
In the left panel of Figure 2, we verify that the value function in value propagation reaches the consensus through message-passing in the end of the training. Particularly, we randomly choose three agent , , and draw their value functions over 20 randomly picked states. It is easy to see that value functions , , over these states are almost same. This is accomplished by the consensus update in value propagation. In the middle and right panel of Figure 2, we compare the result of value propagation with centralized PCL and evaluate the effect of the dual part of value propagation. Particularly, we pick in the experiment, where corresponds to the pure primal formulation. When is too large (), the algorithm would have large variance while the algorithm has some bias. Thus value propagation with has better result. We also see that value propagation () and centralized PCL converge to almost the same value, although there is a gap between centralized and decentralized algorithm. The centralized PCL converges faster than value propagation, since it does not need time to diffuse the reward information over the network.
6.2 Cooperative Navigation task
The aim of this section is to demonstrate that the value propagation outperforms decentralized multi-agent Actor-Critic (MA-AC)(Zhang et al., 2018), independent Q learning (Tan, 1993), the Multi-agent PCL without communication. Here PCL without communication means each agent maintains its own estimation of policy and value function but there is no communication Graph. Notice that this is different from the centralized PCL in Section 6.1, where centralized PCL has a central node to collect all reward information and thus do not need further communication. Note that the original MA-AC is designed for the averaged reward setting thus we adapt it into the discounted case to fit our setting. We test the value propagation in the environment of the Cooperative Navigation task (Lowe et al., 2017), where agents need to reach a set of landmarks through physical movement. We modify this environment to fit our setting. A reward is given when the agent reaches its own landmarks. A penalty is received if agents collide with other agents. Since the position of landmarks are different, the reward function of each agent is different. Here we test the case the state is globally observed and partially observed. In particular, we assume the environment is in a rectangular region with size . There are or agents. Each agent has a single target landmark, i.e., , which is randomly located in the region. Each agent has five actions which corresponds to going up, down, left, right with units 0.1 or staying at the position. The agent has high probability (0.95) to move in the direction following its action and go in other direction randomly otherwise. The maximum length of each epoch is set to be 500 steps. When the agent is close enough to the landmark, e.g., the distance is less than 0.1, we think it reaches the target and gets reward . When two agents are close to each other (with distance less than 0.1), we treat this case as a collision and a penalty is received for each of the agents. The state includes the position of the agents. The communication graph is generated as that in Section 6.1 with connectivity ratio . In the partially observed case, the actor of each agent can only observe its own and neighbors’ states. We report the results in Figure 3.
In the left panel of Figure 3, we see the value function reaches consensus in value propagation. In the middle and right panel of Figure 3, we compare value propagation with PCL without communication, independent Q learning and MA-AC. In PCL without communication, each agent maintains its own policy, value function and dual function, which is trained by the algorithm SBEED (Dai et al., 2018) with . Since there is no communication between agents, intuitively agents may have more collisions in the learning process than those in value propagation. Similar augment holds for the independent Q learning. Indeed, In the middle and right panel, we see value propagation learns the policy much faster than PCL without communication. We also observe that value propagation outperforms MA-AC. One possible reason is that value propagation is an off-policy method thus we can apply experience replay which exploits data more efficiently than the on-policy method MA-AC. We also test the performance of value propagation (result labeled as partial value propagation in Figure 3) when the state information of actor is partially observed. Since the agent has limited information, its performance is worse than the fully observed case. But it is better than the PCL without communication (fully observed state).
- Abadi et al. (2016) Martin Abadi, Andy Chu, Ian Goodfellow, H Brendan McMahan, Ilya Mironov, Kunal Talwar, and Li Zhang. Deep learning with differential privacy. In Proceedings of the 2016 ACM SIGSAC Conference on Computer and Communications Security, pages 308–318. ACM, 2016.
- Adler and Blue (2002) Jeffrey L Adler and Victor J Blue. A cooperative multi-agent transportation management and route guidance system. Transportation Research Part C: Emerging Technologies, 10(5-6):433–454, 2002.
- Antos et al. (2008) András Antos, Csaba Szepesvári, and Rémi Munos. Learning near-optimal policies with bellman-residual minimization based fitted policy iteration and a single sample path. Machine Learning, 71(1):89–129, 2008.
- Boyd et al. (2006) Stephen Boyd, Arpita Ghosh, Balaji Prabhakar, and Devavrat Shah. Randomized gossip algorithms. IEEE transactions on information theory, 52(6):2508–2530, 2006.
- Bu et al. (2008) Lucian Bu, Robert Babu, Bart De Schutter, et al. A comprehensive survey of multiagent reinforcement learning. IEEE Transactions on Systems, Man, and Cybernetics, Part C (Applications and Reviews), 38(2):156–172, 2008.
- Callaway and Hiskens (2011) Duncan S Callaway and Ian A Hiskens. Achieving controllability of electric loads. Proceedings of the IEEE, 99(1):184–199, 2011.
- Cattivelli et al. (2008) Federico S Cattivelli, Cassio G Lopes, and Ali H Sayed. Diffusion recursive least-squares for distributed estimation over adaptive networks. IEEE Transactions on Signal Processing, 56(5):1865–1877, 2008.
- Cortes et al. (2004) Jorge Cortes, Sonia Martinez, Timur Karatas, and Francesco Bullo. Coverage control for mobile sensing networks. IEEE Transactions on robotics and Automation, 20(2):243–255, 2004.
- Dai et al. (2018) Bo Dai, Albert Shaw, Lihong Li, Lin Xiao, Niao He, Zhen Liu, Jianshu Chen, and Le Song. Sbeed: Convergent reinforcement learning with nonlinear function approximation. In International Conference on Machine Learning, pages 1133–1142, 2018.
- Dann et al. (2014) Christoph Dann, Gerhard Neumann, and Jan Peters. Policy evaluation with temporal differences: A survey and comparison. The Journal of Machine Learning Research, 15(1):809–883, 2014.
- DeGroot (1974) Morris H DeGroot. Reaching a consensus. Journal of the American Statistical Association, 69(345):118–121, 1974.
- Duchi et al. (2011) John Duchi, Elad Hazan, and Yoram Singer. Adaptive subgradient methods for online learning and stochastic optimization. Journal of Machine Learning Research, 12(Jul):2121–2159, 2011.
- El-Tantawy et al. (2013) Samah El-Tantawy, Baher Abdulhai, and Hossam Abdelgawad. Multiagent reinforcement learning for integrated network of adaptive traffic signal controllers (marlin-atsc): methodology and large-scale application on downtown toronto. IEEE Transactions on Intelligent Transportation Systems, 14(3):1140–1150, 2013.
- Fax and Murray (2002) J Alexander Fax and Richard M Murray. Information flow and cooperative control of vehicle formations. In IFAC World Congress, volume 22, 2002.
- Foerster et al. (2016) Jakob Foerster, Ioannis Alexandros Assael, Nando de Freitas, and Shimon Whiteson. Learning to communicate with deep multi-agent reinforcement learning. In Advances in Neural Information Processing Systems, pages 2137–2145, 2016.
- Foerster et al. (2018) Jakob N Foerster, Gregory Farquhar, Triantafyllos Afouras, Nantas Nardelli, and Shimon Whiteson. Counterfactual multi-agent policy gradients. In Thirty-Second AAAI Conference on Artificial Intelligence, 2018.
- Ghadimi and Lan (2016) Saeed Ghadimi and Guanghui Lan. Accelerated gradient methods for nonconvex nonlinear and stochastic programming. Mathematical Programming, 156(1-2):59–99, 2016.
- Hong (2016) Mingyi Hong. Decomposing linearly constrained nonconvex problems by a proximal primal dual approach: Algorithms, convergence, and applications. arXiv preprint arXiv:1604.00543, 2016.
- Hong et al. (2017) Mingyi Hong, Davood Hajinezhad, and Ming-Min Zhao. Prox-pda: The proximal primal-dual algorithm for fast distributed nonconvex optimization and learning over networks. In International Conference on Machine Learning, pages 1529–1538, 2017.
- Hu and Wellman (2003) Junling Hu and Michael P Wellman. Nash q-learning for general-sum stochastic games. Journal of machine learning research, 4(Nov):1039–1069, 2003.
- Jiang et al. (2018) Jiechuan Jiang, Chen Dun, and Zongqing Lu. Graph convolutional reinforcement learning for multi-agent cooperation. arXiv preprint arXiv:1810.09202, 2018.
- Kingma and Ba (2014) Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980, 2014.
- Kurakin et al. (2016) Alexey Kurakin, Ian Goodfellow, and Samy Bengio. Adversarial machine learning at scale. arXiv preprint arXiv:1611.01236, 2016.
- Lauer and Riedmiller (2000) Martin Lauer and Martin Riedmiller. An algorithm for distributed reinforcement learning in cooperative multi-agent systems. In In Proceedings of the Seventeenth International Conference on Machine Learning. Citeseer, 2000.
- Littman (1994) Michael L Littman. Markov games as a framework for multi-agent reinforcement learning. In Machine Learning Proceedings 1994, pages 157–163. Elsevier, 1994.
- Lowe et al. (2017) Ryan Lowe, Yi Wu, Aviv Tamar, Jean Harb, OpenAI Pieter Abbeel, and Igor Mordatch. Multi-agent actor-critic for mixed cooperative-competitive environments. In Advances in Neural Information Processing Systems, pages 6379–6390, 2017.
- McMahan et al. (2016) H Brendan McMahan, Eider Moore, Daniel Ramage, Seth Hampson, et al. Communication-efficient learning of deep networks from decentralized data. arXiv preprint arXiv:1602.05629, 2016.
- Mo and Murray (2017) Yilin Mo and Richard M Murray. Privacy preserving average consensus. IEEE Transactions on Automatic Control, 62(2):753–765, 2017.
- Nachum et al. (2017) Ofir Nachum, Mohammad Norouzi, Kelvin Xu, and Dale Schuurmans. Bridging the gap between value and policy based reinforcement learning. In Advances in Neural Information Processing Systems, pages 2775–2785, 2017.
- Rabbat and Nowak (2004) Michael Rabbat and Robert Nowak. Distributed optimization in sensor networks. In Proceedings of the 3rd international symposium on Information processing in sensor networks, pages 20–27. ACM, 2004.
- Raileanu et al. (2018) Roberta Raileanu, Emily Denton, Arthur Szlam, and Rob Fergus. Modeling others using oneself in multi-agent reinforcement learning. arXiv preprint arXiv:1802.09640, 2018.
- Rashid et al. (2018) Tabish Rashid, Mikayel Samvelyan, Christian Schroeder de Witt, Gregory Farquhar, Jakob Foerster, and Shimon Whiteson. Qmix: Monotonic value function factorisation for deep multi-agent reinforcement learning. arXiv preprint arXiv:1803.11485, 2018.
- Shapiro et al. (2009) Alexander Shapiro, Darinka Dentcheva, and Andrzej Ruszczyński. Lectures on stochastic programming: modeling and theory. SIAM, 2009.
- Shi et al. (2015) Wei Shi, Qing Ling, Gang Wu, and Wotao Yin. Extra: An exact first-order algorithm for decentralized consensus optimization. SIAM Journal on Optimization, 25(2):944–966, 2015.
- Tan (1993) Ming Tan. Multi-agent reinforcement learning: Independent vs. cooperative agents. In Proceedings of the tenth international conference on machine learning, pages 330–337, 1993.
- (36) T Tieleman and G Hinton. Divide the gradient by a running average of its recent magnitude. coursera: Neural networks for machine learning. Technical report, Technical Report. Available online: https://zh. coursera. org/learn ….
- Xiao et al. (2005) Lin Xiao, Stephen Boyd, and Sanjay Lall. A scheme for robust distributed sensor fusion based on average consensus. In Information Processing in Sensor Networks, 2005. IPSN 2005. Fourth International Symposium on, pages 63–70. IEEE, 2005.
- Zhang et al. (2018) Kaiqing Zhang, Zhuoran Yang, Han Liu, Tong Zhang, and Tamer Başar. Fully decentralized multi-agent reinforcement learning with networked agents. International Conference on Machine Learning, 2018.
Appendix A Details on the Value Propagation
A.1 Topology of the Graph
Here, we explain the matrix in the Algorithm 3.3 which are closely related to the topology of the Graph, which is left from the main paper due to the limit of the space.
is the degree matrix, with denoting the degree of node .
is the node-edge incidence matrix: if and it connects vertex i and j with , then if , if and otherwise.
The signless incidence matrix , where the absolute value is taken for each component of .
The signless graph Laplacian . By definition if . Notice the non-zeros element in , , the update just depends on each agent itself and its neighbor.
A.2 Practical Acceleration
The algorithm 3.3 trains the agent with vanilla gradient decent method with a extra consensus update. In practice, the adaptive momentum gradient methods including Adagrad Duchi et al. , Rmsprop Tieleman and Hinton and Adam Kingma and Ba  have much better performance in training the deep neural network. We adapt Adam in our setting, and propose algorithm A.2 which has better performance than algorithm 3.3 in practice.
\STATE Input: Environment ENV, learning rate , , discount factor , a mixing matrix , number of step to train dual parameter , replay buffer capacity .
\STATEInitialization of , moment vectors , .
\STATE sample trajectory and add it into the replay buffer.
// Update the dual parameter
Do following update times: \STATERandom sample a mini-batch of transition from the replay buffer.
\STATECalculate the stochastic gradient of w.r.t. .
// update momentum parameters:
// Do consensus update for each agent
// End the update of dual problem
// Update primal parameters . \STATERandom sample a mini-batch of transition from the replay buffer. \FOR agent to \STATECalculate the stochastic gradient , of , w.r.t. ,
//update the momentum parameter:
// Using Adam to update for each agent .
// Do consensus update on for each agent :
Mixing Matrix: In Algorithm A.2, there is a mixing matrix in the consensus update. As its name suggests, it mixes information of the agent and its neighbors. This nonnegative matrix need to satisfy the following condition.
needs to be doubly stochastic, i.e., and .
respects the communication graph , i.e., if .
The spectral norm of is strictly smaller than one.
Here is one particular choice of the mixing matrix used in our work which satisfies above requirement called Metropolis weights Xiao et al. .
where is the set of neighbors of the agent and is the degree of agent . Such mixing matrix is widely used in decentralized and distributed optimization Boyd et al. , Cattivelli et al. . The update rule of the momentum term in Algorithm A.2 is adapted from Adam. The consensus (communication) steps are and .
A.3 Multi-step Extension on value propagation
The temporal consistency can be extended to the multi-step case Nachum et al. , where the following equation holds
Thus in the objective function (3.2), we can replace by and change the estimation of stochastic gradient correspondingly in Algorithm 3.3 and Algorithm A.2 to get the multi-step version of vaue propagation . In practice, the performance of setting is better than which is also observed in single agent case Nachum et al. , Dai et al. . We can tune for each application to get the best performance.
A.4 Implementation details of the experiments
The value function and dual variable are approximated by two hidden-layer neural network with Relu as the activation function where each hidden-layer has hidden units. The policy of each agent is approximated by a one hidden-layer neural network with Relu as the activation function where the number of the hidden units is . The output is the softmax function to approximate . The mixing matrix in Algorithm A.2 is selected as the Metropolis Weights in (10). The graph is generated by randomly placing communication links among agents such that the connectivity ratio is . We set , , learning rate =5e-4. The choice of , are the default value in Adam.
Cooperative Navigation task
The value function is approximated by a two-hidden-layer neural network with Relu as the activation function where inputs are the state information. Each hidden-layer has hidden units. The dual function is also approximated by a two-hidden-layer neural network, where the only difference is that inputs are state-action pairs (s,a). The policy is approximated by a one-hidden-layer neural network with Relu as the activation function. The number of the hidden units is . The output is the softmax function to approximate . In all experiments, we use the multi-step version of value propagation and choose . We choose , . The learning rate of Adam is chosen as 5e-4 and are default value in Adam optimizer. The setting of PCL without communication is exactly same with value propagation except the absence of communication network.
A.5 Consensus update in Algorithm 3.3
We now give details to derive the Consensus Update in Algorithm 3.3 with to ease the exposition. When , we just need to change variable and some notations, the result are almost same. Here we use the primal update as an example, the derivation of the dual update is the same.
In the main paper section 3, we have shown that when , in the primal update, we basically solve following problem.
here for simplicity we assume in the dual optimization, we have already find the optimal solution . It can be any approximated solution of which does not affect the derivation of the update rule in primal optimization. In the later proof, we will show how this approximated solution affects the convergence rate.
When we optimize w.r.t. , we basically we solve a non-convex problem with the following form
Recall the definition of the node-edge incidence matrix : if and it connects vertex and with , then if , if and otherwise. Thus by define we have a equivalent form of (12)
Notice the update of is a special case of above formulation, since we do not have the constraint . Thus in the following, it suffice to analyze above formulation (13). We adapt the Prox-PDA in Hong et al.  to solve above problem. To keep the notation consistent with Hong et al. , we consider a more general problem
In the following we denote where the superscript means transpose. We denote as an estimator of and .
The update rule of Prox-PDA is
where is an estimator of . The signed graph Laplacian matrix is . Now we choose as the signless incidence matrix. Using this choice of , we have which is the signless graph Laplacian whose th diagonal entry is the degree of node and its th entry is 1 if , and 0 otherwise.
where is the degree matrix, with denoting the degree of node .
After simple algebra, we obtain
which is the primal update rule of the consensus step in the algorithm 3.3 (notice here the stepsize is )
Appendix B Convergence Proof of Value Propagation
B.1 Convergence on the primal update
In this section, we first give the convergence analysis of the value propagation (algorithm 3.3) on the primal update. To include the effected of the inexact solution of dual optimization problem, we denote , where is some error terms.
is a zero mean random variable coming from the randomness of the stochastic gradient .
Before we begin the proof, we made some mild assumption on the function .
1. The function f(x) is differentiable and has Lipschitz continuous gradient, i.e.,
2. Further assume that . This assumption is always satisfied by our choice on A and B. We have
3. There exists a constant such that This assumption is satisfied if we require the parameter space is bounded.
Suppose the assumption 2 is satisfied, we have following inequality holds
Using the optimality condition of (14), we obtain
applying equation (15) we have
Note that from the fact that , we have the variable lies in the column space of .
Let denote the smallest non-zero eigenvalue of , we have
Thus we have
where the second inequality holds from the fact that .
Define . Suppose assumptions are satisfied, then the following is true for the algorithm
By the Assumptions , the objective function in (14) is strongly convex with parameter .
Using the optimality condition of and strong convexity, we have for any ,
Now we start to provide a upper bound of .
replace by , we have the result. The inequality (c) holds using the Lemma 1.
Suppose Assumption 2 is satisfied, then the following condition holds.
Using the optimality condition of and in the update rule in (14), we obtain
Replacing by and by , and using the update rule (15)
Now choose in the first inequality and in the second one, adding two inequalities together, we obtain
Rearranging above terms, we have
We first re-express the lhs of above inequality.
Next, we bound the rhs of (31).
where the inequality (a) uses Cauchy-Schwartz inequality, (b) holds from the smoothness assumption on .
Combine all pieces together, we obtain
Same with Hong et al. , we define the potential function
If Assumption 2 holds, we have following
where the second inequality holds from the Cauchy-Schwartz inequality.
We require that
which is satisfied when
We further require
which will be used later in the telescoping.
Thus we require
and choose ∎
Now we do summation over both side of (36) and have
rearrange terms of above inequality.
Next we show is lower bounded
The following lemma is from Lemma 3.5 in Hong , we present here for completeness.
Sum over both side, we obtain
By assumption 2, above sum is lower bounded, which implies that the sum of the potential function is also lower bounded (Recall ). Thus we have
In the next step, we are ready to provide the convergence rate. Following Hong , we define the convergence criteria
It is easy to see, when , and , which are KKT condition of the problem.
Using the proof in Lemma 1, we know there exist two positive constants c1 c2 c3 c4
Using Lemma 4, we know there must exist a constant such that
Divide both side by and take expectation
Now we bound the R.H.S. of above equation.
Recall we choose the mini-batch size , and
Similarly we can bound . Combine all pieces together, we obtain
where are some universal positive constants.
Notice , we have where is a universal positive constant.
B.2 Convergence on the dual update
If the dual objective function is non-convex, we just follow the exact analysis in our proof on the primal problem. Notice the analysis on the dual update is easier than primal one, since we do not have the error term . Therefore, we have the algorithm converges to stationary solution with rate in criteria .