Imitation-Projected Programmatic Reinforcement Learning

  • 2019-11-04 12:32:34
  • Abhinav Verma, Hoang M. Le, Yisong Yue, Swarat Chaudhuri
  • 0

Abstract

We study the problem of programmatic reinforcement learning, in whichpolicies are represented as short programs in a symbolic language. Programmaticpolicies can be more interpretable, generalizable, and amenable to formalverification than neural policies; however, designing rigorous learningapproaches for such policies remains a challenge. Our approach to thischallenge -- a meta-algorithm called PROPEL -- is based on three insights.First, we view our learning task as optimization in policy space, modulo theconstraint that the desired policy has a programmatic representation, and solvethis optimization problem using a form of mirror descent that takes a gradientstep into the unconstrained policy space and then projects back onto theconstrained space. Second, we view the unconstrained policy space as mixingneural and programmatic representations, which enables employingstate-of-the-art deep policy gradient approaches. Third, we cast the projectionstep as program synthesis via imitation learning, and exploit contemporarycombinatorial methods for this task. We present theoretical convergence resultsfor PROPEL and empirically evaluate the approach in three continuous controldomains. The experiments show that PROPEL can significantly outperformstate-of-the-art approaches for learning programmatic policies.

 

Quick Read (beta)

Imitation-Projected Programmatic Reinforcement Learning

Abhinav Verma
Rice University
[email protected]
&Hoang M. Le*
Caltech
[email protected]
&Yisong Yue
Caltech
[email protected]
&Swarat Chaudhuri
Rice University
[email protected]
Equal contribution

Appendix: Imitation-Projected Programmatic Reinforcement Learning

Abhinav Verma
Rice University
[email protected]
&Hoang M. Le*
Caltech
[email protected]
&Yisong Yue
Caltech
[email protected]
&Swarat Chaudhuri
Rice University
[email protected]
Equal contribution
Abstract

We study the problem of programmatic reinforcement learning, in which policies are represented as short programs in a symbolic language. Programmatic policies can be more interpretable, generalizable, and amenable to formal verification than neural policies; however, designing rigorous learning approaches for such policies remains a challenge. Our approach to this challenge — a meta-algorithm called Propel— is based on three insights. First, we view our learning task as optimization in policy space, modulo the constraint that the desired policy has a programmatic representation, and solve this optimization problem using a form of mirror descent that takes a gradient step into the unconstrained policy space and then projects back onto the constrained space. Second, we view the unconstrained policy space as mixing neural and programmatic representations, which enables employing state-of-the-art deep policy gradient approaches. Third, we cast the projection step as program synthesis via imitation learning, and exploit contemporary combinatorial methods for this task. We present theoretical convergence results for Propel and empirically evaluate the approach in three continuous control domains. The experiments show that Propel can significantly outperform state-of-the-art approaches for learning programmatic policies.

 

Imitation-Projected Programmatic Reinforcement Learning


  Abhinav Vermathanks: Equal contribution Rice University [email protected] Hoang M. Le* Caltech [email protected] Yisong Yue Caltech [email protected] Swarat Chaudhuri Rice University [email protected]

\@float

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

1 Introduction

A growing body of work [58, 8, 60] investigates reinforcement learning (RL) approaches that represent policies as programs in a symbolic language, e.g., a domain-specific language for composing control modules such as PID controllers [5]. Short programmatic policies offer many advantages over neural policies discovered through deep RL, including greater interpretability, better generalization to unseen environments, and greater amenability to formal verification. These benefits motivate developing effective approaches for learning such programmatic policies.

However, programmatic reinforcement learning (PRl) remains a challenging problem, owing to the highly structured nature of the policy space. Recent state-of-the-art approaches employ program synthesis methods to imitate or distill a pre-trained neural policy into short programs [58, 8]. However, such a distillation process can yield a highly suboptimal programmatic policy — i.e., a large distillation gap — and the issue of direct policy search for programmatic policies also remains open.

In this paper, we develop Propel (Imitation-Projected Programmatic Reinforcement Learning), a new learning meta-algorithm for PRl, as a response to this challenge. The design of Propel is based on three insights that enables integrating and building upon state-of-the-art approaches for policy gradients and program synthesis. First, we view programmatic policy learning as a constrained policy optimization problem, in which the desired policies are constrained to be those that have a programmatic representation. This insight motivates utilizing constrained mirror descent approaches, which take a gradient step into the unconstrained policy space and then project back onto the constrained space. Second, by allowing the unconstrained policy space to have a mix of neural and programmatic representations, we can employ well-developed deep policy gradient approaches [55, 36, 47, 48, 19] to compute the unconstrained gradient step. Third, we define the projection operator using program synthesis via imitation learning [58, 8], in order to recover a programmatic policy from the unconstrained policy space. Our contributions can be summarized as:

  • We present Propel, a novel meta-algorithm that is based on mirror descent, program synthesis, and imitation learning, for PRl.

  • On the theoretical side, we show how to cast Propel as a form of constrained mirror descent. We provide a thorough theoretical analysis characterizing the impact of approximate gradients and projections. Further, we prove results that provide expected regret bounds and finite-sample guarantees under reasonable assumptions.

  • On the practical side, we provide a concrete instantiation of Propel and evaluate it in three continuous control domains, including the challenging car-racing domain Torcs [59]. The experiments show significant improvements over state-of-the-art approaches for learning programmatic policies.

2 Problem Statement

The problem of programmatic reinforcement learning (PRl) consists of a Markov Decision Process (Mdp) and a programmatic policy class Π. The definition of =(𝒮,𝒜,P,c,p0,γ) is standard [54], with 𝒮 being the state space, 𝒜 the action space, P(s|s,a) the probability density function of transitioning from a state-action pair to a new state, c(s,a) the state-action cost function, p0(s) a distribution over starting states, and γ(0,1) the discount factor. A policy π:𝒮𝒜 (stochastically) maps states to actions. We focus on continuous control problems, so 𝒮 and 𝒜 are assumed to be continuous spaces. The goal is to find a programmatic policy π*Π such that:

π*=argminπΠJ(π),where: J(π)=𝐄[i=0γic(si,aiπ(si))], (1)

with the expectation taken over the initial state distribution s0p0, the policy decisions, and the transition dynamics P. One can also use rewards, in which case (1) becomes a maximization problem.

Programmatic Policy Class. A programmatic policy class Π consists of policies that can be represented parsimoniously by a (domain-specific) programming language. Recent work [58, 8, 60] indicates that such policies can be easier to interpret and formally verify than neural policies, and can also be more robust to changes in the environment.

In this paper, we consider two concrete classes of programmatic policies. The first, a simplification of the class considered in Verma et al. [58], is defined by the modular, high-level language in Figure 1. This language assumes a library of parameterized functions θ representing standard controllers, for instance Proportional-Integral-Derivative (PID) [6] or bang-bang controllers [11]. Programs in the language take states s as inputs and produce actions a as output, and can invoke fully instantiated library controllers along with predefined arithmetic, boolean and relational operators. The second, “lower-level" class, from Bastani et al. [8], consists of decision trees that map states to actions.

π(s) ::= a𝑂𝑝(π1(s),,πk(s))𝐢𝐟b𝐭𝐡𝐞𝐧π1(s)𝐞𝐥𝐬𝐞π2(s)θ(π1(s),,πk(s))
b ::= ϕ(s)𝐵𝑂𝑝(b1,,bk)
Figure 1: A high-level syntax for programmatic policies, inspired by [58]. A policy π(s) takes a state s as input and produces an action a as output. b represents boolean expressions; ϕ is a boolean-valued operator on states; 𝑂𝑝 is an operator that combines multiple policies into one policy; 𝐵𝑂𝑝 is a standard boolean operator; and θ is a “library function" parameterized by θ.
𝐢𝐟(s[𝚃𝚛𝚊𝚌𝚔𝙿𝚘𝚜]<0.011𝐚𝐧𝐝s[𝚃𝚛𝚊𝚌𝚔𝙿𝚘𝚜]>-0.011)𝐭𝐡𝐞𝐧𝙿𝙸𝙳𝚁𝙿𝙼,0.45,3.54,0.03,53.39(s)𝐞𝐥𝐬𝐞𝙿𝙸𝙳𝚁𝙿𝙼,0.39,3.54,0.03,53.39(s)
Figure 2: A programmatic policy for acceleration in Torcs [59], automatically discovered by Propel. s[TrackPos] represents the most recent reading from sensor TrackPos.

Example. Consider the problem of learning a programmatic policy, in the language of Figure 1, that controls a car’s accelerator in the Torcs car-racing environment [59]. Figure 2 shows a program in our language for this task. The program invokes PID controllers 𝙿𝙸𝙳j,θP,θI,θD, where j identifies the sensor (out of 29, in our experiments) that provides inputs to the controller, and θP, θI, and θD are respectively the real-valued coefficients of the proportional, integral, and derivative terms in the controller. We note that the program only uses the sensors 𝚃𝚛𝚊𝚌𝚔𝙿𝚘𝚜 and 𝚁𝙿𝙼. While 𝚃𝚛𝚊𝚌𝚔𝙿𝚘𝚜 (for the position of the car relative to the track axis) is used to decide which controller to use, only the 𝚁𝙿𝙼 sensor is needed to calculate the acceleration.

Learning Challenges. Learning programmatic policies in the continuous RL setting is challenging, as the best performing methods utilize policy gradient approaches [55, 36, 47, 48, 19], but policy gradients are hard to compute in programmatic representations. In many cases, Π may not even be differentiable. For our approach, we only assume access to program synthesis methods that can select a programmatic policy πΠ that minimizes imitation disagreement with demonstrations provided by a teaching oracle. Because imitation learning tends to be easier than general RL in long-horizon tasks [53], the task of imitating a neural policy with a program is, intuitively, significantly simpler than the full programmatic RL problem. This intuition is corroborated by past work on programmatic RL [58], which shows that direct search over programs often fails to meet basic performance objectives.

{algorithm}

[t] Imitation-Projected Programmatic Reinforcement Learning (Propel) \[email protected]@algorithmic[1] \STATEInput: Programmatic & Neural Policy Classes: Π & . \STATEInput: Either initial π0 or initial f0 \STATEDefine joint policy class: Π         //hπ+f defined as h(s)=π(s)+f(s) \IFgiven initial f0 \STATEπ0Project(f0)          //program synthesis via imitation learning \ENDIF\FORt=1,,T \STATEhtUpdate(πt-1,η)          //policy gradient in neural policy space with learning rate η \STATEπtProjectΠ(ht)          //program synthesis via imitation learning \ENDFOR\STATEReturn: Policy πT

3 Learning Algorithm

To develop our approach, we take the viewpoint of (1) being a constrained optimization problem, where Π resides within a larger space of policies . In particular, we will represent Π using a mixing of programmatic policies Π and neural polices . Any mixed policy hπ+f can be invoked as h(s)=π(s)+f(s). In general, we assume that is a good approximation of Π (i.e., for each πΠ there is some f that approximates it well), which we formalize in Section 4.

We can now frame our constrained learning problem as minimizing (1) over Π, that alternate between taking a gradient step in the general space and projecting back down onto Π. This “lift-and-project” perspective motivates viewing our problem via the lens of mirror descent [40]. In standard mirror descent, the unconstrained gradient step can be written as hhprev-ηJ(hprev) for step size η, and the projection can be written as πargminπΠD(π,h) for divergence measure D.

Our approach, Imitation-Projected Programmatic Reinforcement Learning (Propel), is outlined in Algorithm 2 (also see Figure 3). Propel is a meta-algorithm that requires instantiating two subroutines, Update and Project, which correspond to the standard update and projection steps, respectively. Propel can be viewed as a form of functional mirror descent with some notable deviations from vanilla mirror descent.

Figure 3: Depicting the Propel meta-algorithm.

Update. Since policy gradient methods are well-developed for neural policy classes (e.g., [36, 47, 48, 30, 24, 19]) and non-existent for programmatic policy classes Π, Propel is designed to leverage policy gradients in and avoid policy gradients in Π. Algorithm 3 shows one instantiation of Update. Note that standard mirror descent takes unconstrained gradient steps in rather than , and we discuss this discrepancy between Update and Update in Section 4.

ProjectΠ. Projecting onto Π can be implemented using program synthesis via imitation learning, i.e., by synthesizing a πΠ to best imitate demonstrations provided by a teaching oracle h. Recent work [58, 8, 60] has given practical heuristics for this task for various programmatic policy classes. Algorithm 3 shows one instantiation of ProjectΠ (based on DAgger [46]). One complication that arises is that finite-sample runs of such imitation learning approaches only return approximate solutions and so the projection is not exact. We characterize the impact of approximate projections in Section 4.

Practical Considerations. In practice, we often employ multiple gradient steps before taking a projection step (as also described in Algorithm 3), because the step size of individual (stochastic) gradient updates can be quite small. Another issue that arises in virtually all policy gradient approaches is that the gradient estimates can have very high variance [55, 33, 30]. We utilize low-variance policy gradient updates by using the reference π as a proximal regularizer in function space [19].

For the projection step (Algorithm 3), in practice we often retain all previous roll-outs τ from all previous projection steps. It is straightforward to query the current oracle h to provide demonstrations on the states sτ from previous roll-outs, which can lead to substantial savings in sample complexity with regards to executing roll-outs on the environment, while not harming convergence.

{algorithm}

[t] Update: neural policy gradient for mixed policies \[email protected]@algorithmic[1] \STATEInput: Neural Policy Class .   Input: Reference programmatic policy: π \STATEInput: Step size: η.   Input: Regularization parameter: λ \STATEInitialize neural policy: f0         //any standard randomized initialization \FORj=1,,m \STATEfjfj-1-ηλJ(π+λfj-1)         //using DDPG [36], TRPO [47], etc., holding π fixed \ENDFOR\STATEReturn: hπ+λfm

{algorithm}

[t] ProjectΠ: program synthesis via imitation learning \[email protected]@algorithmic[1] \STATEInput: Programmatic Policy Class: Π.   Input: Oracle policy: h \STATERoll-out h on environment, get trajectory: τ0=(s0,h(s0),s1,h(s1),) \STATECreate supervised demonstration set: Γ0={(s,h(s))} from τ0 \STATEDerive π0 from Γ0 via program synthesis          //e.g., using methods in [58, 8] \FORk=1,,M \STATERoll-out πk-1, creating trajectory: τk \STATECollect demonstration data: Γ={(s,h(s))|sτk} \STATEΓkΓΓk-1          //DAgger-style imitation learning [46] \STATEDerive πk from Γk via program synthesis          //e.g., using methods in [58, 8] \ENDFOR\STATEReturn: πM

4 Theoretical Analysis

We start by viewing Propel through the lens of online learning in function space, independent of the specific parametric representation. This start point yields a convergence analysis of Alg. 2 in Section 4.1 under generic approximation errors. We then analyze the issues of policy class representation in Sections 4.2 and 4.3, and connect Algorithms 3 and 3 with the overall performance, under some simplifying conditions. In particular, Section 4.3 characterizes the update error in a possibly non-differentiable setting; to our knowledge, this is the first such analysis of its kind for reinforcement learning.

Preliminaries. We consider Π and to be subspaces of an ambient policy space 𝒰, which is a vector space equipped with inner product ,, induced norm u=u,u, dual norm v*=sup{v,u|u1}, and standard scaling & addition: (au+bv)(s)=au(s)+bv(s) for a,b and u,v𝒰. The cost functional of a policy u is J(u)=𝒮c(s,u(s))𝑑μu(s), where μu is the distribution of states induced by u. The joint policy class is =Π, by ={π+f|πΠ,f}.11 1 The operator is not a direct sum, since Π and are not orthogonal. Note that is a subspace of 𝒰, and inherits its vector space properties. Without affecting the analysis, we simply equate 𝒰 for the remainder of the paper.

We assume that J is convex in , which implies that subgradient J(h) exists (with respect to ) [9]. Where J is differentiable, we utilize the notion of a Fréchet gradient. Recall that a bounded linear operator : is called a Fréchet functional gradient of J at h if limg0J(h+g)-J(h)-J(h),gg=0. By default, (or for emphasis) denotes the gradient with respect to , whereas defines the gradient in the restricted subspace .

4.1 Propel as (Approximate) Functional Mirror Descent

For our analysis, Propel can be viewed as approximating mirror descent in (infinite-dimensional) function space over a convex set Π.22 2 Π can be convexified by considering randomized policies, as stochastic combinations of πΠ (cf. [35]). Similar to the finite-dimensional setting [40], we choose a strongly convex and smooth functional regularizer R to be the mirror map. From the approximate mirror descent perspective, for each iteration t:

  1. 1.

    Obtain a noisy gradient estimate: ^t-1J(πt-1)

  2. 2.

    Update(π) in space: R(ht)=R(πt-1)-η^t-1 (Note UpdateHUpdateF)

  3. 3.

    Obtain approximate projection: πt=ProjectΠR(ht)argminπΠDR(π,ht)

DR(u,v)=R(u)-R(v)-R(u),u-v is a Bregman divergence. Taking R(h)=12h2 will recover projected functional gradient descent in L2-space. Here Update becomes ht=πt-1-η^J(πt-1), and Project solves for argminπΠπ-ht2. While we mainly focus on this choice of R in our experiments, note that other selections of R lead to different Update and Project operators (e.g., minimizing KL divergence if R is negative entropy).

The functional mirror descent scheme above may encounter two additional sources of error compared to standard mirror descent [40]. First, in the stochastic setting (also called bandit feedback [28]), the gradient estimate ^t may be biased, in addition to having high variance. One potential source of bias is the gap between Update and Update. Second, the Project step may be inexact. We start by analyzing the behavior of Propel under generic bias, variance, and projection errors, before discussing the implications of approximating Update and ProjectΠ by Algs. 3 & 3, respectively. Let the bias be bounded by β, i.e., 𝔼[^t|πt]-J(πt)*β almost surely. Similarly let the variance of the gradient estimate be bounded by σ2, and the projection error norm πt-πt*ϵ. We state the expected regret bound below; more details and a proof appear in Appendix A.2.

Theorem 4.1 (Expected regret bound under gradient estimation and projection errors).

Let π1,,πT be a sequence of programmatic policies returned by Algorithm 2, and π* be the optimal programmatic policy. Choosing learning rate η=1σ2(1T+ϵ), we have the expected regret over T iterations:

𝔼[1Tt=1TJ(πt)]-J(π*)=O(σ1T+ϵ+β). (2)

The result shows that error ϵ from Project and the bias β do not accumulate and simply contribute an additive term on the expected regret.33 3 Other mirror descent-style analyses, such as in [52], lead to accumulation of errors over the rounds of learning T. One key difference is that we are leveraging the assumption of convexity of J in the (infinite-dimensional) function space representation. The effect of variance of gradient estimate decreases at a 1/T rate. Note that this regret bound is agnostic to the specific Update and Project operations, and can be applied more generically beyond the specific algorithmic choices used in our paper.

4.2 Finite-Sample Analysis under Vanilla Policy Gradient Update and DAgger Projection

Next, we show how certain instantiations of Update and Project affect the magnitude of errors and influence end-to-end learning performance from finite samples, under some simplifying assumptions on the Update step. For this analysis, we simplify Alg. 3 into the case UpdateUpdate. In particular, we assume programmatic policies in Π to be parameterized by a vector θk, and π is differentiable in θ (e.g., we can view Π where is parameterized in k). We further assume the trajectory roll-out is performed in an exploratory manner, where action is taken uniformly random over finite set of A actions, thus enabling the bound on the bias of gradient estimates via Bernstein’s inequality. The Project step is consistent with Alg. 3, i.e., using DAgger [45] under convex imitation loss, such as 2 loss. We have the following high-probability guarantee:

Theorem 4.2 (Finite-sample guarantee).

At each iteration, we perform vanilla policy gradient estimate of π (over H) using m trajectories and, use DAgger algorithm to collect M roll-outs for the imitation learning projection. Setting the learning rate η=1σ2(1T+HM+log(T/δ)M), after T rounds of the algorithm, we have that:

1Tt=1TJ(πt)-J(π*)O(σ1T+HM+log(T/δ)M)+O(σlog(Tk/δ)m+AHlog(Tk/δ)m)

holds with probability at least 1-δ, with H being the task horizon, A the cardinality of action space, σ2 the variance of policy gradient estimates, and k the dimension Π’s parameterization.

The expanded result and proof are included in Appendix A.3. The proof leverages previous analysis from DAgger [46] and the finite sample analysis of vanilla policy gradient algorithm [32]. The finite-sample regret bound scales linearly with the standard deviation σ of the gradient estimate, while the bias, which is the very last component of the RHS, scales linearly with the task horizon H. Note that the standard deviation σ can be exponential in task horizon H in the worst case [32], and so it is important to have practical implementation strategies to reduce the variance of the Update operation. While conducted in a stylized setting, this analysis provides insight in the relative trade-offs of spending effort in obtaining more accurate projections versus more reliable gradient estimates.

4.3 Closing the gap between Update and Update

Our functional mirror descent analysis rests on taking gradients in : Update(π) involves estimating J(π) in the space. On the other hand, Algorithm 3 performs Update(π) only in the neural policy space . In either case, although J(π) may be differentiable in the non-parametric ambient policy space, it may not be possible to obtain a differentiable parametric programmatic representation in Π. In this section, we discuss theoretical motivations to addressing a practical issue: How do we define and approximate the gradient HJ(π) under a parametric representation? To our knowledge, we are the first to consider such a theoretical question for reinforcement learning.

Defining a consistent approximation of HJ(π). The idea in Update(π) (Line 8 of Alg. 2) is to approximate J(π) by J(f), which has a differentiable representation, at some f close to π (under the norm). Under appropriate conditions on , we show that this approximation is valid.

Proposition 4.3.

Assume that (i) J is Fréchet differentiable on H, (ii) J is also differentiable on the restricted subspace F, and (iii) F is dense in H (i.e., the closure \[email protected]Δ\[email protected]\[email protected]\[email protected]@skewchar\[email protected]@a111F=H). Then for any fixed policy πΠ, define a sequence of policies fkF, k=1,2,), that converges to π: limkfk-π=0. We then have limkFJ(fk)-HJ(π)*=0.

Since the Fréchet gradient is unique in the ambient space , k we have J(fk)=J(fk)J(π) as k (by Proposition 4.3). We thus have an asymptotically unbiased approximation of J(π) via differentiable space as: J(π)J(π)limkJ(fk).44 4 We do not assume J(π) to be differentiable when restricting to the policy subspace Π, i.e., ΠJ(π) may not exist under policy parameterization of Π. Connecting to the result from Theorem 4.1, let σ2 be an upper bound on the policy gradient estimates in the neural policy class F, under an asymptotically unbiased approximation of J(π), the expected regret bound becomes 𝔼[1Tt=1TJ(πt)]-J(π*)=O(σ1T+ϵ).

Bias-variance considerations of UpdateF(π) To further theoretically motivate a practical strategy for Update(π) in Algorithm 3, we utilize an equivalent proximal perspective of mirror descent [10], where Update(π) is equivalent to solving for h=argminhηJ(π),h+DR(h,π).

Proposition 4.4 (Minimizing a relaxed objective).

For a fixed programmatic policy π, with sufficiently small constant λ(0,1), we have that

minhηJ(π),h)+DR(h,π)minfJ(π+λf)-J(π)+J(π),π (3)

Thus, a relaxed Update step is obtained by minimizing the RHS of (3), i.e., minimizing J(π+λf) over f. Each gradient descent update step is now f=f-ηλJ(πt+λf), corresponding to Line 5 of Algorithm 3. For fixed π and small λ, this relaxed optimization problem becomes regularized policy optimization over , which is significantly easier. Functional regularization in policy space around a fixed prior controller π has demonstrated significant reduction in the variance of gradient estimate [19], at the expense of some bias. The below expected regret bound summarizes the impact of this increased bias and reduced variance, with details included in Appendix A.5.

Proposition 4.5 (Bias-variance characterization of Update).

Assuming J(h) is L-strongly smooth over H, i.e., HJ(h) is L-Lipschitz continuous, approximating UpdateH by UpdateF per Alg. 3 leads to the expected regret bound: E[1Tt=1TJ(πt)]-J(π*)=O(λσ1T+ϵ+λ2L2).

Compared to the idealized unbiased approximation in Proposition 4.3, the introduced bias here is related to the inherent smoothness property of cost functional J(h) over the joint policy class , i.e., how close J(π+λf) is to its linear under-approximation J(π)+J(π),λf around π.

5 Experiments

We demonstrate the effectiveness of Propel in synthesizing programmatic controllers in three continuous control environments. For brevity and focus, this section primarily focuses on Torcs 55 5 The code for the Torcs experiments can be found at: https://bitbucket.org/averma8053/propel, a challenging race car simulator environment [59]. Empirical results on two additional classic control tasks, Mountain-Car and Pendulum, are provided in Appendix B; those results follow similar trends as the ones described for Torcs below, and further validate the convergence analysis of Propel.

Figure 4: Median lap-time improvements during multiple iterations of PropelProg over 25 random seeds.

Experimental Setup. We evaluate over five distinct tracks in the Torcs simulator. The difficulty of a track can be characterized by three properties; track length, track width, and number of turns. Our suite of tracks provides environments with varying levels of difficulty for the learning algorithm. The performance of a policy in the Torcs simulator is measured by the lap time achieved on the track. To calculate the lap time, the policies are allowed to complete a three-lap race, and we record the best lap time during this race. We perform the experiments with twenty-five random seeds and report the median lap time over these twenty-five trials. Some of the policies crash the car before completing a lap on certain tracks, even after training for 600 episodes. Such crashes are recorded as a lap time of infinity while calculating the median. If the policy crashes for more than half the seeds, this is reported as Cr in Tables 1 & 2. We choose to report the median because taking the crash timing as infinity, or an arbitrarily large constant, heavily skews other common measures such as the mean.

Figure 5: Median number of crashes during training of Ddpg and PropelProg over 25 random seeds.

Baselines. Among recent state-of-the-art approaches to learning programmatic policies are Ndps [58] for high-level language policies, and Viper [8] for learning tree-based policies. Both Ndps and Viper rely on imitating a fixed (pre-trained) neural policy oracle, and can be viewed as degenerate versions of Propel that only run Lines 4-6 in Algorithm 2. We present two Propel analogues to Ndps and Viper: (i) PropelProg: Propel using the high-level language of Figure 1 as the class of programmatic policies, similar to Ndps. (ii) PropelTree: Propel using regression trees, similar to Viper. We also report results for Prior, which is a (sub-optimal) PID controller that is also used as the initial policy in Propel. In addition, to study generalization ability as well as safety behavior during training, we also include Ddpg, a neural policy learned using the Deep Deterministic Policy Gradients [36] algorithm, with 600 episodes of training for each track. In principle, Propel and its analysis can accommodate different policy gradient subroutines. However, in the Torcs domain, other policy gradient algorithms such as PPO and TRPO failed to learn policies that are able to complete the considered tracks. We thus focus on Ddpg as our main policy gradient component.

Table 1: Performance results in Torcs over 25 random seeds. Each entry is formatted as Lap-time / Crash-ratio, reporting median lap time in seconds over all the seeds (lower is better) and ratio of seeds that result in crashes (lower is better). A lap time of Cr indicates the agent crashed and could not complete a lap for more than half the seeds.
G-Track E-Road Aalborg Ruudskogen Alpine-2
Length 3186m 3260m 2588m 3274m 3774m
Prior 312.92 / 0.0 322.59 / 0.0 244.19 / 0.0 340.29 / 0.0 402.89 / 0.0
Ddpg 78.82 / 0.24 89.71 / 0.28 101.06 / 0.40 Cr / 0.68 Cr / 0.92
Ndps 108.25 / 0.24 126.80 / 0.28 163.25 / 0.40 Cr / 0.68 Cr / 0.92
Viper 83.60 / 0.24 87.53 / 0.28 110.57 / 0.40 Cr / 0.68 Cr / 0.92
PropelProg 93.67 / 0.04 119.17 / 0.04 147.28 / 0.12 124.58 / 0.16 256.59 / 0.16
PropelTree 78.33 / 0.04 79.39 / 0.04 109.83 / 0.16 118.80 / 0.24 236.01 / 0.36

Evaluating Performance. Table 1 shows the performance on the considered Torcs tracks. We see that PropelProg and PropelTree consistently outperform the Ndps [58] and Viper [8] baselines, respectively. While Ddpg outperforms Propel on some tracks, its volatility causes it to be unable to learn in some environments, and hence to crash the majority of the time. Figure 4 shows the consistent improvements made over the prior by PropelProg, over the iterations of the Propel algorithm. Appendix B contains similar results achieved on the two classic control tasks, MountainCar and Pendulum. Figure 5 shows that, compared to Ddpg, our approach suffers far fewer crashes while training in Torcs.

Evaluating Generalization. To compare the ability of the controllers to perform on tracks not seen during training, we executed the learned policies on all the other tracks (Table 2). We observe that Ddpg crashes significantly more often than PropelProg. This demonstrates the generalizability of the policies returned by Propel. Generalization results for the PropelTree policy are given in the appendix. In general, PropelTree policies are more generalizable than Ddpg but less than PropelProg. On an absolute level, the generalization ability of Propel still leaves much room for improvement, which is an interesting direction for future work.

Verifiability of Policies. As shown in prior work [8, 58], parsimonious programmatic policies are more amenable to formal verification than neural policies. Unsurprisingly, the policies generated by PropelTree and PropelProg are easier to verify than Ddpg policies. As a concrete example, we verified a smoothness property of the PropelProg policy using the Z3 SMT-solver [21] (more details in Appendix B). The verification terminated in 0.49 seconds.

Initialization. In principle, Propel can be initialized with a random program, or a random policy trained using Ddpg. In practice, the performance of Propel depends to a certain degree on the stability of the policy gradient procedure, which is Ddpg in our experiments. Unfortunately, Ddpg often exhibits high variance across trials and fares poorly in challenging RL domains. Specifically, in our Torcs experiments, Ddpg fails on a number of tracks (similar phenomena have been reported in previous work that experiments on similar continuous control domains [30, 19, 58]). Agents obtained by initializing Propel with neural policies obtained via Ddpg also fail on multiple tracks. Their performance over the five tracks is reported in Appendix B. In contrast, Propel can often finish the challenging tracks when initialized with a very simple hand-crafted programmatic prior.

Table 2: Generalization results in Torcs, where rows are training and columns are testing tracks. Each entry is formatted as PropelProg / DDPG, and the number reported is the median lap time in seconds over all the seeds (lower is better). Cr indicates the agent crashed and could not complete a lap for more than half the seeds.
G-Track E-Road Aalborg Ruudskogen Alpine-2
G-Track - 124 / Cr Cr / Cr Cr / Cr Cr / Cr
E-Road 102 / 92 - Cr / Cr Cr / Cr Cr / Cr
Aalborg 201 / 91 228 / Cr - 217 / Cr Cr / Cr
Ruudskogen 131 / Cr 135 / Cr Cr / Cr - Cr / Cr
Alpine-2 222 / Cr 231 / Cr 184 / Cr Cr / Cr -

6 Related Work

Program Synthesis. Program synthesis is the problem of automatically searching for a program within a language that fits a given specification [29]. Recent approaches to the problem have leveraged symbolic knowledge about program structure [27], satisfiability solvers [50, 31], and meta-learning techniques [39, 41, 22, 7] to generate interesting programs in many domains [3, 42, 4]. In most prior work, the specification is a logical constraint on the input/output behavior of the target program. However, there is also a growing body of work that considers program synthesis modulo optimality objectives [13, 15, 43], often motivated by machine learning tasks [39, 57, 26, 23, 58, 8, 60]. Synthesis of programs that imitates an oracle has been considered in both the logical [31] and the optimization [58, 8, 60] settings. The projection step in Propel builds on this prior work. While our current implementation of this step is entirely symbolic, in principle, the operation can also utilize contemporary techniques for learning policies that guide the synthesis process [39, 7, 49].

Constrained Policy Learning. Constrained policy learning has seen increased interest in recent years, largely due to the desire to impose side guarantees such as stability and safety on the policy’s behavior. Broadly, there are two approaches to imposing constraints: specifying constraints as an additional cost function [1, 35], and explicitly encoding constraints into the policy class [2, 34, 19, 20, 12]. In some cases, these two approaches can be viewed as duals of each other. For instance, recent work that uses control-theoretic policies as a functional regularizer [34, 19] can be viewed from the perspective of both regularization (additional cost) and an explicitly constrained policy class (a specific mix of neural and control-theoretic policies). We build upon this perspective to develop the gradient update step in our approach.

RL using Imitation Learning. There are two ways to utilize imitation learning subroutines within RL. First, one can leverage limited-access or sub-optimal experts to speed up learning [44, 18, 14, 51]. Second, one can learn over two policy classes (or one policy and one model class) to achieve accelerated learning compared to using only one policy class [38, 17, 52, 16]. Our approach has some stylistic similarities to previous efforts [38, 52] that use a richer policy space to search for improvements before re-training the primary policy to imitate the richer policy. One key difference is that our primary policy is programmatic and potentially non-differentiable. A second key difference is that our theoretical framework takes a functional gradient descent perspective — it would be interesting to carefully compare with previous analysis techniques to find a unifying framework.

RL with Mirror Descent. The mirror descent framework has previously used to analyze and design RL algorithms. For example, Thomas et al. [56] and Mahadevan and Liu [37] use composite objective mirror descent, or Comid [25], which allows incorporating adaptive regularizers into gradient updates, thus offering connections to either natural gradient RL [56] or sparsity inducing RL algorithms [37]. Unlike in our work, these prior approaches perform projection into the same native, differentiable representation. Also, the analyses in these papers do not consider errors introduced by hybrid representations and approximate projection operators. However, one can potentially extend our approach with versions of mirror descent, e.g., Comid, that were considered in these efforts.

7 Conclusion and Future Work

We have presented Propel, a meta-algorithm based on mirror descent, program synthesis, and imitation learning, for programmatic reinforcement learning (PRl). We have presented theoretical convergence results for Propel, developing novel analyses to characterize approximate projections and biased gradients within the mirror descent framework. We also validated Propel empirically, and show that it can discover interpretable, verifiable, generalizable, performant policies and significantly outperform the state of the art in PRl.

The central idea of Propel is the use of imitation learning and combinatorial methods in implementing a projection operation for mirror descent, with the goal of optimization in a functional space that lacks gradients. While we have developed Propel in an RL setting, this idea is not restricted to RL or even sequential decision making. Future work will seek to exploit this insight in other machine learning and program synthesis settings.

Acknowledgements. This work was supported in part by United States Air Force Contract # FA8750-19-C-0092, NSF Award # 1645832, NSF Award # CCF-1704883, the Okawa Foundation, Raytheon, PIMCO, and Intel.

References

  • [1] Joshua Achiam, David Held, Aviv Tamar, and Pieter Abbeel. Constrained policy optimization. In Proceedings of the 34th International Conference on Machine Learning-Volume 70, pages 22–31. JMLR. org, 2017.
  • [2] Mohammed Alshiekh, Roderick Bloem, Rüdiger Ehlers, Bettina Könighofer, Scott Niekum, and Ufuk Topcu. Safe reinforcement learning via shielding. In Thirty-Second AAAI Conference on Artificial Intelligence, 2018.
  • [3] Rajeev Alur, Rastislav Bodík, Eric Dallal, Dana Fisman, Pranav Garg, Garvit Juniwal, Hadas Kress-Gazit, P. Madhusudan, Milo M. K. Martin, Mukund Raghothaman, Shambwaditya Saha, Sanjit A. Seshia, Rishabh Singh, Armando Solar-Lezama, Emina Torlak, and Abhishek Udupa. Syntax-guided synthesis. In Dependable Software Systems Engineering, pages 1–25. 2015.
  • [4] Rajeev Alur, Arjun Radhakrishna, and Abhishek Udupa. Scaling enumerative program synthesis via divide and conquer. In Tools and Algorithms for the Construction and Analysis of Systems - 23rd International Conference, TACAS 2017, Held as Part of the European Joint Conferences on Theory and Practice of Software, ETAPS 2017, Uppsala, Sweden, April 22-29, 2017, Proceedings, Part I, pages 319–336, 2017.
  • [5] Kiam Heong Ang, Gregory Chong, and Yun Li. Pid control system analysis, design, and technology. IEEE transactions on control systems technology, 13(4):559–576, 2005.
  • [6] Karl Johan Åström and Tore Hägglund. Automatic tuning of simple regulators with specifications on phase and amplitude margins. Automatica, 20(5):645–651, 1984.
  • [7] Matej Balog, Alexander L. Gaunt, Marc Brockschmidt, Sebastian Nowozin, and Daniel Tarlow. Deepcoder: Learning to write programs. In 5th International Conference on Learning Representations, ICLR 2017, Toulon, France, April 24-26, 2017, Conference Track Proceedings, 2017.
  • [8] Osbert Bastani, Yewen Pu, and Armando Solar-Lezama. Verifiable reinforcement learning via policy extraction. In Advances in Neural Information Processing Systems, pages 2494–2504, 2018.
  • [9] Heinz H Bauschke, Patrick L Combettes, et al. Convex analysis and monotone operator theory in Hilbert spaces, volume 408. Springer, 2011.
  • [10] Amir Beck and Marc Teboulle. Mirror descent and nonlinear projected subgradient methods for convex optimization. Operations Research Letters, 31(3):167–175, 2003.
  • [11] Richard Bellman, Irving Glicksberg, and Oliver Gross. On the “bang-bang” control problem. Quarterly of Applied Mathematics, 14(1):11–18, 1956.
  • [12] Felix Berkenkamp, Matteo Turchetta, Angela Schoellig, and Andreas Krause. Safe model-based reinforcement learning with stability guarantees. In Advances in neural information processing systems, pages 908–918, 2017.
  • [13] Roderick Bloem, Krishnendu Chatterjee, Thomas A. Henzinger, and Barbara Jobstmann. Better quality in synthesis through quantitative objectives. In Computer Aided Verification, 21st International Conference, CAV 2009, Grenoble, France, June 26 - July 2, 2009. Proceedings, pages 140–156, 2009.
  • [14] Kai-Wei Chang, Akshay Krishnamurthy, Alekh Agarwal, Hal Daumé III, and John Langford. Learning to search better than your teacher. In International Conference on Machine Learning (ICML), 2015.
  • [15] Swarat Chaudhuri, Martin Clochard, and Armando Solar-Lezama. Bridging boolean and quantitative synthesis using smoothed proof search. In POPL, pages 207–220, 2014.
  • [16] Ching-An Cheng, Xinyan Yan, Nathan Ratliff, and Byron Boots. Predictor-corrector policy optimization. In International Conference on Machine Learning (ICML), 2019.
  • [17] Ching-An Cheng, Xinyan Yan, Evangelos Theodorou, and Byron Boots. Accelerating imitation learning with predictive models. In International Conference on Artificial Intelligence and Statistics (AISTATS), 2019.
  • [18] Ching-An Cheng, Xinyan Yan, Nolan Wagener, and Byron Boots. Fast policy learning through imitation and reinforcement. In Uncertainty in artificial intelligence, 2019.
  • [19] Richard Cheng, Abhinav Verma, Gabor Orosz, Swarat Chaudhuri, Yisong Yue, and Joel Burdick. Control regularization for reduced variance reinforcement learning. In International Conference on Machine Learning (ICML), 2019.
  • [20] Gal Dalal, Krishnamurthy Dvijotham, Matej Vecerik, Todd Hester, Cosmin Paduraru, and Yuval Tassa. Safe exploration in continuous action spaces. arXiv preprint arXiv:1801.08757, 2018.
  • [21] Leonardo Mendonça de Moura and Nikolaj Bjørner. Z3: An Efficient SMT Solver. In TACAS, pages 337–340, 2008.
  • [22] Jacob Devlin, Jonathan Uesato, Surya Bhupatiraju, Rishabh Singh, Abdel-rahman Mohamed, and Pushmeet Kohli. Robustfill: Neural program learning under noisy i/o. In Proceedings of the 34th International Conference on Machine Learning-Volume 70, pages 990–998. JMLR. org, 2017.
  • [23] Tao Du, Jeevana Priya Inala, Yewen Pu, Andrew Spielberg, Adriana Schulz, Daniela Rus, Armando Solar-Lezama, and Wojciech Matusik. Inversecsg: automatic conversion of 3d models to CSG trees. ACM Trans. Graph., 37(6):213:1–213:16, 2018.
  • [24] Yan Duan, Xi Chen, Rein Houthooft, John Schulman, and Pieter Abbeel. Benchmarking deep reinforcement learning for continuous control. In International Conference on Machine Learning, pages 1329–1338, 2016.
  • [25] John C Duchi, Shai Shalev-Shwartz, Yoram Singer, and Ambuj Tewari. Composite objective mirror descent. In COLT, pages 14–26, 2010.
  • [26] Kevin Ellis, Daniel Ritchie, Armando Solar-Lezama, and Josh Tenenbaum. Learning to infer graphics programs from hand-drawn images. In Advances in Neural Information Processing Systems, pages 6059–6068, 2018.
  • [27] John K. Feser, Swarat Chaudhuri, and Isil Dillig. Synthesizing data structure transformations from input-output examples. In Proceedings of the 36th ACM SIGPLAN Conference on Programming Language Design and Implementation, Portland, OR, USA, June 15-17, 2015, pages 229–239, 2015.
  • [28] Abraham D Flaxman, Adam Tauman Kalai, and H Brendan McMahan. Online convex optimization in the bandit setting: gradient descent without a gradient. In Proceedings of the sixteenth annual ACM-SIAM symposium on Discrete algorithms, pages 385–394. Society for Industrial and Applied Mathematics, 2005.
  • [29] Sumit Gulwani, Oleksandr Polozov, and Rishabh Singh. Program synthesis. Foundations and Trends in Programming Languages, 4(1-2):1–119, 2017.
  • [30] Peter Henderson, Riashat Islam, Philip Bachman, Joelle Pineau, Doina Precup, and David Meger. Deep reinforcement learning that matters. In Thirty-Second AAAI Conference on Artificial Intelligence, 2018.
  • [31] Susmit Jha, Sumit Gulwani, Sanjit A Seshia, and Ashish Tiwari. Oracle-guided component-based program synthesis. In Proceedings of the 32nd ACM/IEEE International Conference on Software Engineering-Volume 1, pages 215–224. ACM, 2010.
  • [32] Sham Machandranath Kakade et al. On the sample complexity of reinforcement learning. PhD thesis, University of London London, England, 2003.
  • [33] Vijay R Konda and John N Tsitsiklis. Actor-critic algorithms. In Advances in neural information processing systems, pages 1008–1014, 2000.
  • [34] Hoang M. Le, Andrew Kang, Yisong Yue, and Peter Carr. Smooth imitation learning for online sequence prediction. In International Conference on Machine Learning (ICML), 2016.
  • [35] Hoang M Le, Cameron Voloshin, and Yisong Yue. Batch policy learning under constraints. In International Conference on Machine Learning (ICML), 2019.
  • [36] 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. arXiv preprint arXiv:1509.02971, 2015.
  • [37] Sridhar Mahadevan and Bo Liu. Sparse q-learning with mirror descent. In Proceedings of the Twenty-Eighth Conference on Uncertainty in Artificial Intelligence, pages 564–573. AUAI Press, 2012.
  • [38] William H Montgomery and Sergey Levine. Guided policy search via approximate mirror descent. In Advances in Neural Information Processing Systems, pages 4008–4016, 2016.
  • [39] Vijayaraghavan Murali, Swarat Chaudhuri, and Chris Jermaine. Neural sketch learning for conditional program generation. In ICLR, 2018.
  • [40] Arkadii Semenovich Nemirovsky and David Borisovich Yudin. Problem complexity and method efficiency in optimization. 1983.
  • [41] Emilio Parisotto, Abdel-rahman Mohamed, Rishabh Singh, Lihong Li, Dengyong Zhou, and Pushmeet Kohli. Neuro-symbolic program synthesis. arXiv preprint arXiv:1611.01855, 2016.
  • [42] Oleksandr Polozov and Sumit Gulwani. Flashmeta: a framework for inductive program synthesis. In Proceedings of the 2015 ACM SIGPLAN International Conference on Object-Oriented Programming, Systems, Languages, and Applications, OOPSLA 2015, part of SPLASH 2015, Pittsburgh, PA, USA, October 25-30, 2015, pages 107–126, 2015.
  • [43] Veselin Raychev, Pavol Bielik, Martin T. Vechev, and Andreas Krause. Learning programs from noisy data. In Proceedings of the 43rd Annual ACM SIGPLAN-SIGACT Symposium on Principles of Programming Languages, POPL 2016, St. Petersburg, FL, USA, January 20 - 22, 2016, pages 761–774, 2016.
  • [44] Stephane Ross and J Andrew Bagnell. Reinforcement and imitation learning via interactive no-regret learning. arXiv preprint arXiv:1406.5979, 2014.
  • [45] Stéphane Ross, Geoffrey Gordon, and Drew Bagnell. A reduction of imitation learning and structured prediction to no-regret online learning. In Proceedings of the fourteenth international conference on artificial intelligence and statistics, pages 627–635, 2011.
  • [46] Stéphane Ross, Geoffrey J. Gordon, and Drew Bagnell. A reduction of imitation learning and structured prediction to no-regret online learning. In Proceedings of the Fourteenth International Conference on Artificial Intelligence and Statistics, AISTATS 2011, Fort Lauderdale, USA, April 11-13, 2011, pages 627–635, 2011.
  • [47] John Schulman, Sergey Levine, Pieter Abbeel, Michael Jordan, and Philipp Moritz. Trust region policy optimization. In International Conference on Machine Learning, pages 1889–1897, 2015.
  • [48] John Schulman, Filip Wolski, Prafulla Dhariwal, Alec Radford, and Oleg Klimov. Proximal policy optimization algorithms. arXiv preprint arXiv:1707.06347, 2017.
  • [49] Xujie Si, Yuan Yang, Hanjun Dai, Mayur Naik, and Le Song. Learning a meta-solver for syntax-guided program synthesis. In 7th International Conference on Learning Representations, ICLR 2019, New Orleans, LA, USA, May 6-9, 2019, 2019.
  • [50] Armando Solar-Lezama, Liviu Tancau, Rastislav Bodík, Sanjit A. Seshia, and Vijay A. Saraswat. Combinatorial sketching for finite programs. In ASPLOS, pages 404–415, 2006.
  • [51] Wen Sun, J Andrew Bagnell, and Byron Boots. Truncated horizon policy search: Combining reinforcement learning & imitation learning. In International Conference on Learning Representations (ICLR), 2018.
  • [52] Wen Sun, Geoffrey J Gordon, Byron Boots, and J Bagnell. Dual policy iteration. In Advances in Neural Information Processing Systems, pages 7059–7069, 2018.
  • [53] Wen Sun, Arun Venkatraman, Geoffrey J Gordon, Byron Boots, and J Andrew Bagnell. Deeply aggrevated: Differentiable imitation learning for sequential prediction. In International Conference on Machine Learning (ICML), 2017.
  • [54] Richard S Sutton and Andrew G Barto. Reinforcement learning: An introduction. MIT press, 2018.
  • [55] Richard S Sutton, David A McAllester, Satinder P Singh, and Yishay Mansour. Policy gradient methods for reinforcement learning with function approximation. In Advances in neural information processing systems, pages 1057–1063, 2000.
  • [56] Philip S Thomas, William C Dabney, Stephen Giguere, and Sridhar Mahadevan. Projected natural actor-critic. In Advances in neural information processing systems, pages 2337–2345, 2013.
  • [57] Lazar Valkov, Dipak Chaudhari, Akash Srivastava, Charles Sutton, and Swarat Chaudhuri. Houdini: Lifelong learning as program synthesis. In Advances in Neural Information Processing Systems, pages 8687–8698, 2018.
  • [58] Abhinav Verma, Vijayaraghavan Murali, Rishabh Singh, Pushmeet Kohli, and Swarat Chaudhuri. Programmatically interpretable reinforcement learning. In International Conference on Machine Learning, pages 5052–5061, 2018.
  • [59] Bernhard Wymann, Eric Espié, Christophe Guionneau, Christos Dimitrakakis, Rémi Coulom, and Andrew Sumner. TORCS, The Open Racing Car Simulator. http://www.torcs.org, 2014.
  • [60] He Zhu, Zikang Xiong, Stephen Magill, and Suresh Jagannathan. An inductive synthesis framework for verifiable reinforcement learning. In ACM Conference on Programming Language Design and Implementation (SIGPLAN), 2019.

Appendix A Theoretical Analysis

A.1 Preliminaries and Notations

We formally define an ambient control policy space 𝒰 to be a vector space equipped with inner product ,:𝒰×𝒰, which induces a norm u=u,u, and its dual norm defined as v*=sup{v,u|u1}. While multiple ways to define the inner product exist, for concreteness we can think of the example of square-integrable stationary policies with u,v=𝒮u(s)v(s)𝑑s. The addition operator + between two policies u,v𝒰 is defined as (u+v)(s)=u(s)+v(s) for all state s𝒮. Scaling λu+κv is defined similarly for scalar λ,κ.

The cost functional of a control policy u is defined as J(u)=0c(s(τ),u(τ))𝑑τ, or J(u)=𝒮c(s,u(s))𝑑μu(s), where μu is the distribution of states induced by policy u. This latter example is equivalent to the standard notion of value function in reinforcement learning.

Separate from the parametric representation issues, both programmatic policy class Π and neural policy class , and by extension - the joint policy class , are considered to live in the ambient vector space 𝒰. We thus have a common and well-defined notion of distance between policies from different classes.

We make an important distinction between differentiability of J(h) in the ambient policy space (non-parametric), versus differentiability in parameterization (parametric). For example, if Π is a class of decision-tree based policy, policies in Π may not be differentiable under representation. However, policies πΠ might still be differentiable when considered as points in the ambient vector space 𝒰.

We will use the following standard notion of gradient and differentiability from functional analysis:

Definition A.1 (Subgradients).

The subgradient of J at h, denoted J(h), is the non-empty set {g|j:j-h,g+J(h)J(j)}

Definition A.2 (Fréchet gradient).

A bounded linear operator : is called Fréchet functional gradient of J at h if limg0J(h+g)-J(h)-J(h),gg=0

The notions of convexity, smoothness and Bregman divergence are analogous to finite-dimensional setting:

Definition A.3 (Strong convexity).

A differentiable function R is α-strongly convex w.r.t norm if R(y)R(x)+R(x),y-x+α2y-x2

Definition A.4 (Lipschitz continuous gradient smoothness).

A differentiable function R is LR-strongly smooth w.r.t norm if R(x)-R(y)*LRx-y

Definition A.5 (Bregman Divergence).

For a strongly convex regularizer R, DR(x,y)=R(x)-R(y)-R(y),x-y is the Bregman divergence between x and y (not necessarily symmetric)

The following standard result for Bregman divergence will be useful:

Lemma A.1.

[10] For all x,y,z we have the identity R(x)-R(y),x-z=DR(x,y)+DR(z,x)-DR(z,y). Since Bregman divergence is non-negative, a consequence of this identity is that DR(z,x)-DR(z,y)R(x)-R(y),z-x

A.2 Expected Regret Bound under Noisy Policy Gradient Estimates and Projection Errors

In this section, we show regret bound for the performance of the sequence of returned programs π1,,πT of the algorithm. The analysis here is agnostic to the particular implementation of algorithm 3 and algorithm 3.

Let R be a α-strongly convex and LR-smooth functional with respect to norm on . The steps from algorithm 2 can be described as follows.

  • Initialize π0Π. For each iteration t:

    1. 1.

      Obtain a noisy estimate of the gradient ^J(πt-1)J(πt-1)

    2. 2.

      Update in the space: R(ht)=R(πt-1)-η^J(πt-1)

    3. 3.

      Obtain approximate projection πt=ProjectπR(ht)argminpiΠDR(π,ht)

This procedure is an approximate functional mirror descent scheme under bandit feedback. We will develop the following result, which is a more detailed version of 4.1 in the main paper.

In the statement below, D is the diameter on Π with respect to defined norm (i.e., D=supπ-π). LJ is the Lipschitz constant of the functional J on . β,σ2 are the bound on the bias and variance of the gradient estimate at each iteration, respectively. α and ŁR are the strongly convex and smooth coefficients of the functional regularizer R. Finally, ϵ is the bound on the projection error with respect to the same norm .

Theorem A.2 (Regret bound of returned policies).

Let π1,,πT be a sequence of programmatic policies returned by algorithm 2 and π* be the optimal programmatic policy. We have the expected regret bound:

𝔼[1Tt=1TJ(πt)]-J(π*)LRD2ηT+ϵLRDη+η(σ2+LJ2)α+βD

In particular, choosing the learning rate η=1T+ϵσ2, the expected regret is simplified into:

(4)
𝔼[1Tt=1TJ(πt)]-J(π*)=O(σ1T+ϵ+β)
Proof.

At each round t, let \[email protected]Δ\[email protected]\[email protected]\[email protected]@skewchar\[email protected]@a111t=𝐄[^t|πt] be the conditional expectation of the gradient estimate. We will use the shorthand notation t=J(πt). Denote the upper-bound on the bias of the estimate by βt, i.e., \[email protected]Δ\[email protected]\[email protected]\[email protected]@skewchar\[email protected]@a111t-t*βt almost surely. Denote the noise of the gradient estimate by ξt=\[email protected]Δ\[email protected]\[email protected]\[email protected]@skewchar\[email protected]@a111t-^t, and σt2=𝐄[^t-\[email protected]Δ\[email protected]\[email protected]\[email protected]@skewchar\[email protected]@a111t*2] is the variance of gradient estimate ^t.

The projection operator is ϵ-approximate in the sense that πt-ProjectΠR(ft)=Project^ΠR(ht)-ProjectΠR(ht)ϵ with some constant ϵ, which reflects the statistical error of the imitation learning procedure. This projection error in general is independent of the choice of function classes Π and .We will use the shorthand notation πt*=ProjectΠR(ft) for the true Bregman projection of ht onto Π.

Due to convexity of J over the space (which includes Π), we have for all πΠ:

J(πt)-J(π)t,πt-π

We proceed to bound the RHS, starting with bounding the inner product where the actual gradient is replaced by the estimated gradient.

^t,πt-π=1ηtR(πt)-R(ht+1),πt-π (5)
=1ηt(DR(π,πt)-DR(π,ht+1)+DR(πt,ht+1)) (6)
1ηt(DR(π,πt)-DR(π,πt+1*)-DR(πt+1*,ht+1)+DR(πt,ht+1)) (7)
=1ηt(DR(π,πt)-DR(π,πt+1)telescoping+DR(π,πt+1)-DR(π,πt+1*)projection error-DR(πt+1*,ht+1)+DR(πt,ht+1)relative improvement) (8)

Equation (5) is due to the gradient update rule in space. Equation (6) is derived from definition of Bregman divergence. Equation (7) is due to the generalized Pythagorean theorem of Bregman projection DR(x,y)DR(x,ProjectΠR(x))+DR(ProjectΠR(x),y). The RHS of equation (7) are decomposed into three components that will be bounded separately.

Bounding projection error. By lemma (A.1) we have

DR(π,πt+1)-DR(π,πt+1*)R(πt+1)-R(πt+1*),π-πt+1 (9)
R(πt+1)-R(πt+1*)π-πt+1* (10)
LRπt+1-πt+1*DϵLRD (11)

Equation (10) is due to Cauchy–Schwarz. Equation (11) is due to Lipschitz smoothness of R and definition of ϵ-approximate projection.

Bounding relative improvement. This follows standard argument from analysis of mirror descent algorithm.

DR(πt,ht+1)-DR(πt+1*,ht+1)=R(πt)-R(πt+1*)+R(ht+1),πt+1*-πt (12)
R(πt),πt-πt+1*-α2πt+1*-πt*2+R(ht+1),πt+1*-πt (13)
=-ηt^t,πt+1*-πt-α2πt+1*-πt2 (14)
ηt22α^t*2ηt2α(σt2+LJ2) (15)

Equation (13) is from the α-strong convexity property of regularizer R. Equation (14) is by definition of the gradient update. Combining the bounds on the three components and taking expectation, we thus have

𝔼[^t,πt-π]1ηt(DR(π,πt)-DR(π,πt+1)+ϵLRD+ηt2α(σt2+LJ2)) (16)

Next, the difference between estimated gradient ^t and actual gradient t factors into the bound via Cauchy-Schwarz:

𝔼[t-^t,πt-π]t-𝔼[^t]*πt-πβtD (17)

The results can be deduced from equations (16) and (17).

Unbiased gradient estimates. For the case when the gradient estimate is unbiased, assume the variance of the noise of gradient estimates is bounded by σ2, we have the expected regret bound for all piΠ

𝔼[1Tt=1TJ(πt)]-J(π)LRD2ηT+ϵLRDη+η(σ2+LJ2)α (18)

here to clarify, LR is the smoothness coefficient of regularizer R (i.e., the gradient of R is LR-Lipschitz, LJ is Lipschitz constant of J, D is the diameter of Π under norm , σ2 is the upper-bound on the variance of gradient estimates, and ϵ is the error from the projection procedure (i.e., imitation learning loss).

We can set learning rate η=1T+ϵσ2 to observe that the expected regret is bounded by O(σ1T+ϵ).

Biased gradient estimates. Assume that the bias of gradient estimate at each round is upper-bounded by βtβ. Similar to before, combining inequalities from (16) and (17), we have

𝔼[1Tt=1TJ(πt)]-J(π)LRD2ηT+ϵLRDη+η(σ2+LJ2)α+βD (19)

Similar to before, we can set learning rate η=1T+ϵσ2 to observe that on the expected regret is bounded by O(σ1T+ϵ+β). Compared to the bound on (18), in the biased case, the extra regret incurred per bound is simply a constant, and does not depend on T. ∎

A.3 Finite-Sample Analysis

In this section, we provide overall finite-sample analysis for Propel under some simplifying assumptions. We first consider the case where exact gradient estimate is available, before extending the result to the general case of noisy policy gradient update. Combining the two steps will give us the proof for the following statement (theorem 4.2 in the main paper)

Theorem A.3 (Finite-sample guarantee).

At each iteration, we perform vanilla policy gradient estimate of π (over H) using m trajectories and use DAgger algorithm to collect M roll-outs. Setting the learning rate η=1σ2(1T+HM+log(T/δ)M), after T rounds of the algorithm, we have that

1Tt=1TJ(πt)-J(π*)O(σ1T+HM+log(T/δ)M)+O(σlog(Tk/δ)m+AHlog(Tk/δ)m)

holds with probability at least 1-δ, with H the task horizon, A the cardinality of action space, σ2 the variance of policy gradient estimates, and k the dimension Π’s parameterization.

Exact gradient estimate case. Assuming that the policy gradients can be calculated exactly, it is straight-forward to provide high-probability guarantee for the effect of the projection error. We start with the following result, adapted from [45] for the case of projection error bound. In this version of DAgger, we assume that we only collect a single (state, expert action) pair from each trajectory roll-out. Result is similar, with tighter bound, when multiple data points are collected along the trajectory.

Lemma A.4 (Projection error bound from imitation learning procedure).

Using DAgger as the imitation learning sub-routine for our Project operator in algorithm 3, let M be the number of trajectories rolled-out for learning, and H be the horizon of the task. With probability at least 1-δ, we have

DR(π,π*)O~(1/M)+2max(1+H)M+2maxlog(1/δ))M

where π is the result of Project, π* is the true Bregman projection of h onto Π, and max is the maximum value of the imitation learning loss function DR(,)

The bound in lemma A.4 is simpler than previous imitation learning results with cost information ([44, 45]. The reason is that the goal of the Project operator is more modest. Since we only care about the distance between the empirical projection π and the true projection π*, the loss objective in imitation learning is simplified (i.e., this is only a regret bound), and we can disregard how well policies in Π can imitate the expert h, as well as the performance of J(π) relative to the true cost from the environment J(h).

A consequence of this lemma is that for the number of trajectories at each round of imitation learning M=O(log1/δϵ2)+O(Hϵ), we have DR(πt,πt*)ϵ with probability at least 1-δ. Applying union bound across T rounds of learning, we obtain the following guarantee (under no gradient estimation error)

Proposition A.5 (Finite-sample Projection Error Bound).

To simplify the presentation of the result, we consider LR,D,L,α to be known constants. Using DAgger algorithm to collect M=O(logT/δϵ2)+O(Hϵ) roll-outs at each iteration, we have the following regret guarantee after T rounds of our main algorithm:

1Tt=1TJ(πt)-J(π*)O(1ηT+ϵη+η)

with probability at least 1-δ. Consequently, setting η=1T+HM+log(T/δ)M, we have that

1Tt=1TJ(πt)-J(π*)O(1T+HM+log(T/δ)M)

with probability at least 1-δ

Note that the dependence on the time horizon of the task is sub-linear. This is different from standard imitation learning regret bounds, which are often at least linear in the task horizon. The main reason is that our comparison benchmark π* does live in the space Π, whereas for DAgger, the expert policy may not reside in the same space.

Noisy gradient estimate case. We now turn to the issue of estimating the gradient of J(π). We make the following simplifying assumption about the gradient estimation:

  • The π is parameterized by vector θk (such as a neural network). The parameterization is differentiable with respect to θ (Alternatively, we can view Π as a differentiable subspace of , in which case we have =)

  • At each Update loop, the policy is rolled out m times to collect the data, each trajectory has horizon length H

  • For each visited state sdh, the policy takes a uniformly random action a. The action space is finite with cardinality A.

  • The gradient hθ is bounded by B

The gradient estimate is performed consistent with a generic policy gradient scheme, i.e.,

^J(θ)=Ami=1Hj=1mπθ(aij|sij,θ)Q^ij

where Q^ij is the estimated cost-to-go [55].

Taking uniform random exploratory actions ensures that the samples are i.i.d. We can thus apply Bernstein’s inequality to obtain the bound between estimated gradient and the true gradient. Indeed, with probability at least 1-δ, we have that the following bound on the bias component-wise:

^J(θ)-J(θ)β when m(2σ2+2AHBβ3)logkδβ2

which leads to similar bound with respect to * (here we leverage the equivalence of norms in finite dimensional setting):

t-^t*β when m=O((σ2+AHBβ)logkδβ2)

Applying union bound of this result over T rounds of learning, and combining with the result from proposition (A.5), we have the following finite-sample guarantee in the simplifying policy gradient update. This is also the more detailed statement of theorem 4.2 in the main paper.

Proposition A.6 (Finite-sample Guarantee under Noisy Gradient Updates and Projection Error).

At each iteration, we perform policy gradient estimate using m=O((σ2+AHBβ)logTkδβ2) trajectories and use DAgger algorithm to collect M=O(logT/δϵ2)+O(Hϵ) roll-outs. Setting the learning rate η=1σ2(1T+HM+log(T/δ)M), after T rounds of the algorithm, we have that

1Tt=1TJ(πt)-J(π*)O(σ1T+HM+log(T/δ)M)+β

with probability at least 1-δ.

Consequently, we also have the following regret bound:

1Tt=1TJ(πt)-J(π*)O(σ1T+HM+log(T/δ)M)+O(σlog(Tk/δ)m+AHlog(Tk/δ)m)

holds with probability at least 1-δ, where again H is the task horizon, A is the cardinality of action space, and k is the dimension of function class Π’s parameterization.

Proof.

(For both proposition (A.6) and (A.5)). The results follow by taking the inequality from equation (19), and by solving for ϵ and β explicitly in terms of relevant quantities. Based on the specification of M and m, we obtain the necessary precision for each round of learning in terms of number of trajectories:

β =O(σlog(k/δ)m+AHBlog(k/δ)m)
ϵ =O(HM+log(1/δ)M)

Setting the learning rate η=1σ2(1T+ϵ) and rearranging the inequalities lead to the desired bounds. ∎

The regret bound depends on the variance σ2 of the policy gradient estimates. It is well-known that vanilla policy gradient updates suffer from high variance. We instead use functional regularization technique, based on CORE-RL, in the practical implementation of our algorithm. The CORE-RL subroutine has been demonstrated to reduce the variance in policy gradient updates [19].

A.4 Defining a consistent approximation of J(π) - Proof of Proposition 4.3

We are using the notion of Fréchet derivative to define gradient of differentiable functional. Note that while Gateaux derivative can also be utilized, Fréchet derivative ensures continuity of the gradient operator that would be useful for our analysis.

Definition A.6 (Fréchet gradient).

A bounded linear operator : is called Fréchet functional gradient of J at h if limg0J(h+g)-J(h)-J(h),gg=0

We make the following assumption about and . One interpretation of this assumption is that the space of policies Π and that we consider have the property that a programmatic policy πΠ can be well-approximated by a large space of neural policies f.

Assumption 1.

J is Fréchet differentiable on H. J is also differentiable on the restricted subspace F. And F is dense in H (i.e., the closure \[email protected]Δ\[email protected]\[email protected]\[email protected]@skewchar\[email protected]@a111F=H)

It is then clear that f the Fréchet gradient J(f), restricted to the subspace is equal to the gradient of f in the ambient space (since Fréchet gradient is unique). In general, given πΠ and f, π+f is not necessarily in . However, the restricted gradient on subspace of J(π+f) can be defined asymptotically.

Proposition A.7.

Fixing a policy πΠ, define a sequence of policies fkF, k=1,2, that converges to π: limkfk-g=0, we then have limkFJ(fk)-HJ(π)*=0

Proof.

Since Fréchet derivative is a continuous linear operator, we have limkJ(fk)-J(π)*=0. By the reasoning above, for f, the gradient J(f) defined via restriction to the space does not change compared to J(f), the gradient defined over the ambient space . Thus we also have limkJ(fk)-J(π)*=0. By the same argument, we also have that for any given πΠ and f, even if π+f, the gradient J(π+f) with respect to the can be approximated similarly. ∎

Note that we are not assuming J(π) to be differentiable when restricting to the policy subspace Π.

A.5 Theoretical motivation for Algorithm 3 - Proof of Proposition 4.4 and 4.5

We consider the case where Π is not differentiable by parameterization. Note that this does not preclude J(π) for πΠ to be differentiable in the non-parametric function space. Two complications arise compared to our previous approximate mirror descent procedure. First, for each πΠ, estimating the gradient J(π) (which may not exist under certain parameterization, per section 4.3) can become much more difficult. Second, the update rule R(π)-J(π) may not be in the dual space of , as in the simple case where Π, thus making direct gradient update in the space inappropriate.

Assumption 2.

J is convex in H.

By convexity of J in , sub-gradients J(h) exists for all h. In particular, J(π) exists for all πΠ. Note that J(π) reflects sub-gradient of π with respect to the ambient policy space .

We will make use of the following equivalent perspective to mirror descent[10], which consists of two-step process for each iteration t

  1. 1.

    Solve for ht+1=argminhηJ(πt),h+DR(h,πt)

  2. 2.

    Solve for πt+1=argminπΠDR(π,ht+1)

We will show how this version of the algorithm motivates our main algorithm. Consider step 1 of the main loop of Propel, where given a fixed πΠ, the optimization problem within is

(OBJECTIVE_1)=minhηJ(π),h+DR(h,π) (20)

Due to convexity of and the objective, problem (OBJECTIVE_1) is equivalent to:

(OBJECTIVE_1)= minJ(π),h (21)
s.t. DR(h,π)τ (22)

where τ depends on η. Since π is fixed, this optimization problem can be relaxed by choosing λ[0,1], and a set of candidate policies h=π+λf, for all f, such that DR(h,π)τ is satisfied (Selection of λ is possible with bounded spaces). Since this constraint set is potentially a restricted set compared to the space of policies satisfying inequality (22), the optimization problem (20) is relaxed into:

(OBJECTIVE_1)(OBJECTIVE_2)=minfJ(π),π+λf (23)

Due to convexity property of J, we have

J(π),λf=J(π),π+λf-π)J(π+λf)-J(π) (24)

The original problem OBJECTIVE_1 is thus upper bounded by:

minhηJ(π),h)+DR(h,π)minfJ(π+λf)-J(π)+J(π),π

Thus, a relaxed version of original optimization problem OBJECTIVE_1 can be obtained by miniziming J(π+λf) over f (note that π is fixed). This naturally motivates using functional regularization technique, such as CORE-RL algorithm [19], to update the parameters of differentiable function f via policy gradient descent update:

f=f-ηλλJ(π+λf)

where the gradient of J is taken with respect to the parameters of f (neural networks). This is exactly the update step in algorithm 3 (also similar to iterative updte of CORE-RL algorithm), where the neural network policy is regularized by a prior controller π.

Statement and Proof of Proposition 4.5

Proposition A.8 (Regret bound for the relaxed optimization objective).

Assuming J(h) is L-strongly smooth over H, i.e., HJ(h) is L-Lipschitz continuous, approximating UpdateH by UpdateF per Alg. 3 leads to the expected regret bound: E[1Tt=1TJ(πt)]-J(π*)=O(λσ1T+ϵ+λ2L2)

Proof.

Instead of focusing on the bias of the gradient estimate J(π), we will shift our focus on the alternative proximal formulation of mirror descent, under optimization and projection errors. In particular, at each iteration t, let ht+1*=argminhηJ(πt),h+DR(h,πt) and let the optimization error be defined as βt where R(ht+1)=R(ht+1*)+βt. Note here that this is different from (but related to) the notion of bias from gradient estimate of J(π) used in theorem 4.1 and theorem A.2. The projection error from imitation learning procedure is defined similarly to theorem 4.1: πt+1*=argminπΠDR(π,ht+1) is the true projection, and πt+1-πt+1*ϵ.

We start with similar bounding steps to the proof of theorem 4.1:

J(πt),πt-π =1ηR(ht+1*)-R(πt),πt-π
=1η(R(ht+1)-R(πt),πt-π-βt,πt-π)
=1η(DR(π,πt)-DR(π,ht+1)+DR(πt,ht+1))component_1+1ηβt,πt-πcomponent_2 (25)

As seen from the proof of theorem A.2, component_1 can be upperbounded by: 1η(DR(π,πt)-DR(π,πt+1)telescoping+DR(π,πt+1)-DR(π,πt+1*)projection error-DR(πt+1*,ht+1)+DR(πt,ht+1)relative improvement) The bound on projection error is identical to theorem A.2:

DR(π,πt)-DR(π,πt+1*)ϵLRD (26)

The bound on relative improvement is slightly different:

DR(πt,ht+1)-DR(πt+1*,ht+1)=R(πt)-R(πt+1*)+R(ht+1),πt+1*-πt
=R(πt)-R(πt+1*+R(ht+1*),πt+1*-πt)+βt,πt+1*-πt
R(πt),πt-πt+1*-α2πt+1*-πt2+R(ht+1*),πt+1*-πt+βt,πt+1*-πt
=-ηJ(πt),πt+1*-πt-α2πt+1*-πt2+βt,πt+1*-πt (27)
η22αJ(πt)*2+βt,πt+1*-πt
η22αLJ2+βt,πt+1*-πt (28)

Note here that the gradient J(πt) is not the result of estimation. Combining equations (25), (26), (27), (28), we have:

J(πt),πt-π1η(DR(π,πt)-DR(π,πt+1)+ϵLRD+η22αLJ2+βt,πt+1*-π) (29)

Next, we want to bound βt. Choose regularizer R to be 122 (consistent with the pseudocode in algorithm 3). We have that:

ht+1*=argminhηJ(πt),h+12h-πt2

which is equivalent to:

ht+1*=πt+argminfηJ(πt),f+12f2

Let ft+1*=argminfηJ(πt),f+12f2. Taking the gradient over f, we can see that ft+1*=-ηJ(πt). Let ft+1 be the minimizer of minfJ(πt+λf). We then have ht+1*=πt+ft+1* and ht+1=π+λft+1. Thus βt=ht+1-ht+1*=ft+1-ft+1*.

On one hand, we have

J(πt+λft+1)J(πt+ωft+1*)J(πt)+J(πt),ωft+1*+L2ωft+1*2

due to optimality of ft+1 and strong smoothness property of J. On the other hand, since J is convex, we also have the first-order condition:

J(πt+λft+1)J(πt)+J(πt),λft+1

Combine with the inequality above, and subtract J(πt) from both sides, and using the relationship ft+1*=-ηJ(πt), we have that:

-1ηft+1*,λft+1-1ηft+1*,ωft+1*+Lω22ft+1*2

Since this is true ω, rearrange and choose ω such that ωη-Lω22=-λ2η, namely ω=1-1-ληLLη, and complete the square, we can establish the bound that:

ft+1-ft+1*η(λL)2B (30)

for B the upperbound on ft+1. We thus have βt=O(η(λL)2). Plugging the result from equation 30 into RHS of equation 29, we have:

J(πt),πt-π1η(DR(π,πt)-DR(π,πt+1)+ϵLRD+η22αLJ2)+(η(λL)2B) (31)

Since J is convex in , we have J(πt)-J(π)J(πt),πt-π. Similar to theorem 4.1, setting η=1λ2σ2(1T+ϵ) and taking expectation on both sides, we have:

𝔼[1Tt=1TJ(πt)]-J(π*)=O(λσ1T+ϵ+λ2L2) (32)

Note that unlike regret bound from theorem 4.1 under general bias, variance of gradient estimate and projection error, σ2 here explicitly refers to the bound on neural-network based policy gradient variance. The variance reduction of λσ, at the expense of some bias, was also similarly noted in a recent functional regularization technique for policy gradient [19]. ∎

Appendix B Additional Experimental Results and Details

B.1 TORCS

We generate controllers for cars in The Open Racing Car Simulator (Torcs) [59]. In its full generality Torcs provides a rich environment with input from up to 89 sensors, and optionally the 3D graphic from a chosen camera angle in the race. The controllers have to decide the values of 5 parameters during game play, which correspond to the acceleration, brake, clutch, gear and steering of the car.

Apart from the immediate challenge of driving the car on the track, controllers also have to make race-level strategy decisions, like making pit-stops for fuel. A lower level of complexity is provided in the Practice Mode setting of TORCS. In this mode all race-level strategies are removed. Currently, so far as we know, state-of-the-art Drl models are capable of racing only in Practice Mode, and this is also the environment that we use. Here we consider the input from 29 sensors, and decide values for the acceleration, steering, and braking actions.

We chose a suite of tracks that provide varying levels of difficulty for the learning algorithms. In particular, for the tracks Ruudskogen and Alpine-2, the Ddpg agent is unable to reliably learn a policy that would complete a lap. We perform the experiments with twenty-five random seeds and report the median lap time over these twenty-five trials. However we note that the Torcs simulator is not deterministic even for a fixed random seed. Since we model the environment as a Markov Decision Process, this non-determinism is consistent with our problem statement.

For our Deep Reinforcement Learning agents we used standard open source implementations (with pre-tuned hyper-parameters) for the relevant domain.

All experiments were conducted on standard workstation with a 2.5 GHz Intel Core i7 CPU and a GTX 1080 Ti GPU card.

The code for the Torcs experiments can be found at: https://bitbucket.org/averma8053/propel

In Table 3 we show the lap time performance and crash ratios of Propel agents initialized with neural policies obtained via Ddpg. As discussed in Section 5, Ddpg often exhibits high variance across trials and this adversely affects the performance of the Propel agents when they are initialized via Ddpg. In Table 4 we show generalization results for the PropelTree agent. As noted in Section 5, the generalization results for PropelTree are in between those of Ddpg and PropelProg.

Verified Smoothness Property. For the program given in Figure 2 we proved using symbolic verification techniques, that k,i=kk+5𝐩𝐞𝐞𝐤(s[𝚁𝙿𝙼],i+1)-𝐩𝐞𝐞𝐤(s[𝚁𝙿𝙼],i)<0.003𝐩𝐞𝐞𝐤(a[𝙰𝚌𝚌𝚎𝚕],k+1)-𝐩𝐞𝐞𝐤(a[𝙰𝚌𝚌𝚎𝚕],k)<0.63. Here the function 𝐩𝐞𝐞𝐤(.,i) takes in a history/sequence of sensor or action values and returns the value at position i, . Intuitively, the above logical implication means that if the sum of the consecutive differences of the last six 𝚁𝙿𝙼 sensor values is less than 0.003, then the acceleration actions calculated at the last and penultimate step will not differ by more than 0.63.

Table 3: Performance results in Torcs of Propel agents initialized with neural policies obtained via Ddpg, over 25 random seeds. Each entry is formatted as Lap-time / Crash-ratio, reporting median lap time in seconds over all the seeds (lower is better) and ratio of seeds that result in crashes (lower is better). A lap time of Cr indicates the agent crashed and could not complete a lap for more than half the seeds.
G-Track E-Road Aalborg Ruudskogen Alpine-2
Length 3186m 3260m 2588m 3274m 3774m
PropelProg-Ddpg 97.76/.12 108.06/.08 140.48/.48 Cr / 0.68 Cr / 0.92
PropelTree-Ddpg 78.47/0.16 85.46/.04 Cr / 0.56 Cr / 0.68 Cr / 0.92
Table 4: Generalization results in Torcs for PropelTree, where rows are training and columns are testing tracks. Each entry is formatted as PropelProg / DDPG, and the number reported is the median lap time in seconds over all the seeds (lower is better). Cr indicates the agent crashed and could not complete a lap for more than half the seeds.
G-Track E-Road Aalborg Ruudskogen Alpine-2
G-Track - 95 Cr Cr Cr
E-Road 84 - Cr Cr Cr
Aalborg 111 Cr - Cr Cr
Ruudskogen 154 Cr Cr - Cr
Alpine-2 Cr 276 Cr Cr -
Table 5: Performance results in Classic Control problems. Higher scores are better.
MountainCar Pendulum
Prior 00.59±0.00 -875.53±0.00
Ddpg 97.16±3.21 -132.70±6.44
Trpo 93.03±1.86 -131.54±4.56
Ndps 66.98±3.11 -435.71±4.83
Viper 64.86±3.28 -394.11±4.97
PropelProg 95.63±1.02 -187.71±2.35
PropelTree 96.56±2.81 -139.09±3.31

B.2 Classic Control

We present results from two classic control problems, Mountain-Car (with continuous actions) and Pendulum, in Table 5. We use the OpenAI Gym implementations of these environments. More information about these environments can be found at the links: MountainCar and Pendulum.

In Mountain-Car the goal is to drive an under-powered car up the side of a mountain in as few time-steps as possible. In Pendulum, the goal is to swing a pendulum up so that it stays upright. In both the environments an episode terminates after a maximum of 200 time-steps.

In Table 5 we report the mean and standard deviation, over twenty-five random seeds, of the average scores over 100 episodes for the listed agents and environments. In Figure 7 and Figure 7 we show the improvements made over the prior by the PropelProg agent in MountainCar and Pendulum respectively, with each iteration of the Propel algorithm.

Figure 6: Score improvements in the MountainCar environment over iterations of PropelProg.
Figure 7: Score improvements in the Pendulum environment over iterations of PropelProg.