Distributional Reward Decomposition for Reinforcement Learning

  • 2019-11-06 02:13:50
  • Zichuan Lin, Li Zhao, Derek Yang, Tao Qin, Guangwen Yang, Tie-Yan Liu
  • 0

Abstract

Many reinforcement learning (RL) tasks have specific properties that can beleveraged to modify existing RL algorithms to adapt to those tasks and furtherimprove performance, and a general class of such properties is the multiplereward channel. In those environments the full reward can be decomposed intosub-rewards obtained from different channels. Existing work on rewarddecomposition either requires prior knowledge of the environment to decomposethe full reward, or decomposes reward without prior knowledge but with degradedperformance. In this paper, we propose Distributional Reward Decomposition forReinforcement Learning (DRDRL), a novel reward decomposition algorithm whichcaptures the multiple reward channel structure under distributional setting.Empirically, our method captures the multi-channel structure and discoversmeaningful reward decomposition, without any requirements on prior knowledge.Consequently, our agent achieves better performance than existing methods onenvironments with multiple reward channels.

 

Quick Read (beta)

Distributional Reward Decomposition for Reinforcement Learning

Zichuan Lin
Tsinghua University
[email protected]
&Li Zhao
Microsoft Research
[email protected]
&Derek Yang
UC San Diego
[email protected]
\ANDTao Qin
Microsoft Research
[email protected]
&Guangwen Yang
Tsinghua University
[email protected]
&Tie-Yan Liu
Microsoft Research
[email protected]
Contributed during internship at Microsoft Research.
Abstract

Many reinforcement learning (RL) tasks have specific properties that can be leveraged to modify existing RL algorithms to adapt to those tasks and further improve performance, and a general class of such properties is the multiple reward channel. In those environments the full reward can be decomposed into sub-rewards obtained from different channels. Existing work on reward decomposition either requires prior knowledge of the environment to decompose the full reward, or decomposes reward without prior knowledge but with degraded performance. In this paper, we propose Distributional Reward Decomposition for Reinforcement Learning (DRDRL), a novel reward decomposition algorithm which captures the multiple reward channel structure under distributional setting. Empirically, our method captures the multi-channel structure and discovers meaningful reward decomposition, without any requirements on prior knowledge. Consequently, our agent achieves better performance than existing methods on environments with multiple reward channels.

 

Distributional Reward Decomposition for Reinforcement Learning


  Zichuan Linthanks: Contributed during internship at Microsoft Research. Tsinghua University [email protected] Li Zhao Microsoft Research [email protected] Derek Yang UC San Diego [email protected] Tao Qin Microsoft Research [email protected] Guangwen Yang Tsinghua University [email protected] Tie-Yan Liu Microsoft Research [email protected]

\@float

noticebox[b]33rd Conference on Neural Information Processing Systems (NeurIPS 2019), Vancouver, Canada.\[email protected]

1 Introduction

Reinforcement learning has achieved great success in decision making problems since Deep Q-learning was proposed by mnih2015human. While general RL algorithms have been deeply studied, here we focus on those RL tasks with specific properties that can be utilized to modify general RL algorithms to achieve better performance. Specifically, we focus on RL environments with multiple reward channels, but only the full reward is available.

Reward decomposition has been proposed to investigate such properties. For example, in Atari game Seaquest, rewards of environment can be decomposed into sub-rewards of shooting sharks and those of rescuing divers. Reward decomposition views the total reward as the sum of sub-rewards that are usually disentangled and can be obtained independently (sprague2003multiple; russell2003q; van2017hybrid; grimm2019learning), and aims at decomposing the total reward into sub-rewards. The sub-rewards may further be leveraged to learn better policies.

van2017hybrid propose to split a state into different sub-states, each with a sub-reward obtained by training a general value function, and learn multiple value functions with sub-rewards. The architecture is rather limited due to requiring prior knowledge of how to split into sub-states. grimm2019learning propose a more general method for reward decomposition via maximizing disentanglement between sub-rewards. In their work, an explicit reward decomposition is learned via maximizing the disentanglement of two sub-rewards estimated with action-value functions. However, their work requires that the environment can be reset to arbitrary state and can not apply to general RL settings where states can hardly be revisited. Furthermore, despite the meaningful reward decomposition they achieved, they fail to utilize the reward decomposition into learning better policies.

In this paper, we propose Distributional Reward Decomposition for Reinforcement Learning (DRDRL), an RL algorithm that captures the latent multiple-channel structure for reward, under the setting of distributional RL. Distributional RL differs from value-based RL in that it estimates the distribution rather than the expectation of returns, and therefore captures richer information than value-based RL. We propose an RL algorithm that estimates distributions of the sub-returns, and combine the sub-returns to get the distribution of the total returns. In order to avoid naive decomposition such as 0-1 or half-half, we further propose a disentanglement regularization term to encourage the sub-returns to be diverged. To better separate reward channels, we also design our network to learn different state representations for different channels.

We test our algorithm on chosen Atari Games with multiple reward channels. Empirically, our method has following achievements:

  • Discovers meaningful reward decomposition.

  • Requires no external information.

  • Achieves better performance than existing RL methods.

2 Background

We consider a general reinforcement learning setting, in which the interaction of the agent and the environment can be viewed as a Markov Decision Process (MDP). Denote the state space by 𝒳, action space by A, the state transition function by P, the action-state dependent reward function by R and γ the discount factor, we write this MDP as (𝒳,A,R,P,γ).

Given a fixed policy π, reinforcement learning estimates the action-value function of π, defined by Qπ(x,a)=t=0γtrt(xt,at) where (xt,at) is the state-action pair at time t, x0=x,a0=a and rt is the corresponding reward. The Bellman equation characterizes the action-value function by temporal equivalence, given by

Qπ(x,a)=R(x,a)+γ𝔼x,a[Qπ(x,a)]

where xP(|x,a),aπ(|x). To maximize the total return given by 𝔼x0,a0[Qπ(x0,a0)], one common approach is to find the fixed point for the Bellman optimality operator

Q*(x,a)=𝒯Q*(x,a)=R(x,a)+γ𝔼x[maxaQ*(x,a)]

with the temporal difference (TD) error, given by

δt2=[rt+γmaxa𝒜Q(xt+1,a)-Q(xt,at)]2

over the samples (xt,at,st,xt+1) along the trajectory. mnih2015human propose Deep Q-Networks (DQN) that learns the action-value function with a neural network and achieves human-level performance on the Atari-57 benchmark.

2.1 Reward Decomposition

Studies for reward decomposition also leads to state decomposition (laversanne2018curiosity; thomas2017independently), where state decomposition is leveraged to learn different policies. Extending their work, grimm2019learning explore the decomposition of the reward function directly, which is considered to be most related to our work. Denote the i-th (i=1,2,…,N) sub-reward function at state-action pair (x,a) as ri(x,a), the complete reward function is given by

r=i=1Nri

For each sub-reward function, consider the sub-value function Uiπ and corresponding policy πi:

Uiπ(x0,a0)=𝔼xt,at[t=0γtri(xt,at)]
πi=argmaxπUiπ

where xtP(|π,x0,a0),atπ(|xt).

In their work, reward decomposition is considered meaningful if each reward is obtained independently (i.e. πi should not obtain rj) and each reward is obtainable.

Two evaluate the two desiderata, the work proposes the following values:

Jindependent(r1,,rn)=𝔼sμ[ijαi,j(s)Uiπj*(s)], (1)
Jnontrivial(r1,,rn)=𝔼sμ[i=1nαi,i(s)Uiπi*(s)], (2)

where αi,j is for weight control and for simplicity set to 1 in their work. During training, the network would maximize Jnontrivial-Jindependent to achieve the desired reward decomposition.

2.2 Distributional Reinforcement Learning

In most reinforcement learning settings, the environment is not deterministic. Moreover, in general people train RL models with an ϵ-greedy policy to allow exploration, making the agent also stochastic. To better analyze the randomness under this setting, bellemare2017distributional propose C51 algorithm and conduct theoretical analysis of distributional RL.

In distributional RL, reward Rt is viewed as a random variable, and the total return is defined by Z=t=0γtRt. The expectation of Z is the traditional action-value Q and the distributional Bellman optimality operator is given by

𝒯Z(x,a):=DR(x,a)+γZ(x,argmaxa𝒜𝔼Z(x,a))

where if random variable A and B satisfies A=DB then A and B follow the same distribution.

Random variable is characterized by a categorical distribution over a fixed set of values in C51, and it outperforms all previous variants of DQN on Atari domain.

3 Distributional Reward Decomposition for Reinforcement Learning

3.1 Distributional Reward Decomposition

Figure 1: (a) Distributional reward decomposition network architecture. (b) Examples of multiple reward channels in Atari games: the top row shows examples of Seaquest in which the submarine receives rewards from both shooting sharks and rescuing divers; the bottom row shows examples of Hero where the hero receives rewards from both shooting bats and rescuing people.

In many reinforcement learning environments, there are multiple sources for an agent to receive reward as shown in Figure 1. Our method is mainly designed for environments with such property.

Under distributional setting, we will assume reward and sub-rewards are random variables and denote them by R and Ri respectively. In our architecture, the categorical distribution of each sub-return Zi=t=0γtRi,t is the output of a network, denoted by i(x,a). Note that in most cases, sub-returns are not independent, i.e. P(Zi=v)!=P(Zi=v|Zj). So theoretically we need ij(x,a) for each i and j to obtain the distribution of the full return. We call this architecture as non-factorial model or full-distribution model. The non-factorial model architecture is shown in appendix. However, experiment shows that using an approximation form of P(Zi=v)P(Zi=v|Zj) so that only i(x,a) is required performs much better than brutally computing ij(x,a) for all i,j, we believe that is due to the increased sample number. In this paper, we approximate the conditional probability P(Zi=v|Zj) with P(Zi=v).

Consider categorical distribution function i and j with same number of atoms K, the k-th atom is denoted by ak with value ak=a0+kl,1kK where l is a constant. Let random variable Zii and Zjj, from basic probability theory we know that the distribution function of Z=Zi+Zj is the convolution of Fi and Fj

(v)=P(Zi+Zj=v)=k=1KP(Zi=ak)P(Zj=v-ak|Zi=ak)k=1KP(Zi=ak)P(Zj=v-ak)=i(v)*j(v). (3)

When we use N sub-returns, the distribution function of the total return is then given by =1*2**N where * denotes linear 1D-convolution.

While reward decomposition is not explicitly done in our algorithm, we can derive the decomposed reward with using trained agents. Recall that total return Z=i=1NZi follows bellman equation, so naturally we have

𝒯Z=D𝒯(i=1NZi)=DR+γZ=(i=1NRi)+γ(i=1NZi) (4)

where Zi represents sub-return on the next state-action pair. Note that we only have access to a sample of the full reward R, the sub-rewards Ri are arbitrary and for better visualization a direct way of deriving them is given by

Ri=Zi-γZi (5)

In the next section we will present an example of those sub-rewards by taking their expectation 𝔼(Ri). Note that our reward decomposition is latent and we do not need Ri for our algorithm, Eq. 5 only provides an approach to visualize our reward decomposition.

3.2 Disentangled Sub-returns

To obtain meaningful reward decomposition, we want the sub-rewards to be disentangled. Inspired by grimm2019learning, we compute the disentanglement of distributions of two sub-returns Fi and Fj on state x with the following value:

Jdisentangij=DKL(x,argmaxa𝔼(Zi)||x,argmaxa𝔼(Zj)), (6)

where DKL denotes the cross-entropy term of KL divergence.

Intuitively, Jdisentangij estimates the disentanglement of sub-returns Zi and Zj by first obtaining actions that maximize 𝔼(Zi) and 𝔼(Zj) respectively, and then computes the KL divergence between the two estimated total returns of the actions. If Zi and Zj are independent, the action maximizing two sub-returns would be different and such difference would reflect in the estimation for total return. Through maximizing this value, we can expect a meaningful reward decomposition that learns independent rewards.

3.3 Projected Bellman Update with Regularization

Following the work of C51 (bellemare2017distributional), we use projected Bellman update for our algorithm. When applied with the Bellman optimality operator, the atoms of 𝒯Z is shifted by rt and shrank by γ. However to compute the loss, usually KL divergence between Z and 𝒯Z, it is required that the two categorical distributions are defined on the same set of atoms, so the target distribution 𝒯Z would need to be projected to the original set of atoms before Bellman update. Consider a sample transition (x,a,r,x), the projection operator Φ proposed in C51 is given by

(Φ𝒯Z(x,a))i=j=0M-1[1-|[r+γaj]VminVmax-ai|l]01x,a(aj), (7)

where M is the number of atoms in C51 and []ab bounds its argument in [a,b]. The sample loss for (x,a,r,x) would be given by the cross-entropy term of KL divergence of Z and Φ𝒯Z

x,a,r,x=DKL(Φ𝒯Z(x,a)||Z(x,a)). (8)

Let θ be a neural network parameterized by θ, we combine distributional TD error and disentanglement to jointly update θ. For each sample transition (x,a,r,x), θ is updated by minimizing the following objective function:

x,a,r,x-λij!=iJdisentangij, (9)

where λ denotes learning rate.

3.4 Multi-channel State Representation

One complication of our approach outlined above is that very often the distribution i cannot distinguish itself from other distributions (e.g., j,ji) during learning since they all depend on the same state feature input. This brings difficulties in maximizing disentanglement by jointly training as different distribution functions are exchangeable. A naive idea is to split the state feature ψ(x) into N pieces (e.g., ψ(x)1, ψ(x)2, …, ψ(x)N) so that each distribution depends on different sub-state-features. However, we empirically found that this method is not enough to help learn good disentangled sub-returns.

To address this problem, we utilize an idea similar to universal value function approximation (UVFA) (schaul2015universal). The key idea is to take one-hot embedding as an additional input to condition the categorical distribution function, and apply the element-wise multiplication ψϕ, to force interaction between state features and the one-hot embedding feature:

i(x,a)=θi(ψ(x)ϕ(ei))a, (10)

where ei denotes the one-hot embedding where the i-th element is one and ϕ denotes one-layer non-linear neural network that is updated by backpropagation during training.

In this way, the agent explicitly learns different distribution functions for different channels. The complete network architecture is shown in Figure 1.

4 Experiment Results

We tested our algorithm on the games from Arcade Learning Environment (ALE; bellemare2013arcade). We conduct experiments on six Atari games, some with complicated rules and some with simple rules. For our study, we implemented our algorithm based on Rainbow (hessel2018rainbow) which is an advanced variant of C51 (bellemare2017distributional) and achieved state-of-the-art results in Atari games domain. We replace the update rule of Rainbow by Eq. 9 and network architecture of Rainbow by our convoluted architecture as shown in Figure 1. In Rainbow, the Q-value is bounded by [Vmin,Vmax] where Vmax=-Vmin=10. In our method, we bound the categorical distribution of each sub-return Zi(i=1,2,,N) by a range of [VminN,VmaxN]. Rainbow uses a categorical distribution with M=51 atoms. For fair comparison, we assign K=MN atoms for the distribution of each sub-return, which results in the same network capacity as the original network architecture.

Our code is built upon dopamine framework (castro18dopamine). We use the default well-tuned hyper-parameter setting in dopamine. For our updating rule in Eq. 9, we set λ=0.0001. We run our agents for 100 epochs, each with 0.25 million of training steps and 0.125 million of evaluation steps. For evaluation, we follow common practice in van2016deep, starting the game with up to 30 no-op actions to provide random starting positions for the agent. All experiments are performed on NVIDIA Tesla V100 16GB graphics cards.

4.1 Comparison with Rainbow

To verify that our architecture achieves reward decomposition without degraded performance, we compare our methods with Rainbow. However we are not able to compare our method with van2017hybrid and grimm2019learning since they require either pre-defined state pre-processing or specific-state resettable environments. We test our reward decomposition (RD) with 2 and 3 channels (e.g., RD(2), RD(3)). The results are shown in Figure 2. We found that our methods perform significantly better than Rainbow on the environments that we tested. This implies that our distributional reward decomposition method can help accelerate the learning process. We also discover that on some environments, RD(3) performs better than RD(2) while in the rest the two have similar performance. We conjecture that this is due to the intrinsic settings of the environments. For example, in Seaquest and UpNDown, the rules are relatively complicated, so RD(3) characterizes such complex reward better. However in simple environments like Gopher and Asterix, RD(2) and RD(3) obtain similar performance, and sometimes RD(2) even outperforms RD(3).

Figure 2: Performance comparison with Rainbow. RD(N) represents using N-channel reward decomposition. Each training curve is averaged by three random seeds.

4.2 Reward Decomposition Analysis

Here we use Seaquest to illustrate our reward decomposition. Figure 3 shows the sub-rewards obtained by taking the expectation of the LHS of Eq.5 and the original reward along an actual trajectory. We observe that while r1=𝔼(R1) and r2=𝔼(R2) basically add up to the original reward r, r1 dominates as the submarine is close to the surface, i.e. when it rescues the divers and refills oxygen. When the submarine scores by shooting sharks, r2 becomes the main source of reward. We also monitor the distribution of different sub-returns when the agent is playing game. In Figure 4 (a), the submarine floats to the surface to rescue the divers and refill oxygen and Z1 has higher values. While in Figure 4 (b), as the submarine dives into the sea and shoots sharks, expected values of Z2 (orange) are higher than Z1 (blue). This result implies that the reward decomposition indeed captures different sources of returns, in this case shooting sharks and rescuing divers/refilling oxygen. We also provide statistics on actions for quantitative analysis to support the argument. In Figure 6(a), we count the occurrence of each action obtained with argmaxa𝔼(Z1) and argmaxa𝔼(Z2) in a single trajectory, using the same policy as in Figure 4. We see that while Z1 prefers going up, Z2 prefers going down with fire.

Figure 3: Reward decomposition along the trajectory. While sub-rewards r1 and r2 usually adds up to the original reward r, we see that the proportion of sub-rewards greatly depends on how the original reward is obtained.
Figure 4: An illustration of how the sub-returns discriminates at different stage of the game. In figure (a), the submarine is refilling oxygen while in figure (b) the submarine is shooting sharks.

4.3 Visualization by saliency maps

Figure 5: Sub-distribution saliency maps on the Atari game Seaquest, for a trained DRDRL of two channels (N=2). One channel learns to pay attention to the oxygen, while another channel learns to pay attention to the sharks.

To better understand the roles of different sub-rewards, we train a DRDRL agent with two channels (N=2) and compute saliency maps (simonyan2013deep). Specifically, to visualize the salient part of the images as seen by different sub-policies, we compute the absolute value of the Jacobian |xQi(x,argmaxaQ(x,a))| for each channel. Figure 5 shows that visualization results. We find that channel 1 (red region) focuses on refilling oxygen while channel 2 (green region) pays more attention to shooting sharks as well as the positions where sharks are more likely to appear.

4.4 Direct Control using Induced Sub-policies

We also provide videos11 1 https://sites.google.com/view/drdpaper of running sub-policies defined by πi=argmaxa𝔼(Zi). To clarify, the sub-policies are never rolled out during training or evaluation and are only used to compute Jdisentangij in Eq. 6. We execute these sub-policies and observe its difference with the main policy π=argmaxa𝔼(i=1MZi) to get a better visual effect of the reward decomposition. Take Seaquest in Figure 6(b) as an example, two sub-policies show distinctive preference. As Z1 mainly captures the reward for surviving and rescuing divers, π1 tends to stay close to the surface. However Z2 represents the return gained from shooting sharks, so π2 appears much more aggressive than π1. Also, without π1 we see that π2 dies quickly due to out of oxygen.

Figure 6: (a) Action statistics in an example trajectory of Seaquest. (b) Direct controlling using two induced sub-policies π1=argmaxaE(Z1),π2=argmaxaE(Z2): the top picture shows that π1 prefers to stay at the top to keep agent alive; the bottom picture shows that π2 prefers aggressive action of shooting sharks.

5 Related Work

Our method is closely related to previous work of reward decomposition. Reward function decomposition has been studied among others by russell2003q and sprague2003multiple. While these earlier works mainly focus on how to achieve optimal policy given the decomposed reward functions, there have been several recent works attempting to learn latent decomposed rewards. van2017hybrid construct an easy-to-learn value function by decomposing the reward function of the environment into n different reward functions. To ensure the learned decomposition is non-trivial, van2017hybrid proposed to split a state into different pieces following domain knowledge and then feed different state pieces into each reward function branch. While such method can accelerate learning process, it always requires many pre-defined preprocessing techniques. There has been other work that explores learn reward decomposition network end-to-end. grimm2019learning investigates how to learn independently-obtainable reward functions. While it learns interesting reward decomposition, their method requires that the environments be resettable to specific states since it needs multiple trajectories from the same starting state to compute their objective function. Besides, their method aims at learning different optimal policies for each decomposed reward function. Different from the works above, our method can learn meaningful implicit reward decomposition without any requirements on prior knowledge. Also, our method can leverage the decomposed sub-rewards to find better behaviour for a single agent.

Our work also relates to Horde (sutton2011horde). The Horde architecture consists of a large number of ‘sub-agents’ that learn in parallel via off-policy learning. Each demon trains a separate general value function (GVF) based on its own policy and pseudo-reward function. A pseudo-reward can be any feature-based signal that encodes useful information. The Horde architecture is focused on building up general knowledge about the world, encoded via a large number of GVFs. UVFA (schaul2015universal) extends Horde along a different direction that enables value function generalizing across different goals. Our method focuses on learning implicit reward decomposition in order to more efficiently learn a control policy.

6 Conclusion

In this paper, we propose Distributional Reward Decomposition for Reinforcement Learning (DRDRL), a novel reward decomposition algorithm which captures the multiple reward channel structure under distributional setting. Our algorithm significantly outperforms state-of-the-art RL methods RAINBOW on Atari games with multiple reward channels. We also provide interesting experimental analysis to get insight for our algorithm. In the future, we might try to develop reward decomposition method based on quantile networks (dabney2018implicit; dabney2018distributional).

Acknowledgments

This work was supported in part by the National Key Research & Development Plan of China (grant No. 2016YFA0602200 and 2017YFA0604500), and by Center for High Performance Computing and System Simulation, Pilot National Laboratory for Marine Science and Technology (Qingdao).

References

Appendix

Non-factorial model

To show that P(Zi=v)=P(Zi=v|Zj) is a reasonable assumption, we also implemented a non-factorial version of DRDRL, i.e. assuming that P(Zi=v)!=P(Zi=v|Zj), with 3 channels. The architecture of the non-factorial model is shown in figure 7. The three sub-distributions have K atoms, K2 atoms and K3 atoms respectively and they multiply to form a full distribution with K3 atoms. To maintain similar network capacity as C51, we set K=4 to form a full distribution of 64 atoms.

Figure 7: Model architecture of full distribution (non-factorial) method mentioned in section 3.1.

Ablative results of the non-factorial model are shown in the following section together with other tricks.

Ablative Analysis

Figure 8 shows how each component of our method contributes to DRDRL. Specifically, we test the performance of ‘1D convolution’, ‘1D convolution + KL’, ‘1D-convolution + onehot’ and ‘1D-convolution + KL + onehot’ as well as the non-factorial model.

One may argue that the superior performance of DRDRL might come from the fact that fitting N simple sub-distributions with M/N categories is easier than fitting one with M categories. We include another set of experiments to justify this argument by using M atoms instead of M/N atoms for each sub-distribution.

Figure 8: Training curves.

Results show that every part of DRDRL is important. To our surprise, the KL term not only allows us to perform reward decomposition, we also observe significant differences between curves with and without the KL term. This suggests that learning decomposed reward can greatly boost the performance of RL algorithms. Together with the degraded performance of M atoms instead of M/N atoms in each sub-distribution, it is suffice to suggest that the DRDRL’s success does not come from fitting easier distributions with M/N atoms.