Title: HIERARCHICAL REINFORCEMENT LEARNING VIA ADVANTAGE-WEIGHTED INFORMATION MAXIMIZATION

Abstract: Real-world tasks are often highly structured. Hierarchical reinforcement learning (HRL) has attracted research interest as an approach for leveraging the hierarchical structure of a given task in reinforcement learning (RL). However, identifying the hierarchical policy structure that enhances the performance of RL is not a trivial task. In this paper, we propose an HRL method that learns a latent variable of a hierarchical policy using mutual information maximization. Our approach can be interpreted as a way to learn a discrete and latent representation of the state-action space. To learn option policies that correspond to modes of the advantage function, we introduce advantage-weighted importance sampling. In our HRL method, the gating policy learns to select option policies based on an option-value function, and these option policies are optimized based on the deterministic policy gradient method. This framework is derived by leveraging the analogy between a monolithic policy in standard RL and a hierarchical policy in HRL by using a deterministic option policy. Experimental results indicate that our HRL approach can learn a diversity of options and that it can enhance the performance of RL in continuous control tasks.

1. INTRODUCTION:

Reinforcement learning (RL) has been successfully applied to a variety of tasks, including board games (Silver et al., 2016) , robotic manipulation tasks (Levine et al., 2016) , and video games (Mnih et al., 2015) . Hierarchical reinforcement learning (HRL) is a type of RL that leverages the hierarchical structure of a given task by learning a hierarchical policy (Sutton et al., 1999; Dietterich, 2000) . Past studies in this field have shown that HRL can solve challenging tasks in the video game domain (Vezhnevets et al., 2017; Bacon et al., 2017) and robotic manipulation (Daniel et al., 2016; Osa et al., 2018b) . In HRL, lower-level policies, which are often referred to as option policies, learn different behavior/control patterns, and the upper-level policy, which is often referred to as the gating policy, learns to select option policies. Recent studies have developed HRL methods using deep learning (Goodfellow et al., 2016) and have shown that HRL can yield impressive performance for complex tasks (Bacon et al., 2017; Frans et al., 2018; Vezhnevets et al., 2017; Haarnoja et al., 2018a) . However, identifying the hierarchical policy structure that yields efficient learning is not a trivial task, since the problem involves learning a sufficient variety of types of behavior to solve a given task.

In this study, we present an HRL method via the mutual information (MI) maximization with advantage-weighted importance, which we refer to as adInfoHRL. We formulate the problem of learning a latent variable in a hierarchical policy as one of learning discrete and interpretable repre-sentations of states and actions. Ideally, each option policy should be located at separate modes of the advantage function. To estimate the latent variable that corresponds to modes of the advantage function, we introduce advantage-weighted importance weights. Our approach can be considered to divide the state-action space based on an information maximization criterion, and it learns option policies corresponding to each region of the state-action space. We derive adInfoHRL as an HRL method based on deterministic option policies that are trained based on an extension of the deterministic policy gradient (Silver et al., 2014; Fujimoto et al., 2018) . The contributions of this paper are twofold:

1. We propose the learning of a latent variable of a hierarchical policy as a discrete and hidden representation of the state-action space. To learn option policies that correspond to the modes of the advantage function, we introduce advantage-weighted importance. 2. We propose an HRL method, where the option policies are optimized based on the deterministic policy gradient and the gating policy selects the option that maximizes the expected return. The experimental results show that our proposed method adInfoHRL can learn a diversity of options on continuous control tasks. Moreover, our approach can improve the performance of TD3 on such tasks as the Walker2d and Ant tasks in OpenAI Gym with MuJoco simulator.

2. BACKGROUND:

In this section, we formulate the problem of HRL in this paper and describe methods related to our proposal.

2.1. HIERARCHICAL REINFORCEMENT LEARNING:

We consider tasks that can be modeled as a Markov decision process (MDP), consisting of a state space S, an action space A, a reward function r : S × A → R, an initial state distribution ρ(s 0 ), and a transition probability p(s t+1 |s t , a t ) that defines the probability of transitioning from state s t and action a t at time t to next state s t+1 . The return is defined as R t = T i=t γ i-t r(s i , a i ), where γ is a discount factor, and policy π(a|s) is defined as the density of action a given state s. Let d π (s) = T t=0 γ t p(s t = s) denote the discounted visitation frequency induced by the policy π. The goal of reinforcement learning is to learn a policy that maximizes the expected return J(π) = E s0,a0,... [R 0 ] where s 0 ∼ ρ(s 0 ), a ∼ π and s t+1 ∼ p(s t+1 |s t , a t ). By defining the Q-function as Q π (s, a) = E s0,a0,... [R t |s t = s, a t = a], the objective function of reinforcement learning can be rewritten as follows:

EQUATION

Herein, we consider hierarchical policy π(a|s) = o∈O π(o|s)π(a|s, o), where o is the latent variable and O is the set of possible values of o. Many existing HRL methods employ a policy structure of this form (Frans et al., 2018; Vezhnevets et al., 2017; Bacon et al., 2017; Florensa et al., 2017; Daniel et al., 2016) . In general, latent variable o can be discrete (Frans et al., 2018; Bacon et al., 2017; Florensa et al., 2017; Daniel et al., 2016; Osa & Sugiyama, 2018) or continuous (Vezhnevets et al., 2017) . π(o|s) is often referred to as a gating policy (Daniel et al., 2016; Osa & Sugiyama, 2018) , policy over options (Bacon et al., 2017) , or manager (Vezhnevets et al., 2017) . Likewise, π(a|s, o) is often referred to as an option policy (Osa & Sugiyama, 2018) , sub-policy (Daniel et al., 2016) , or worker (Vezhnevets et al., 2017) . In HRL, the objective function is given by

EQUATION

As discussed in the literature on inverse RL (Ziebart, 2010) , multiple policies can yield equivalent expected returns. This indicates that there exist multiple solutions to latent variable o that maximizes the expected return. To obtain the preferable solution for o, we need to impose additional constraints in HRL. Although prior work has employed regularizers (Bacon et al., 2017) and constraints (Daniel et al., 2016) to obtain various option policies, the method of learning a good latent variable o that improves sample-efficiency of the learning process remains unclear. In this study we propose the learning of the latent variable by maximizing MI between latent variables and state-action pairs.

2.2. DETERMINISTIC POLICY GRADIENT:

The deterministic policy gradient (DPG) algorithm was developed for learning a monolithic deterministic policy µ θ (s) : S → A by Silver et al. (2014) . In off-policy RL, the objective is to maximize the expectation of the return, averaged over the state distribution induced by a behavior policy β(a|s):

EQUATION

When a policy is deterministic, the objective becomes J(π) = d β (s)Q π s, µ θ (s) ds. Silver et al. (2014) have shown that the gradient of a deterministic policy is given by

EQUATION

The DPG algorithm has been extended to the deep deterministic policy gradient (DDPG) for continuous control problems that require neural network policies (Lillicrap et al., 2016) . Twin Delayed Deep Deterministic policy gradient algorithm (TD3) proposed by Fujimoto et al. (2018) is a variant of DDPG that outperforms the state-of-the-art on-policy methods such as TRPO (Schulman et al., 2017a) and PPO (Schulman et al., 2017b) in certain domains. We extend this deterministic policy gradient to learn a hierarchical policy.

2.3. REPRESENTATION LEARNING VIA INFORMATION MAXIMIZATION:

Recent studies such as those by Chen et al. (2016) ; Hu et al. (2017) ; Li et al. (2017) have shown that an interpretable representation can be learned by maximizing MI. Given a dataset X = (x 1 , ..., x n ), regularized information maximization (RIM) proposed by Gomes et al. (2010) involves learning a conditional model p(y|x; η) with parameter vector η that predicts a label y. The objective of RIM is to minimize

EQUATION

where (η) is the regularization term, I η (x, y) is MI, and λ is a coefficient. MI can be decomposed as I η (x, y) = H(y) -H(y|x) where H(y) is entropy and H(y|x) the conditional entropy. Increasing H(y) conduces the label to be uniformly distributed, and decreasing H(y|x) conduces to clear cluster assignments. Although RIM was originally developed for unsupervised clustering problems, the concept is applicable to various problems that require learning a hidden discrete representation.

In this study, we formulate the problem of learning the latent variable o of a hierarchical policy as one of learning a latent representation of the state-action space.

3. LEARNING OPTIONS VIA ADVANTAGE-WEIGHTED INFORMATION MAXIMIZATION:

In this section, we propose a novel HRL method based on advantage-weighted information maximization. We first introduce the latent representation learning via advantage-weighted information maximization, and we then describe the HRL framework based on deterministic option policies.

3.1. LATENT REPRESENTATION LEARNING VIA ADVANTAGE-WEIGHTED INFORMATION MAXIMIZATION:

Although prior work has often considered H(o|s) or I(s, o), which results in a division of the state space, we are interested in using I (s, a), o for dividing the state-action space instead. A schematic sketch of our approach is shown in Figure 1 . As shown in the left side of Figure 1 , the advantage function often has multiple modes. Ideally, each option policies should correspond to separate modes of the advantage function. However, it is non-trivial to find the modes of the advantage function in practice. For this purpose, we reduce the problem of finding modes of the advantage function to that of finding the modes of the probability density of state action pairs.

We consider a policy based on the advantage function of the form Figure 1 : Schematic sketch of our HRL approach. By using the advantage-weighted importance, the problem of finding the modes of the advantage-function can be reduced to that of finding the modes of the density of state action pairs.

EQUATION

where

A π (s, a) = Q π (s, a) -V π (s) is the advantage function, V π (s)

is the state value function, and Z is the partition function. f (•) is a functional, which is a function of a function. f (•) is a monotonically increasing function with respect to the input variable and always satisfies f (•) > 0.

In our implementation we used the exponential function f (•) = exp(•). When following such a policy, an action with the larger advantage is drawn with a higher probability. Under this assumption, finding the modes of the advantage function is equivalent to finding modes of the density induced by π Ad . Thus, finding the modes of the advantage function can be reduced to the problem of clustering samples induced by π Ad .

Following the formulation of RIM introduced in Section 2.3, we formulate the problem of clustering samples induced by π Ad as the learning of discrete representations via MI maximization. For this purpose, we consider a neural network that estimates p(o|s, a; η) parameterized with vector η, which we refer to as the option network. We formulate the learning of the latent variable o as minimizing L option (η) = (η) -λI o, (s, a); η ,

where I(o, (s, a)) = Ĥ(o|s, a; η) -Ĥ(o; η), and (η) is the regularization term. In practice, we need to approximate the advantage function, and we learn the discrete variable o that corresponds to the modes of the current estimate of the advantage function. For regularization, we used a simplified version of virtual adversarial training (VAT) proposed by Miyato et al. (2016) . Namely, we set (η) = D KL p(o|s noise , a noise ; η)||p(o|s, a; η) where s noise = s + s , a noise = a + a , s and a denote white noise. This regularization term penalizes dissimilarity between an original state-action pair and a perturbed one, and Hu et al. (2017) empirically show that this regularization improves the performance of learning latent discrete representations.

When computing MI, we need to compute p(o) and H(o|s, a) given by 

EQUATION

Thus, the probability density of (s, a) induced by π Ad is necessary for computing MI for our purpose. To estimate the probability density of (s, a) induced by π Ad , we introduce the advantageweighted importance in the next section.

3.2. IMPORTANCE WEIGHTS FOR MUTUAL INFORMATION ESTIMATION:

Although we show that the problem of finding the modes of the advantage function can be reduced to MI maximization with respect to the samples induced by π Ad , samples induced by π Ad are not available in practice. While those induced during the learning process are available, a discrete representation obtained from such samples does not correspond to the modes of the advantage function.

To estimate the density induced by π Ad , we employ an importance sampling approach.

We assume that the change of the state distribution induced by the policy update is sufficiently small, namely, d πAd (s) ≈ d β (s). Then, the importance weight can be approximated as

EQUATION

)

and the normalized importance weight is given gy

EQUATION

)

As the partition function Z is canceled, we do not need to compute Z when computing the importance weight in practice. We call this importance weight W the advantage-weighted importance and employ it to compute the objective function used to estimate the latent variable.

This advantage-weighted importance is used to compute the entropy terms for computing MI in Equation ( 7). The empirical estimate of the entropy H(o) is given by

EQUATION

)

where the samples (s i , a i ) are drawn from p β (s, a) induced by a behavior policy β(a|s). Likewise, the empirical estimate of the conditional entropy H(o|s, a) is given by

EQUATION

The derivations of Equations ( 12) and ( 13) are provided in Appendix A. To train the option network, we store the samples collected by the M most recent behavior policies, to which we refer as onpolicy buffer D on . Although the algorithm works with entire samples stored in the replay buffer, we observe that the use of the on-policy buffer for latent representation learning exhibits better performance. For this reason, we decided to use the on-policy buffer in our implementation. Therefore, while the algorithm is off-policy in the sense that the option is learned from samples collected by behavior policies, our implementation is "semi"on-policy in the sense that we use samples collected by the most recent behavior policies.

4. HRL OBJECTIVE WITH DETERMINISTIC OPTION POLICIES:

Instead of stochastic option policies, we consider deterministic option policies and model them using separate neural networks. We denote by π(a|s, o) = µ o θ (s) deterministic option policies parameterized by vector θ. The objective function of off-policy HRL with deterministic option policies can then be obtained by replacing π(a|s) with o∈O π(o|s)π(a|s, o) in Equation (3):

EQUATION

)

where Q π (s, a; w) is an approximated Q-function parameterized using vector w. This form of the objective function is analogous to Equation (3). Thus, we can extend standard RL techniques to the learning of the gating policy π(o|s) in HRL with deterministic option policies.

In HRL, the goal of the gating policy is to generate a value of o that maximizes the conditional expectation of the return:

EQUATION

which is often referred to as the option-value function (Sutton et al., 1999) . When option policies are stochastic, it is often necessary to approximate the option-value function Q π Ω (s, o) in addition to the action-value function Q π (s, a). However, in our case, the option-value function for deterministic option policies is given by 

EQUATION

w target ← τ w +(1-τ )w target , θ target ← τ θ +(1-τ )θ target end if end for until the convergence return θ

which we can estimate using the deterministic option policy µ o θ (s) and the approximated actionvalue function Q π (s, a; w). In this work we employ the softmax gating policy of the form

EQUATION

)

which encodes the exploration in its form (Daniel et al., 2016) . The state value function is given as

EQUATION

which can be computed using Equation ( 17). We use this state-value function when computing the advantage-weighted importance as A(s, a) = Q(s, a) -V (s). In this study, the Q-function is trained in a manner proposed by Fujimoto et al. (2018) . Two neural networks (Q π w1 , Q π w2 ) are trained to estimate the Q-function, and the target value of the Q-function is computed as y i = r i + γ min 1,2 Q(s i , a i ) for sample (s i , a i , a i , r i ) in a batch sampled from a replay buffer, where r i = r(s i , a i ). In this study, the gating policy determines the option once every N time steps, i.e., t = 0, N, 2N, . . .

Neural networks that model µ o:

θ (a|s) for o = 1, ..., O, which we refer to as option-policy networks, are trained separately for each option. In the learning phase, p(o|s, a) is estimated by the option network. Then, samples are assigned to option o * = arg max o p(o|s, a; η) and are used to update the option-policy network that corresponds to o * . When performing a rollout, o is drawn by following the gating policy in Equation ( 17), and an action is generated by the selected option-policy network.

Differentiating the objective function in Equation ( 14), we obtain the deterministic policy gradient of our option-policy µ o θ (s) given by

EQUATION

The procedure of adInfoHRL is summarized by Algorithm 1. As in TD3 (Fujimoto et al., 2018) , we employed the soft update using a target value network and a target policy network. 

5. EXPERIMENTS:

We evaluated the proposed algorithm adInfoHRL on the OpenAI Gym platform (Brockman et al., 2016) with the MuJoCo Physics simulator (Todorov et al., 2012) . We compared its performance with that of PPO implemented in OpenAI baselines (Dhariwal et al., 2017) and TD3. Henderson et al. (2018) have recently claimed that algorithm performance varies across environment, there is thus no clearly best method for all benchmark environments, and off-policy and on-policy methods have advantages in different problem domains. To analyze the performance of adInfoHRL, we compared it with state-of-the-art algorithms for both on-policy and off-policy methods, although we focused on the comparison with TD3, as our implementation of adInfoHRL is based on it. To determine the effect of learning the latent variable via information maximization, we used the same network architectures for the actor and critic in adInfoHRL and TD3. In addition, to evaluate the benefit of the advantage-weighted importance, we evaluated a variant of adInfoHRL, which does not use the advantage-weighted importance for computing mutual information.We refer to this variant of adInfoHRL as infoHRL. The gating policy updated variable o once every three time steps. We tested the performance of adInfoHRL with two and four options.

The activation of options over time and snapshots of the learned option policies on the Walker2d task are shown in Figure 2 , which visualizes the result from adInfoHRL with four options. One can see that the option policies are activated in different phases of locomotion. While the option indicated by yellow in Figure 2 corresponds to the phase for kicking the floor, the option indicated by blue corresponds to the phase when the agent was on the fly. Visualization of the options learned on the HalfCheetah and Ant tasks are shown in Appendix D.

The averaged return of five trials is reported in Figure 3 (a)-(d). AdIfoHRL yields the best performance on Ant1 and Walker2d, whereas the performance of TD3 and adInfoHRL was comparable on HalfCheetah and Hopper, and PPO outperformed the other methods on Hopper. Henderson et al. (2018) claimed that on-policy methods show their superiority on tasks with unstable dynamics, and our experimental results are in line with such previous studies. AdinfoHRL outperformed infoHRL, which isthe variant of adInfoHRL without the advantage-weighted importance on all the tasks. This result shows that the adavatage-weighted importance enhanced the performance of learning options.

AdInfoHRL exhibited the sample efficiency on Ant and Walker2d in the sense that it required fewer samples than TD3 to achieve comparable performance on those tasks. The concept underlying ad-InfoHRL is to divide the state-action space to deal with the multi-modal advantage function and learn option policies corresponding to separate modes of the advantage function. Therefore, adIn-foHRL shows its superiority on tasks with the multi-modal advantage function and not on tasks with a simple advantage function. Thus, it is natural that the benefit of adInfoHRL is dependent on the characteristics of the task.

The outputs of the option network and the activation of options on Walker2d are shown in 

6. RELATED WORK AND DISCUSSION:

Past studies have proposed several ways to deal with the latent variable in HRL. The recent work by Smith et al. (2018) proposed inferred option policy gradients (IOPG), which is derived as an extension of policy gradient to the option framework. Nachum et al. (2018) recently proposed off-policy target correction for HRL on goal-oriented tasks, where a higher-level policy instructs a lower-level policy by generating the goal signal instead of an inferred latent variable. A popular approach for learning the latent variable in HRL is the variational approach. The recent work by Haarnoja et al. (2018a) is based on soft actor critic (Haarnoja et al., 2018b) , and the latent variable is inferred using the variational approach. The work by Hausman et al. (2018) is also closely related to the variational approach, and they proposed a method for learning a latent variable of a hierarchical policy via a variational bound. On the contrary, our method learns the latent variable by maximizing MI with advantage-weighted importance. Recent studies by Gregor et al. (2016) ; Florensa et al. (2017) ; Eysenbach et al. (2018) also considered the MI in their formulation. In these methods, MI between the state and the latent variable is considered so as to obtain diverse behaviors. Our approach is different from the previous studies in the sense that we employ MI between the latent variable and the state-action pairs, which leads to the division of the state-action space instead of considering only the state space. We think that dividing the state-action space is an efficient approach when the advantage function is multi-modal, as depicted in Figure 1 . InfoGAIL proposed by Li et al. (2017) learns the interpretable representation of the state-action space via MI maximization. InfoGAIL can be interpreted as a method that divides the state-action space based on the density induced by an expert's policy by maximizing the regularized MI objective. In this sense, it is closely related to our method, although their problem setting is imitation learning (Osa et al., 2018a) , which is different from our HRL problem setting.

The use of the importance weight based on the value function has appeared in previous studies (Dayan & Hinton, 1997; Kober & Peters, 2011; Neumann & Peters, 2009; Osa & Sugiyama, 2018) . For example, the method proposed by Neumann & Peters (2009) employs the importance weight based on the advantage function for learning a monolithic policy, while our method uses a similar importance weight for learning a latent variable of a hierarchical policy. Although Osa & Sugiyama (2018) proposed to learn a latent variable in HRL with importance sampling, their method is limited to episodic settings where only a single option is used in an episode.

Our method can be interpreted as an approach that divides the state-action space based on the MI criterion. This concept is related to that of Divide and Conquer (DnC) proposed by Ghosh et al. (2018) , although DnC clusters the initial states and does not consider switching between option policies during the execution of a single trajectory.

In this study we developed adInfoHRL based on deterministic option policies. However, the concept of dividing the state-action space via advantage-weighted importance can be applied to stochastic policy gradients as well. Further investigation in this direction is necessary in future work.

7. CONCLUSIONS:

We proposed a novel HRL method, hierarchical reinforcement learning via advantage-weighted information maximization. In our framework, the latent variable of a hierarchical policy is learned as a discrete latent representation of the state-action space. Our HRL framework is derived by considering deterministic option policies and by leveraging the analogy between the gating policy for HRL and a monolithic policy for the standard RL. The results of the experiments indicate that adInfoHRL can learn diverse options on continuous control tasks. Our results also suggested that our approach can improve the performance of TD3 in certain problem domains. 

C EXPERIMENTAL DETAILS:

We performed evaluations using benchmark tasks in the OpenAI Gym platform (Brockman et al., 2016) with Mujoco physics simulator (Todorov et al., 2012) . Hyperparameters of reinforcement learning methods used in the experiment are shown in Tables 1 2 3 . For exploration, both adInfoHRL and TD3 used the clipped noise drawn from the normal distribution as ∼ clip N (0, σ), -c, c , where σ = 0.2 and c = 0.5. For hyperparameters of PPO, we used the default values in OpenAI baselines (Dhariwal et al., 2017) . For the Walker2d, HalfCheetah, and Hopper tasks, we used the Walker2d-v1, HalfCHeetah-v1, and Hopper-v1 in the OpenAI Gym, respectively. For the Ant task, we used the AntEnv implemented in the rllab (Duan et al., 2016) . When training a policy with AdInfoHRL, infoHRL, and TD3, critics are trained once per time step, and actors are trained once every after two updates of the critics. The source code is available at https://github.com/ TakaOsa/adInfoHRL.

We performed the experiments five times with different seeds, and reported the averaged test return where the test return was computed once every 5000 time steps by executing 10 episodes without exploration. When performing the learned policy without exploration, the option was drawn as

EQUATION

instead of following the stochastic gating policy in Equations ( 17). 

D ADDITIONAL INFORMATION ON EXPERIMENTAL RESULTS:

On the HalfCheetah task, adInfoHRL delivered the best performance with two options. The distribution of options on HalfCheetah0v1 after one million steps is shown in Figure 4 . Although the state-action space is evenly divided, the options are not evenly activated. This behavior can occur because the state-action space is divided based on the density induced by the behavior policy while the activation of options is determined based on the quality of the option policies in a given state. Moreover, an even division in the action-state space is not necessarily the even division in the state space.

The activation of the options over time is shown in Figure 5 . It is clear that one of the option corresponds to the stable running phase and the other corresponds to the phase for recovering from unstable states. Since the Ant task implemented in rllab is known to be harder than the Ant-v1 implemented in the OpenAI gym, we reported the result of the Ant task in rllab in the main manuscript. Here, we report the result of the Ant-v1 task implemented in the OpenAI gym. On the Ant-v1 task, adInfoHRL yielded the best performance with two options. The performance of adInfoHRL with two options is comparable to that of TD3 on Ant-v1. This result indicates that the Ant-v1 task does not require a hierarchical policy structure, while a hierarchical policy improves the performance of learning on Ant-rllab. The distribution of options on Ant-v1 task after one million steps is shown in Figure 8 . The activation of the options over time is shown in Figure 9 . It is evident that two option policies on the Ant-v1 task corresponded to different postures of the agent.

A recent study on HRL by Smith et al. (2018) reported the performance of IOPG on Walker2d-v1, Hopper-v1, and HalfCheetah-v1. The study by Haarnoja et al. (2018a) reported the performance of SAC-LSP on Walker2d-v1, Hopper-v1, HalfCheetah-v1, and Ant-rllab. A comparison of performance between our method, IOPG, and SAC-LSP is summarized in Table 4 . We report the performance after 1 million steps. It is worth noting that adInfoHRL outperformed IOPG on these tasks in terms of the achieved return, although we are aware that the qualitative performance is also important in HRL. AdInfoHRL outperformed SAC-LSP on Walker2d-v1 and Ant-rllab, and SAC-LSP shows its superiority on HalfCheetah-v1 and Hopper-v1. However, the results of SAC-LSP were obtained by using reward scaling, which was not used in the evaluation of adInfoHRL. Therefore, further experiments are necessary for fair comparison under the same condition. Table 4 : Comparison of performance with existing methods after 1 million steps as reported in the literature. For adInfoHRL, we show the mean and the standard deviation of the results from the 10 different seeds. The performance of IOPG and SAC-LSP is from the original papers (Smith et al., 2018) and (Haarnoja et al., 

:

We report the result on the Ant task implemented in rllab(Duan et al.,

2016) instead of Ant-v1 implemented in the OpenAI gym, since the Ant task in the rllab is known to be harder than the Ant-v1 in the OpenAI gym. Results on Ant-v1 in the OpenAI gym is reported in Appendix D.

Figures and Tables:

) H(o|s, a) = E (s,a)∼p π Ad (s,a) [p(o|s, a; η) log p(o|s, a; η)] .

Figure 2: Activation of the four options over time steps on the Walker2d task.

Figure  3(e)-(f), which visualize the result from adInfoHRL with four options. For visualization, the dimensionality was reduced using t-SNE (van derMaaten & Hinton, 2008). The state-action space

Output of the option network in the state-action space on Walker2d. Activation of options in the state space on Walker2d.

Figure 3: Performance of adInfoHRL. (a)-(d) show comparison with baseline methods. (e) and (f) show the output of the option network and the activation of options on Walker2d, respectively.

Output of the option network in the stateaction space on HalfCheetah-v1. Activation of options in the state space on HalfCheetah-v1.

Figure 4: Distribution of options on the HalfCheetah-v1 task using adInfoHRL with two options. The dimensionality is reduced by t-SNE for visualization.

Figure 5: Activation of options over time steps on the HalfCheetah-v1 task using adInfoHRL with two options.

Output of the option network in the stateaction space on the Ant-rllab task.

Activation of options in the state space on the Ant-rllab task.

Figure 6: Distribution of options on Ant-rllab task using adInfoHRL with four options. The dimensionality is reduced by t-SNE for visualization.

Figure 7: Activation of the options over time steps on Ant-rllab task. Four options are learned.

Averaged return on the Ant-v1 task in the OpenAI gym.

Output of the option network in the state-action space on Ant-v1.

Activation of options in the state space on Ant-v1.

Figure 8: Distribution of options on the Ant-v1 task using adInfoHRL with two options. The dimensionality is reduced by t-SNE for visualization.

Figure 9: Activation of options over time steps on the Ant-v1 task using adInfoHRL with two options.

Algorithm 1 HRL via Advantage-Weighted Information Maximization (adInfoHRL) Input: Number of options O, size of on-policy buffer Initialize: Replay buffer D R , on-policy buffer D on , network parameters η, θ, w, θ target , w target repeat for t = 0 to t = T do Draw an option for a given s by following Equation 17: o ∼ π(o|s) Draw an action a ∼ β(a|s, o) = µ o θ (s) + Record a data sample (s, a, r, s ) Aggregate the data in D R and D on if the on-policy buffer is full then Update the option network by minimizing Equation (7) for samples in D on Clear the on-policy buffer D on end if Sample a batch D batch ∈ D R Update the Q network parameter w if t mod d then Estimate p(o|s i , a i ) for (s i , a i ) ∈ D batch using the option network Assign samples (s i , a i ) ∈ D batch to the option o

Table content: <table/>

Hyperparameters of adInfoHRL used in the experiment.

Table content: <table><tr><td>Description</td><td colspan="2">Symbol Value</td></tr><tr><td>Coefficient for updating the target network</td><td>τ</td><td>0.005</td></tr><tr><td>Discount factor</td><td>γ</td><td>0.99</td></tr><tr><td>Learning rate for actor</td><td/><td>0.001</td></tr><tr><td>Learning rate for critic</td><td/><td>0.001</td></tr><tr><td>Batch size for critic</td><td/><td>100</td></tr><tr><td>Total batch size for all option policies</td><td/><td>200 (option num=2), 400 (option num=4)</td></tr><tr><td>Batch size for the option network</td><td/><td>50</td></tr><tr><td>Size of the on-policy buffer</td><td/><td>5000</td></tr><tr><td>Number of epochs for training the option network</td><td/><td>40</td></tr><tr><td>Number of units in hidden layers</td><td/><td>(400, 300)</td></tr><tr><td>Activation function</td><td/><td>Relu, Relu, tanh</td></tr><tr><td>optimizer</td><td/><td>Adam</td></tr><tr><td>noise clip threshold</td><td>c</td><td>0.5</td></tr><tr><td>noise for exploration</td><td/><td>0.1</td></tr><tr><td>action noise for the critic update</td><td/><td>0.2</td></tr><tr><td>variance of the noise for MI regularization</td><td/><td>0.04</td></tr><tr><td>coefficient for the MI term</td><td>λ</td><td>0.1</td></tr></table>

Hyperparameters of TD3 used in the experiment.

Table content: <table><tr><td>Description</td><td colspan="2">Symbol Value</td></tr><tr><td>Coefficient for updating the target network</td><td>τ</td><td>0.005</td></tr><tr><td>Discount factor</td><td>γ</td><td>0.99</td></tr><tr><td>Learning rate for actor</td><td/><td>0.001</td></tr><tr><td>Learning rate for critic</td><td/><td>0.001</td></tr><tr><td>Batch size</td><td/><td>100</td></tr><tr><td>Number of units in hidden layers</td><td/><td>(400, 300)</td></tr><tr><td>Activation function</td><td/><td>Relu, Relu, tanh</td></tr><tr><td>optimizer</td><td/><td>Adam</td></tr><tr><td>noise clip threshold</td><td>c</td><td>0.5</td></tr><tr><td>noise for exploration</td><td/><td>0.1</td></tr><tr><td>action noise for the critic update</td><td/><td>0.2</td></tr></table>

Hyperparameters of PPO used in the experiment. We tuned hyperparameters for our tasks, which are defferent from the default parameters in OpenAI baselines(Dhariwal et al., 2017).

Table content: <table><tr><td>Description</td><td colspan="2">Symbol Value</td></tr><tr><td colspan="2">Coefficient for updating the target network τ</td><td>0.001</td></tr><tr><td>Discount factor</td><td>γ</td><td>0.99</td></tr><tr><td>Batch size</td><td/><td>2048</td></tr><tr><td>Number of units in hidden layers</td><td/><td>(64, 64)</td></tr><tr><td>Clipping parameter</td><td/><td>0.15</td></tr><tr><td>Initial learning rate</td><td/><td>0.0005</td></tr><tr><td>Learning rate schedule</td><td/><td>linear</td></tr></table>

Peer review (iclr2019):
 The authors propose an HRL algorithm that attempts to learn options that maximize their mutual information with the state-action density under the optimal policy. Several key terms are used in ways that differ from the rest of the literature. The authors claim options are learned in an "unsupervised" manner, but it is unclear what this means. Previous work (none of which is cited) has dealt with unsupervised option discovery in the context of mutual information maximization (Variational intrinsic control, diversity is all you need, etc), but they do so in the absence of reward, unlike this paper. "Optimal policy" is similarly abused, with it appearing to mean optimal from the perspective of the current model parameters, rather than optimal in any global sense. Or at least I think that is what the authors intend. If they do mean the globally optimal policy, then its unclear how to interpret Equation 8, with its reference to a behavior policy and an advantage function, neither of which would be available if meant to represent the global optimum. Equation 10 comes out of nowhere. One must assume they meant "maximize mutual information" and not "minimize", but who knows. Why is white-noise being added to the states and actions? Is this some sort of noise-contrastive estimation approach to mutual information estimation? It doesn*t appear to be, but it is unclear what else could motivate it. Even the appendices fail to shine light on this equation. The algorithm block isn*t terribly helpful. The "t" variable is used outside of its for loop, which draws into question the exact nesting structure of the underlying algorithm (which isn*t obvious for HRL methods). There aren*t any equations referenced, with the option policy network*s update not even referencing the loss nor data over which the loss would be evaluated. Some of the experimental results show promise, but the PPO Ant result raises some questions. Clearly the OpenAI implementation of PPO used would have tuned for the OpenAI gym Ant implementation, and the appendix shows it getting decent results. But it never takes off in the harder RlLab version -- were the hyper-parameters adjusted for this new environment? It is also odd that no other HRL approaches are evaluated against, given the number cited. Running these methods might be too costly, but surely a table comparing results reported in those papers should be included. A minor point: another good baseline would be TD3 with the action repeat adjusted to be inline with the gating policy. I apologise if this review came off as too harsh -- I believe a good paper can be made of this with extensive rewrites and additional experiments. But the complete lack of clarity makes it feel like it was rushed out prematurely. EDIT: Now this is a paper that makes sense! With the terminology cleared up and the algorithm fully unpacked, this approach seems quite interesting. The experimental results could always be stronger, but no longer have any holes in them. Score 3-->6