# modelaugmented_actorcritic_backpropagating_through_paths__5d9b00a6.pdf

Under review as a conference paper at ICLR 2020

MODEL-AUGMENTED ACTOR-CRITIC: BACKPROPAGATING THROUGH PATHS

Anonymous authors Paper under double-blind review

Current model-based reinforcement learning approaches use the model simply as a learned black-box simulator to augment the data for policy optimization or value function learning. In this paper, we show how to make more effective use of the model by exploiting its differentiability. We construct a policy optimization algorithm that uses the pathwise derivative of the learned model and policy across future timesteps. Instabilities of learning across many timesteps are prevented by using a terminal value function, learning the policy in an actor-critic fashion. Furthermore, we present a derivation on the monotonic improvement of our objective in terms of the gradient error in the model and value function. We show that our approach (i) is consistently more sample efficient than existing state-of-the-art model-based algorithms, (ii) matches the asymptotic performance of model-free algorithms, and (iii) scales to long horizons, a regime where typically past model-based approaches have struggled.

1 INTRODUCTION

Model-based reinforcement learning (RL) offers the potential to be a general-purpose tool for learning complex policies while being sample efficient. When learning in real-world physical systems, data collection can be an arduous process. Contrary to model-free methods, model-based approaches are appealing due to their comparatively fast learning. By first learning the dynamics of the system in a supervised learning way, it can exploit off-policy data. Then, model-based methods use the model to derive controllers from it either parametric controllers (Luo et al., 2019; Buckman et al., 2018; Janner et al., 2019) or non-parametric controllers (Nagabandi et al., 2017; Chua et al., 2018).

Current model-based methods learn with an order of magnitude less data than their model-free counterparts while achieving the same asymptotic convergence. Tools like ensembles, probabilistic models, planning over shorter horizons, and meta-learning have been used to achieved such performance (Kurutach et al., 2018; Chua et al., 2018; Clavera et al., 2018). However, the model usage in all of these methods is the same: simple data augmentation. They use the learned model as a black-box simulator generating samples from it. In high-dimensional environments or environments that require longer planning, substantial sampling is needed to provide meaningful signal for the policy. Can we further exploit our learned models?

In this work, we propose to estimate the policy gradient by backpropagating its gradient through the model using the pathwise derivative estimator. Since the learned model is differentiable, one can link together the model, reward function, and policy to obtain an analytic expression for the gradient of the returns with respect to the policy. By computing the gradient in this manner, we obtain an expressive signal that allows rapid policy learning. We avoid the instabilities that often result from back-propagating through long horizons by using a terminal Q-function. This scheme fully exploits the learned model without harming the learning stability seen in previous approaches (Kurutach et al., 2018; Heess et al., 2015). The horizon at which we apply the terminal Q-function acts as a hyperparameter between model-free (when fully relying on the Q-function) and model-based (when using a longer horizon) of our algorithm.

The main contribution of this work is a model-based method that significantly reduces the sample complexity compared to state-of-the-art model-based algorithms (Janner et al., 2019; Buckman et al., 2018). For instance, we achieve a 10k return in the half-cheetah environment in just 50 trajectories. We theoretically justify our optimization objective and derive the monotonic improvement of our

Under review as a conference paper at ICLR 2020

learned policy in terms of the Q-function and the model error. Furtermore, we experimentally analyze the theoretical derivations. Finally, we pinpoint the importance of our objective by ablating all the components of our algorithm. The results are reported in four model-based benchmarking environments (Wang et al., 2019; Todorov et al., 2012). The low sample complexity and high performance of our method carry high promise towards learning directly on real robots.

2 RELATED WORK

Model-Based Reinforcement Learning. Learned dynamics models offer the possibility to reduce sample complexity while maintaining the asymptotic performance. For instance, the models can act as a learned simulator on which a model-free policy is trained on (Kurutach et al., 2018; Luo et al., 2019; Janner et al., 2019). The model can also be used to improve the target value estimates (Feinberg et al., 2018) or to provide additional context to a policy (Du & Narasimhan, 2019). Contrary to these methods, our approach uses the model in a different way: we exploit the fact that the learned simulator is differentiable and optimize the policy with the analytical gradient. Long term predictions suffer from a compounding error effect in the model, resulting in unrealistic predictions. In such cases, the policy tends to overfit to the deficiencies of the model, which translates to poor performance in the real environment; this problem is known as model-bias (Deisenroth & Rasmussen, 2011). The model-bias problem has motivated work that uses meta-learning (Clavera et al., 2018), interpolation between different horizon predictions (Buckman et al., 2018; Janner et al., 2019), and interpolating between model and real data (Kalweit & Boedecker, 2017). To prevent model-bias, we exploit the model for a short horizon and use a terminal value function to model the rest of the trajectory. Finally, since our approach returns a stochastic policy, dynamics model, and value function could use model-predictive control (MPC) for better performance at test time, similar to (Lowrey et al., 2018; Hong et al., 2019). MPC methods (Nagabandi et al., 2017) have shown to be very effective when the uncertainty of the dynamics is modelled (Chua et al., 2018; Wang & Ba, 2019).

Differentable Planning. Previous work has used backpropagate through learned models to obtain the optimal sequences of actions. For instance, Levine & Abbeel (2014) learn linear local models and obtain the optimal sequences of actions, which is then distilled into a neural network policy. The planning can be incorporated into the neural network architecture (Okada et al., 2017; Tamar et al., 2016; Srinivas et al., 2018; Karkus et al., 2019) or formulated as a differentiable function (Pereira et al., 2018; Amos et al., 2018). Planning sequences of actions, even when doing model-predictive control (MPC), does not scale well to high-dimensional, complex domains Janner et al. (2019). Our method, instead learns a neural network policy in an actor-critic fashion aided with a learned model. In our study, we evaluate the benefit of carrying out MPC on top of our learned policy at test time, Section 5.4. The results suggest that the policy captures the optimal sequence of action, and re-planning does not result in significant benefits.

Policy Gradient Estimation. The reinforcement learning objective involves computing the gradient of an expectation (Schulman et al., 2015a). By using Gaussian processes (Deisenroth & Rasmussen, 2011), it is possible to compute the expectation analytically. However, when learning expressive parametric non-linear dynamical models and policies, such closed form solutions do not exist. The gradient is then estimated using Monte-Carlo methods (Mohamed et al., 2019). In the context of model-based RL, previous approaches mostly made use of the score-function, or REINFORCE estimator (Peters & Schaal, 2006; Kurutach et al., 2018). However, this estimator has high variance and extensive sampling is needed, which hampers its applicability in high-dimensional environments. In this work, we make use of the pathwise derivative estimator (Mohamed et al., 2019). Similar to our approach, Heess et al. (2015) uses this estimator in the context of model-based RL. However, they just make use of real-world trajectories that introduces the need of a likelihood ratio term for the model predictions, which in turn increases the variance of the gradient estimate. Instead, we entirely rely on the predictions of the model, removing the need of likelihood ratio terms.

Actor-Critic Methods. Actor-critic methods alternate between policy evaluation, computing the value function for the policy; and policy improvement using such value function (Sutton & Barto, 1998; Barto et al., 1983). Actor-critic methods can be classified between on-policy and off-policy. On-policy methods tend to be more stable, but at the cost of sample efficiency (Sutton, 1991; Mnih et al., 2016). On the other hand, off-policy methods offer better sample complexity (Lillicrap et al., 2015). Recent work has significantly stabilized and improved the performance of off-policy methods using maximum-entropy objectives (Haarnoja et al., 2018a) and multiple value functions (Fujimoto

Under review as a conference paper at ICLR 2020

et al., 2018). Our method combines the benefit of both. By using the learned model we can have a learning that resembles an on-policy method while still being off-policy.

3 BACKGROUND

In this section, we present the reinforcement learning problem, two different lines of algorithms that tackle it, and a summary on Monte-Carlo gradient estimators.

3.1 REINFORCEMENT LEARNING

A discrete-time finite Markov decision process (MDP) M is defined by the tuple (S, A, f, r, γ, p0, T). Here, S is the set of states, A the action space, st+1 f(st, at) the transition distribution, r : S A R is a reward function, p0 : S R+ represents the initial state distribution, γ the discount factor, and T is the horizon of the process. We define the return as the sum of rewards r(st, at) along a trajectory τ := (s0, a0, ..., s T 1, a T 1, s T ). The goal of reinforcement learning is to find a policy πθ : S A R+ that maximizes the expected return, i.e., maxθ J(θ) = maxθ E[P

t γtr(st, at)].

Actor-Critic. In actor-critic methods, we learn a function ˆQ (critic) that approximates the expected return conditioned on a state s and action a, E[P

t γtr(st, at)|s0 = s, a0 = a]. Then, the learned Q-function is used to optimize a policy π (actor). Usually, the Q-function is learned by iteratively minimizing the Bellman residual:

JQ = E[( ˆQ(st, at) (r(st, at) + γ ˆQ(st+1, at+1)))2]

The above method is referred as one-step Q-learning, and while a naive implementation often results in unstable behaviour, recent methods have succeeded in stabilizing the Q-function training (Fujimoto et al., 2018). The actor then can be trained to maximize the learned ˆQ function Jπ = E h ˆQ(s, π(s)) i . The benefit of this form of actor-critic method is that it can be applied in an off-policy fashion, sampling random mini-batches of transitions from an experience replay buffer (Lin, 1992).

Model-Based RL. Model-based methods, contrary to model-free RL, learn the transition distribution from experience. Typically, this is carried out by learning a parametric function approximator ˆfφ, known as a dynamics model. We define the state predicted by the dynamics model as ˆst+1, i.e., ˆst+1 ˆfφ(st, at). The models are trained via maximum likelihood: maxφ Jf(φ) = maxφ E[log p(ˆst+1|st, at)]

3.2 MONTE-CARLO GRADIENT ESTIMATORS

In order to optimize the reinforcement learning objective, it is needed to take the gradient of an expectation. In general, it is not possible to compute the exact expectation so Monte-Carlo gradient estimators are used instead. These are mainly categorized into three classes: the pathwise, score function, and measure-valued gradient estimator (Mohamed et al., 2019). In this work, we use the pathwise gradient estimator, which is also known as the re-parameterization trick (Kingma & Welling, 2013). This estimator is derived from the law of the unconscious statistician (LOTUS) (Grimmett & Stirzaker, 2001)

Epθ(x)[f(x)] = Ep(ϵ)[f(gθ(ϵ)]

Here, we have stated that we can compute the expectation of a random variable x without knowing its distribution, if we know its corresponding sampling path and base distribution. A common case, and the one used in this manuscript, θ parameterizes a Gaussian distribution: x pθ = N(µθ, σ2 θ), which is equivalent to x = µθ + ϵσθ for ϵ N(0, 1).

4 POLICY GRADIENT VIA MODEL-AUGMENTED PATHWISE DERIVATIVE

Exploiting the full capability of learned models has the potential to enable complex and highdimensional real robotics tasks while maintaining low sample complexity. Our approach, modelaugmented actor-critic (MAAC), exploits the learned model by computing the analytic gradient

Under review as a conference paper at ICLR 2020

of the returns with respect to the policy. In contrast to sample-based methods, which one can think of as providing directional derivatives in trajectory space, MAAC computes the full gradient, providing a strong learning signal for policy learning, which further decreases the sample complexity.

r0 r1 r2 r H Q

Figure 1: Stochastic computation graph of the proposed objective Jπ. The stochastic nodes are represented by circles and the deterministic ones by squares.

In the following, we present our policy optimization scheme and describe the full algorithm.

4.1 MODEL-AUGMENTED ACTOR-CRITIC OBJECTIVE

Among model-free methods, actor-critic methods have shown superior performance in terms of sample efficiency and asymptotic performance (Haarnoja et al., 2018a). However, their sample efficiency remains worse than modelbased approaches, and fully off-policy methods still show instabilities comparing to on-policy algorithms (Mnih et al., 2016). Here, we propose a modification of the Q-function parametrization by using the model predictions on the first time-steps after the action is taken. Specifically, we do policy optimization by maximizing the following objective:

t=0 γtr(st) + γH ˆQ(s H, a H)

whereby, st+1 ˆf(st, at) and at πθ(st). Note that under the true dynamics and Q-function, this objective is the same as the RL objective. Contrary to previous reinforcement learning methods, we optimize this objective by back-propagation through time. Since the learned dynamics model and policy are parameterized as Gaussian distributions, we can make use of the pathwise derivative estimator to compute the gradient, resulting in an objective that captures uncertainty while presenting low variance. The computational graph of the proposed objective is shown in Figure 1.

While the proposed objective resembles n-step bootstrap (Sutton & Barto, 1998), our model usage fundamentally differs from previous approaches. First, we do not compromise between being offpolicy and stability. Typically, n-step bootstrap is either on-policy, which harms the sample complexity, or its gradient estimation uses likelihood ratios, which presents large variance and results in unstable learning (Heess et al., 2015). Second, we obtain a strong learning signal by backpropagating the gradient of the policy across multiple steps using the pathwise derivative estimator, instead of the REINFORCE estimator (Mohamed et al., 2019; Peters & Schaal, 2006). And finally, we prevent the exploding and vanishing gradients effect inherent to back-propagation through time by the means of the terminal Q-function (Kurutach et al., 2018).

The horizon H in our proposed objective allows us to trade off between the accuracy of our learned model and the accuracy of our learned Q-function. Hence, it controls the degree to which our algorithm is model-based or well model-free. If we were not to trust our model at all (H = 0), we would end up with a model-free update; for H = , the objective results in a shooting objective. Note that we will perform policy optimization by taking derivatives of the objective, hence we require accuracy on the derivatives of the objective and not on its value. The following lemma provides a bound on the gradient error in terms of the error on the derivatives of the model, the Q-function, and the horizon H. Lemma 4.1 (Gradient Error). Let ˆf and ˆQ be the learned approximation of the dynamics f and Q-function Q, respectively. Assume that Q and ˆQ have Lq/2-Lipschitz continuous gradient and f and ˆf have Lf/2-Lipschitz continuous gradient. Let ϵf = maxt ˆf(ˆst, ˆat) f(st, at) 2 be the error on the model derivatives and ϵQ = ˆQ(ˆs H, ˆa H) Q(s H, a H) 2 the error on the Q-function derivative. Then the error on the gradient between the learned objective and the true objective can be bounded by:

E h θJπ θ ˆJπ 2 i c1(H)ϵf + c2(H)ϵQ

Proof. See Appendix.

Under review as a conference paper at ICLR 2020

The result in Lemma 4.1 stipulates the error of the policy gradient in terms of the maximum error in the model derivatives and the error in the Q derivatives. The functions c1 and c2 are functions of the horizon and depend on the Lipschitz constants of the model and the Q-function. Note that we are just interested in the relation between both sources of error, since the gradient magnitude will be scaled by the learning rate, or by the optimizer, when applying it to the weights.

4.2 MONOTONIC IMPROVEMENT

In the previous section, we presented our objective and the error it incurs in the policy gradient with respect to approximation error in the model and the Q function. However, the error on the gradient is not indicative of the effect of the desired metric: the average return. Here, we quantify the effect of the modeling error on the return. First, we will bound the KL-divergence between the policies resulting from taking the gradient with the true objective and the approximated one. Then we will bound the performance in terms of the KL.

Lemma 4.2 (Total Variation Bound). Under the assumptions of the Lemma 4.1, let θ = θo +α θJπ be the parameters resulting from taking a gradient step on the exact objective, and ˆθ = θo + α θ ˆJπ the parameters resulting from taking a gradient step on approximated objective, where α R+. Then the following bound on the total variation distance holds

max s DTV(πθ||πˆθ) αc3(ϵfc1(H) + ϵQc2(H))

Proof. See Appendix.

The previous lemma results in a bound on the distance between the policies originated from taking a gradient step using the true dynamics and Q-function, and using its learned counterparts. Now, we can derive a similar result from Kakade & Langford (2002) to bound the difference in average returns.

Theorem 4.1 (Monotonic Improvement). Under the assumptions of the Lemma 4.1, be θ and ˆθ as defined in Lemma 4.2, and assuming that the reward is bounded by rmax. Then the average return of the πˆθ satisfies

Jπ(ˆθ) Jπ(θ) 2αrmax

1 γ αc3(ϵfc1(H) + ϵQc2(H))

Proof. See Appendix.

Hence, we can provide explicit lower bounds of improvement in terms of model error and function error. Theorem 4.1 extends previous work of monotonic improvement for model-free policies (Schulman et al., 2015b; Kakade & Langford, 2002), to the model-based and actor critic set up by taking the error on the learned functions into account. From this bound one could, in principle, derive the optimal horizon H that minimizes the gradient error. However, in practice, approximation errors are hard to determine and we treat H as an extra hyper-parameter. In section 5.2, we experimentally analyze the error on the gradient for different estimators and values of H.

4.3 ALGORITHM

Based on the previous sections, we develop a new algorithm that explicitly optimizes the modelaugmented actor-critic (MAAC) objective. The overall algorithm is divided into three main steps: model learning, policy optimization, and Q-function learning.

Model learning. In order to prevent overfitting and overcome model-bias (Deisenroth & Rasmussen, 2011), we use a bootstrap ensemble of dynamics models { ˆfφ1, ..., ˆfφM }. Each of the dynamics models parameterizes the mean and the covariance of a Gaussian distribution with diagonal covariance. The bootstrap ensemble captures the epistemic uncertainty, uncertainty due to the limited capacity or data, while the probabilistic models are able to capture the aleatoric uncertainty (Chua et al., 2018), inherent uncertainty of the environment. We denote by fφ the transitions dynamics resulting from φU, where U U[M] is uniform random variable on {1, ..., M}. The dynamics models are trained via maximum likelihood with early stopping on a validation set.

Under review as a conference paper at ICLR 2020

Algorithm 1 MAAC

1: Initialize the policy πθ, model ˆfφ, ˆQψ, Denv , and the model dataset Dmodel 2: repeat 3: Sample trajectories from the real environment with policy πθ. Add them to Denv. 4: for i = 1 . . . G1 do 5: φ φ βf φJf(φ) using data from Denv. 6: end for 7: Sample trajectories T from ˆfφ. Add them to Dmodel. 8: D Dmodel Denv 9: for i = 1 . . . G2 do 10: Update θ θ + βπ θJπ(θ) using data from D 11: Update ψ ψ βQ ψJQ(ψ) using data from D 12: end for 13: until the policy performs well in the real environment 14: return Optimal parameters θ

Policy Optimization. We extend the MAAC objective with an entropy bonus (Haarnoja et al., 2018b), and perform policy learning by maximizing

t=0 γtr(ˆst) + γHQψ(ˆs H, a H)

where ˆst+1 fφ(ˆst, at), ˆs0 D, a πθ. We learn the policy by using the pathwise derivative of the model through H steps and the Q-function by sampling multiple trajectories from the same ˆs0. Hence, we learn a maximum entropy policy using pathwise derivative of the model through H steps and the Q-function. We compute the expectation by sampling multiple actions and states from the policy and learned dynamics, respectively.

Q-function Learning. In practice, we train two Q-functions (Fujimoto et al., 2018) since it has been experimentally proven to yield better results. We train both Q functions by minimizing the Bellman error (Section 3.1):

JQ(ψ) = E[(Qψ(st, at) (r(st, at) + γQψ(st+1, at+1)))2]

Similar to (Janner et al., 2019), we minimize the Bellman residual on states previously visited and imagined states obtained from unrolling the learned model. Finally, the value targets are obtained in the same fashion the Stochastic Ensemble Value Expansion (Buckman et al., 2018), using H as a horizon for the expansion. In doing so, we maximally make use of the model by not only using it for the policy gradient step, but also for training the Q-function.

Our method, MAAC, iterates between collecting samples from the environment, model training, policy optimization, and Q-function learning. A practical implementation of our method is described in Algorithm 1. First, we obtain trajectories from the real environment using the latest policy available. Those samples are appended to a replay buffer Denv, on which the dynamics models are trained until convergence. The third step is to collect imaginary data from the models: we collect k-step transitions by unrolling the latest policy from a randomly sampled state on the replay buffer. The imaginary data constitutes the Dmodel, which together with the replay buffer, is used to learn the Q-function and train the policy.

Our algorithm consolidates the insights built through the course of this paper, while at the same time making maximal use of recently developed actor-critic and model-based methods. All in all, it consistently outperforms previous model-based and actor-critic methods.

Our experimental evaluation aims to examine the following questions: 1) How does MAAC compares against state-of-the-art model-based and model-free methods? 2) Does the gradient error correlate

Under review as a conference paper at ICLR 2020

with the derived bound?, 3) Which are the key components of its performance?, and 4) Does it benefit from planning at test time?

In order to answer the posed questions, we evaluate our approach on model-based continuous control benchmark tasks in the Mu Jo Co simulator (Todorov et al., 2012; Wang et al., 2019).

5.1 COMPARISON AGAINST STATE-OF-THE-ART

We compare our method on sample complexity and asymptotic performance against state-of-the-art model-free (MF) and model-based (MB) baselines. Specifically, we compare against the model-free soft actor-critic (SAC) (Haarnoja et al., 2018a), which is an off-policy algorithm that has been proven to be sample efficient and performant; as well as two state-of-the-art model-based baselines: modelbased policy-optimization (MBPO) (Janner et al., 2019) and stochastic ensemble value expansion (STEVE) (Buckman et al., 2018). The original STEVE algorithm builds on top of the model-free algorithm DDPG (Lillicrap et al., 2015), however this algorithm is outperformed by SAC. In order to remove confounding effects of the underlying model-free algorithm, we have implemented the STEVE algorithm on top of SAC. We also add SVG(1) (Heess et al., 2015) to comparison, which similar to our method uses the derivative of dynamic models to learn the policy.

The results, shown in Fig. 2, highlight the strength of MAAC in terms of performance and sample complexity. MAAC scales to higher dimensional tasks while maintaining its sample efficiency and asymptotic performance. In all the four environments, our method learns faster than previous MB and MF methods. We are able to learn near-optimal policies in the half-cheetah environment in just over 50 rollouts, while previous model-based methods need at least the double amount of data. Furthermore, in complex environments, such as ant, MAAC achieves near-optimal performance within 150 rollouts while other take orders of magnitudes more data.

Figure 2: Comparison against state-of-the-art model-free and model-based baselines in four different Mu Jo Co environments. Our method, MAAC, attains better sample efficiency and asymptotic performance than previous approaches. The gap in performance between MAAC and previous work increases as the environments increase in complexity.

5.2 GRADIENT ERROR

Here, we investigate how the bounds obtained relate to the empirical performance. In particular, we study the effect of the horizon of the model predictions on the gradient error. In order to do so, we

Under review as a conference paper at ICLR 2020

construct a double integrator environment; since the transitions are linear and the cost is quadratic for a linear policy, we can obtain the analytic gradient of the expect return.

Figure 3: L1 error on the policy gradient when using the proposed objective for different values of the horizon H as well as the error obtained when using the true dynamics. The results correlate with the assumption that the error in the learned dynamics is lower than the error in the Q-function, as well as they correlate with the derived bounds.

Figure 3 depicts the L1 error of the MAAC objective for different values of the horizon H as well as what would be the error using the true dynamics. As expected, using the true dynamics yields to lower gradient error since the only source comes from the learned Q-function that is weighted down by γH. The results using learned dynamics correlate with our assumptions and the derived bounds: the error from the learned dynamics is lower than the one in the Q-funtion, but it scales poorly with the horizon. For short horizons the error decreases as we increase the horizon. However, large horizons is detrimental since it magnifies the error on the models.

5.3 ABLATIONS

In order to investigate the importance of each of the components of our overall algorithm, we carry out an ablation test. Specifically, we test three different components: 1) not using the model to train the policy, i.e., set H = 0, 2) not using the STEVE targets for training the critic, and 3) using a single sample estimate of the path-wise derivative.

The ablation test is shown in Figure 4. The test underpins the importance of backpropagating through the model: setting H to be 0 inflicts a severe drop in the algorithm performance. On the other hand, using the STEVE targets results in slightly more stable training, but it does not have a significant effect. Finally, while single sample estimates can be used in simple environments, they are not accurate enough in higher dimensional environments such as ant.

Figure 4: Ablation test of our method. We test the importance of several components of our method: not using the model to train the policy (H = 0), not using the STEVE targets for training the Q-function (-STEVE), and using a single sample estimate of the pathwise derivative. Using the model is the component that affects the most the performance, highlighting the importance of our derived estimator.

5.4 MODEL PREDICTIVE CONTROL

One of the key benefits of methods that combine model-based reinforcement learning and actor-critic methods is that the optimization procedure results in a stochastic policy, a dynamics model and a Q-function. Hence, we have all the components for, at test time, refine the action selection by the

Under review as a conference paper at ICLR 2020

means of model predictive control (MPC). Here, we investigate the improvement in performance of planning at test time. Specifically, we use the cross-entropy method with our stochastic policy as our initial distributions. The results, shown in Table 2, show benefits in online planning in complex domains; however, its improvement gains are more timid in easier domains, showing that the learned policy has already interiorized the optimal behaviour.

Ant Env Half Cheetah Env Hopper Env Walker2d Env

MAAC+MPC 3.97e3 1.48e3 1.09e4 94.5 2.8e3 11 1.76e3 78 MAAC 3.06e3 1.45e3 1.07e4 253 2.77e3 3.31 1.61e3 404

Table 1: Performance at test time with (maac+mpc) and without (maac) planning of the converged policy using the MAAC objective.

6 CONCLUSION

In this work, we present model-augmented actor-critic, MAAC, a reinforcement learning algorithm that makes use of a learned model by using the pathwise derivative across future timesteps. We prevent instabilities arisen from backpropagation through time by the means of a terminal value function. The objective is theoretically analyzed in terms of the model and value error, and we derive a policy improvement expression with respect to those terms. Our algorithm that builds on top of MAAC is able to achieve superior performance and sample efficiency than state-of-the-art model-based and model-free reinforcement learning algorithms. For future work, it would be enticing to deploy the presented algorithm on a real-robotic agent.

Brandon Amos, Ivan Dario Jimenez Rodriguez, Jacob Sacks, Byron Boots, and J. Zico Kolter. Differentiable mpc for end-to-end planning and control, 2018.

A. G. Barto, R. S. Sutton, and C. W. Anderson. Neuronlike adaptive elements that can solve difficult learning control problems. IEEE Transactions on Systems, Man, and Cybernetics, SMC-13(5): 834 846, Sep. 1983. doi: 10.1109/TSMC.1983.6313077.

Jacob Buckman, Danijar Hafner, George Tucker, Eugene Brevdo, and Honglak Lee. Sample-efficient reinforcement learning with stochastic ensemble value expansion. Co RR, abs/1807.01675, 2018. URL http://arxiv.org/abs/1807.01675.

Kurtland Chua, Roberto Calandra, Rowan Mc Allister, and Sergey Levine. Deep reinforcement learning in a handful of trials using probabilistic dynamics models. ar Xiv preprint ar Xiv:1805.12114, 2018.

Ignasi Clavera, Jonas Rothfuss, John Schulman, Yasuhiro Fujita, Tamim Asfour, and Pieter Abbeel. Model-based reinforcement learning via meta-policy optimization. Co RR, abs/1809.05214, 2018.

Marc Deisenroth and Carl E Rasmussen. Pilco: A model-based and data-efficient approach to policy search. In Proceedings of the 28th International Conference on machine learning (ICML-11), pp. 465 472, 2011.

Yilun Du and Karthik Narasimhan. Task-agnostic dynamics priors for deep reinforcement learning. Co RR, abs/1905.04819, 2019. URL http://arxiv.org/abs/1905.04819.

Vladimir Feinberg, Alvin Wan, Ion Stoica, Michael I Jordan, Joseph E Gonzalez, and Sergey Levine. Model-based value estimation for efficient model-free reinforcement learning. ar Xiv preprint ar Xiv:1803.00101, 2018.

Scott Fujimoto, Herke van Hoof, and David Meger. Addressing function approximation error in actor-critic methods. ar Xiv preprint ar Xiv:1802.09477, 2018.

Under review as a conference paper at ICLR 2020

G.R. Grimmett and D.R. Stirzaker. Probability and random processes, volume 80. Oxford university press, 2001. URL http://scholar.google.com/scholar.bib?q=info: xz St ZXK20Nk J:scholar.google.com/&output=citation&hl=en&as_sdt=0, 5&ct=citation&cd=0.

Tuomas Haarnoja, Aurick Zhou, Pieter Abbeel, and Sergey Levine. Soft actor-critic: Off-policy maximum entropy deep reinforcement learning with a stochastic actor. ar Xiv preprint ar Xiv:1801.01290, 2018a.

Tuomas Haarnoja, Aurick Zhou, Kristian Hartikainen, George Tucker, Sehoon Ha, Jie Tan, Vikash Kumar, Henry Zhu, Abhishek Gupta, Pieter Abbeel, and Sergey Levine. Soft actor-critic algorithms and applications. Co RR, abs/1812.05905, 2018b.

Nicolas Heess, Gregory Wayne, David Silver, Timothy Lillicrap, Tom Erez, and Yuval Tassa. Learning continuous control policies by stochastic value gradients. In Advances in Neural Information Processing Systems, pp. 2944 2952, 2015.

Zhang-Wei Hong, Joni Pajarinen, and Jan Peters. Model-based lookahead reinforcement learning. Ar Xiv, abs/1908.06012, 2019.

Michael Janner, Justin Fu, Marvin Zhang, and Sergey Levine. When to trust your model: Model-based policy optimization. Co RR, abs/1906.08253, 2019.

Sham Kakade and John Langford. Approximately optimal approximate reinforcement learning. In IN PROC. 19TH INTERNATIONAL CONFERENCE ON MACHINE LEARNING, pp. 267 274, 2002.

Gabriel Kalweit and Joschka Boedecker. Uncertainty-driven imagination for continuous deep reinforcement learning. In Sergey Levine, Vincent Vanhoucke, and Ken Goldberg (eds.), Proceedings of the 1st Annual Conference on Robot Learning, volume 78 of Proceedings of Machine Learning Research, pp. 195 206. PMLR, 13 15 Nov 2017.

Péter Karkus, Xiao Ma, David Hsu, Leslie Pack Kaelbling, Wee Sun Lee, and Tomás Lozano-Pérez. Differentiable algorithm networks for composable robot learning. Co RR, abs/1905.11602, 2019. URL http://arxiv.org/abs/1905.11602.

Diederik P Kingma and Max Welling. Auto-encoding variational bayes. ar Xiv preprint ar Xiv:1312.6114, 2013.

Thanard Kurutach, Ignasi Clavera, Yan Duan, Aviv Tamar, and Pieter Abbeel. Model-ensemble trust-region policy optimization. ar Xiv preprint ar Xiv:1802.10592, 2018.

Sergey Levine and Pieter Abbeel. Learning neural network policies with guided policy search under unknown dynamics. In Advances in Neural Information Processing Systems, pp. 1071 1079, 2014.

Timothy P Lillicrap, Jonathan J Hunt, Alexander Pritzel, Nicolas Heess, Tom Erez, Yuval Tassa, David Silver, and Daan Wierstra. Continuous control with deep reinforcement learning. ar Xiv preprint ar Xiv:1509.02971, 2015.

Long-Ji Lin. Self-improving reactive agents based on reinforcement learning, planning and teaching. Machine Learning, 8(3):293 321, May 1992. ISSN 1573-0565. doi: 10.1007/BF00992699.

Kendall Lowrey, Aravind Rajeswaran, Sham M. Kakade, Emanuel Todorov, and Igor Mordatch. Plan online, learn offline: Efficient learning and exploration via model-based control. Co RR, abs/1811.01848, 2018.

Yuping Luo, Huazhe Xu, Yuanzhi Li, Yuandong Tian, Trevor Darrell, and Tengyu Ma. Algorithmic framework for model-based deep reinforcement learning with theoretical guarantees. ICLR, 2019.

Volodymyr Mnih, Adria Puigdomenech Badia, Mehdi Mirza, Alex Graves, Timothy Lillicrap, Tim Harley, David Silver, and Koray Kavukcuoglu. Asynchronous methods for deep reinforcement learning. In International Conference on Machine Learning, pp. 1928 1937, 2016.

Shakir Mohamed, Mihaela Rosca, Michael Figurnov, and Andriy Mnih. Monte carlo gradient estimation in machine learning, 2019.

Under review as a conference paper at ICLR 2020

Anusha Nagabandi, Gregory Kahn, Ronald S Fearing, and Sergey Levine. Neural network dynamics for model-based deep reinforcement learning with model-free fine-tuning. ar Xiv preprint ar Xiv:1708.02596, 2017.

Masashi Okada, Luca Rigazio, and Takenobu Aoshima. Path integral networks: End-to-end differentiable optimal control, 2017.

Marcus Pereira, David D. Fan, Gabriel Nakajima An, and Evangelos A. Theodorou. Mpc-inspired neural network policies for sequential decision making. Co RR, abs/1802.05803, 2018. URL http://arxiv.org/abs/1802.05803.

J. Peters and S. Schaal. Policy gradient methods for robotics. In 2006 IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 2219 2225, Oct 2006. doi: 10.1109/IROS.2006. 282564.

John Schulman, Nicolas Heess, Theophane Weber, and Pieter Abbeel. Gradient estimation using stochastic computation graphs. Co RR, abs/1506.05254, 2015a. URL http://arxiv.org/ abs/1506.05254.

John Schulman, Sergey Levine, Pieter Abbeel, Michael Jordan, and Philipp Moritz. Trust region policy optimization. In Proceedings of the 32nd International Conference on Machine Learning (ICML-15), pp. 1889 1897, 2015b.

Aravind Srinivas, Allan Jabri, Pieter Abbeel, Sergey Levine, and Chelsea Finn. Universal planning networks. ar Xiv preprint ar Xiv:1804.00645, 2018.

Richard S Sutton. Planning by incremental dynamic programming. In Machine Learning Proceedings 1991, pp. 353 357. Elsevier, 1991.

Richard S. Sutton and Andrew G. Barto. Introduction to Reinforcement Learning. MIT Press, Cambridge, MA, USA, 1st edition, 1998. ISBN 0262193981.

Aviv Tamar, Sergey Levine, and Pieter Abbeel. Value iteration networks. Co RR, abs/1602.02867, 2016. URL http://arxiv.org/abs/1602.02867.

Emanuel Todorov, Tom Erez, and Yuval Tassa. Mujoco: A physics engine for model-based control. In Intelligent Robots and Systems (IROS), 2012 IEEE/RSJ International Conference on, pp. 5026 5033. IEEE, 2012.

Tingwu Wang and Jimmy Ba. Exploring model-based planning with policy networks. Co RR, abs/1906.08649, 2019. URL http://arxiv.org/abs/1906.08649.

Tingwu Wang, Xuchan Bao, Ignasi Clavera, Jerrick Hoang, Yeming Wen, Eric Langlois, Shunshi Zhang, Guodong Zhang, Pieter Abbeel, and Jimmy Ba. Benchmarking model-based reinforcement learning. Co RR, abs/1907.02057, 2019.

Under review as a conference paper at ICLR 2020

Here we prove the lemmas and theorems stated in the manuscript.

A.1 PROOF OF LEMMA 4.1

Let Jπ(θ) and ˆJπ(ˆθ) be the expected return of the policy πθ under our objective and under the RL objective, respectively. Then, we can write the MSE of the gradient as

E[ θJπ(θ) θ ˆJπ(θ) 2] = E[ θ(M ˆ M) + | θγH(Q ˆQ) 2]

E[ θ(M ˆ M) 2] + E[ θγH(Q ˆQ) 2]

whereby, M = PH t=0 γtr(st) and ˆ M = PH t=0 γtr(ˆst).

We will denote as the gradient w.r.t the inputs of network, xt = (st, at) and ˆxt = (ˆst, ˆat); where ˆat πθ(ˆst). Notice that since f ˆf and π are Lipschitz and their gradient is Lipschitz as well, we have that θˆxt Kt, where K depends on the Lipschitz constants of the model and the policy. Without loss of generality, we assume that K is larger than 1. Now, we can bound the error on the Q as

θ(Q ˆQ) 2 = Q θx H ˆQ θˆx H 2

= ( Q ˆQ) θx H ˆQ( θˆx H θx H) 2

Q ˆQ 2 θx H 2 + ˆQ 2 θˆx H θx H 2 ϵQ θx H 2 + LQ θˆx H θx H 2 ϵQKH + LQ θˆx H θx H 2

Now, we will bound the term θˆst+1 θst+1 2:

θˆst+1 θst+1 2 = s ˆf θˆst + a ˆf θˆat sf θst af θat 2

s ˆf θˆst sf θst 2 + a ˆf θˆat af θat 2 ϵf θˆst 2 + Lf θˆst θst 2 + Lf θˆat θat + ϵf θˆat 2 ϵf θˆst 2 + (Lf + Lf Lπ) θˆst θst 2 + ϵf θˆat = ϵf θˆxt 2 + (Lf + Lf Lπ) θˆst θst 2

Hence, applying this recursion we obtain that

θˆxt+1 θxt+1 2 ϵf

k=0 (Lf + Lf Lπ)t k θˆxk 2 ϵf Lt+1 1

where L = Lf + Lf Lπ. Then, the error in the gradient in the previous term is bounded by

θ(Q ˆQ) 2 ϵQKH + LQϵf LH 1

In order to bound the model term we need first to bound the rewards since

t=0 γt θ(r(st) r(ˆst)) 2

Similar to the previous bounds, we can bound now each reward term by

θ(r(st) r(ˆst)) 2 ϵf Lr Lt+1 1

With this result we can bound the total error in models

t=0 γtϵf Lr Lt+1 1

L 1 Kt = Lϵf (L 1)

γKL 1 (γK)H 1

Under review as a conference paper at ICLR 2020

Then, the gradient error has the form

E[ θJπ(θ) θ ˆJπ(θ) 2] Lϵf (L 1)

γKL 1 (γK)H 1

+ ϵQ(γK)H + LQϵf LH 1

= ϵfc1(H) + ϵQc2(H)

A.2 PROOF OF LEMMA 4.2

The total variation distance can be bounded by the KL-divergence using the Pinsker s inequality

DT V (πθ πˆθ)

DKL(πθ πˆθ)

Then if we assume third order smoothness on our policy, by the Fisher information metric theorem then

DKL(πθ πˆθ) = c θ ˆθ 2 2 + O( θ ˆθ 3 2)

Given that θ ˆθ 2 = α θJπ θ ˆJπ 2, for a small enough step the following inequality holds

DKL(πθ πˆθ) α2 c(ϵfc1(H) + ϵQc2(H))2 =

Combining this bound with the Pinsker inequality

DT V (πθ πˆθ) α

c 2(ϵfc1(H) + ϵQc2(H)) = αc3(ϵfc1(H) + ϵQc2(H))

A.3 PROOF OF THEOREM 4.1

Given the bound on the total variation distance, we can now make use of the monotonic improvement theorem to establish an improvement bound in terms of the gradient error. Let Jπ(θ) and Jπ(ˆθ) be the expected return of the policy πθ and πˆθ under the true dynamics. Let ρ and ˆρ be the discounted state marginal for the policy πθ and πˆθ, respectively

|Jπ(θ) Jπ(ˆθ)| =| X

s,a ρ(s)πθr(s, a) ˆρ(s)πˆθr(s, a)|

s,a ρ(s)πθ(a|s)r(s, a) ˆρ(s)πˆθ(a|s)r(s, a)|

s,a ρ(s)πθ(a|s) ˆρ(s)πˆθ(a|s)|

a |πθ(a|s) πˆθ(a|s)|

1 γ max s DT V (πθ πˆθ)

Then, combining the results from Lemma 4.2 we obtain the desired bound.

A.4 ABLATIONS

In order to show the significance of each component of MAAC, we conducted more ablation studies. The results are shown in Figure 5. Here, we analyze the effect of training the Q-function with data coming from just the real environment, not learning a maximum entropy policy, and increasing the batch size instead of increasing the amount of samples to estimate the value function.

Under review as a conference paper at ICLR 2020

Figure 5: We further test the significance of some components of our method: not use the dynamics to generate data, and only use real data sampled from environments to train policy and Q-functions (real_data), remove entropy from optimization objects (no_entropy), and using a single sample estimate of the pathwise derivative but increase the batch size accordingly (5x batch size). Considering entropy and using dynamic models to augment data set are both very important.

A.5 EXECUTION TIME COMPARISON

Iteration (s) Training Model (s) Optimization (s) MBPO Iteration (s)

Half Cheetah Env 1312 486 738 708 Hopper Env 845 209 517 723

Table 2: This table shows the time that different parts of MAAC need to train for one iteration after 6000 time steps, averaged across 4 seeds. We also add the time needed for MBPO for one iteration here for comparison.