Contrastive Behavioral Similarity
Embeddings for Generalization in
Reinforcement Learning

Abstract

Reinforcement learning methods trained on few environments rarely learn policies that generalize to unseen environments. To improve generalization, we incorporate the inherent sequential structure in reinforcement learning into the representation learning process. This approach is orthogonal to recent approaches, which rarely exploit this structure explicitly. Specifically, we introduce a theoretically motivated policy similarity metric (PSM) for measuring behavioral similarity between states. PSM assigns high similarity to states for which the optimal policies in those states as well as in future states are similar. We also present a contrastive representation learning procedure to embed any state similarity metric, which we instantiate with PSM to obtain policy similarity embeddings (PSEs). We demonstrate that PSEs improve generalization on diverse benchmarks, including LQR with spurious correlations, a jumping task from pixels, and Distracting DM Control Suite. Source code is available at agarwl.github.io/pse.

1 Introduction

Figure 1: Jumping task: The agent (white block), learning from pixels, needs to jump over an obstacle (grey square). The challenge is to generalize to unseen obstacle positions and floor heights in test tasks using a small number of training tasks. We show the agent’s trajectories using faded blocks.

Current reinforcement learning (RL) approaches often learn policies that do not generalize to environments different than those the agent was trained on, even when these environments are semantically equivalent (Tachet des Combes et al., 2018; Song et al., 2019; Cobbe et al., 2019). For example, consider a jumping task where an agent, learning from pixels, needs to jump over an obstacle (Figure 1). Deep RL agents trained on a few of these jumping tasks with different obstacle positions struggle to successfully jump in test tasks where obstacles are at previously unseen locations.

Recent solutions to circumvent poor generalization in RL are adapted from supervised learning, and, as such, largely ignore the sequential aspect of RL. Most of these solutions revolve around enhancing the learning process, including data augmentation (e.g., Kostrikov et al., 2020; Lee et al., 2020a), regularization (Cobbe et al., 2019; Farebrother et al., 2018), noise injection (Igl et al., 2019), and diverse training conditions (Tobin et al., 2017); they rarely exploit properties of the sequential decision making problem such as similarity in actions across temporal observations.

Instead, we tackle generalization by incorporating properties of the RL problem into the representation learning process. Our approach exploits the fact that an agent, when operating in environments with similar underlying mechanics, exhibits at least short sequences of behaviors that are similar across these environments. Concretely, the agent is optimized to learn an embedding in which states are close when the agent’s optimal policies in these states and future states are similar. This notion of proximity is general and it is applicable to observations from different environments.

Specifically, inspired by bisimulation metrics (Castro, 2020; Ferns et al., 2004), we propose a novel policy similarity metric (PSM). PSM (Section 3) defines a notion of similarity between states originated from different environments by the proximity of the long-term optimal behavior from these states. PSM is reward-agnostic, making it more robust for generalization compared to approaches that rely on reward information. We prove that PSM yields an upper bound on suboptimality of policies transferred from one environment to another (Theorem 1), which is not attainable with bisimulation.

We employ PSM for representation learning and introduce policy similarity embeddings (PSEs) for deep RL. To do so, we present a general contrastive procedure (Section 4) to learn an embedding based on any state similarity metric. PSEs are the instantiation of this procedure with PSM. PSEs are appealing for generalization as they encode task-relevant invariances by putting behaviorally equivalent states together. This is unlike prior approaches, which rely on capturing such invariances without being explicitly trained to do so, for example, through value function similarities across states (e.g., Castro & Precup, 2010), or being robust to fixed transformations of the observation space (e.g., Kostrikov et al., 2020; Laskin et al., 2020a).

PSEs lead to better generalization while being orthogonal to how most of the field has been tackling generalization. We illustrate the efficacy and broad applicability of our approach on three existing benchmarks specifically designed to test generalization: (i) jumping task from pixels (Tachet des Combes et al., 2018) (Section 5), (ii) LQR with spurious correlations (Song et al., 2019) (Section 6.1), and (iii) Distracting DM Control Suite (Stone et al., 2021) (Section 6.2). Our approach improves generalization compared to a wide variety of approaches including standard regularization (Farebrother et al., 2018; Cobbe et al., 2019), bisimulation (Castro & Precup, 2010; Castro, 2020; Zhang et al., 2021), out-of-distribution generalization (Arjovsky et al., 2019) and state-of-the-art data augmentation (Kostrikov et al., 2020; Laskin et al., 2020a; Lee et al., 2020a).

2 Preliminaries

We describe an environment as a Markov decision process (MDP) (Puterman, 1994) $ℳ=(𝒳,𝒜,R,P,\gamma )$, with a state space $𝒳$, an action space $𝒜$, a reward function $R$, transition dynamics $P$, and a discount factor $\gamma \in [0,1)$. A policy $\pi (\cdot |x)$ maps states $x\in 𝒳$ to distributions over actions. Whenever convenient, we abuse notation and write $\pi (x)$ to describe the probability distribution $\pi (\cdot |x)$, treating $\pi (x)$ as a vector. In RL, the goal is to find an optimal policy ${\pi }^{\ast }$ that maximizes the cumulative expected return ${𝔼}_{a_t\sim \pi (\cdot |x_t)}[{\sum }_t{\gamma }^{t}R(x_t,a_t)]$ starting from an initial state $x_0$.

We are interested in learning a policy that generalizes across related environments. We formalize this by considering a collection $\rho$ of MDPs, sharing an action space $𝒜$ but with disjoint state spaces. We use $𝒳$ and $𝒴$ to denote the state spaces of specific environments, and write $R_{𝒳}$, $P_{𝒳}$ for the reward and transition functions of the MDP whose state space is $𝒳$, and ${\pi }_{𝒳}^{\ast }$ for its optimal policy, which we assume unique without loss of generality. For a given policy $\pi$, we further specialize these into $R_{𝒳}^{\pi }$ and $P_{𝒳}^{\pi }$, the reward and state-to-state transition dynamics arising from following $\pi$ in that MDP.

We write $𝒮$ for the union of the state spaces of the MDPs in $\rho$. Concretely, different MDPs correspond to specific scenarios in a problem class (Figure 1), and $𝒮$ is the space of all possible configurations. Used without subscripts, $R$, $P$, and $\pi$ refer to the reward and transition function of this “union MDP”, and a policy defined over $𝒮$; this notation simplifies the exposition. We measure distances between states across environments using pseudometrics on $𝒮$; the set of all such pseudometrics is $𝕄$, and ${𝕄}_p$ is the set of metrics on probability distributions over $𝒮$.

In our setting, the learner has access to a collection of training MDPs ${\{{ℳ}_i\}}_{i=1}^N$, drawn from $\rho$. After interacting with these environments, the learner must produce a policy $\pi$ over the entire state space $𝒮$, which is then evaluated on unseen MDPs from $\rho$. Similar in spirit to the setting of transfer learning (Taylor & Stone, 2009), here we evaluate the policy’s zero-shot performance on $\rho$.

Our policy similarity metric (Section 3) builds on the concept of $𝝅$-bisimulation (Castro, 2020). Under the $\pi$-bisimulation metric, the distance between two states, $x$ and $y$, is defined in terms of the difference between the expected rewards obtained when following policy $\pi$. The $\pi$-bisimulation metric $d_{\pi }$ satisfies a recursive equation based on the 1-Wasserstein metric ${𝒲}_1:𝕄\to {𝕄}_p$, where ${𝒲}_1(d)(A,B)$ is the minimal cost of transporting probability mass from $A$ to $B$ (two probability distributions on $𝒮$) under the base metric $d$ (Villani, 2008). The recursion is

$d_{\pi }(x,y)=|R^{\pi }(x)-R^{\pi }(y)|+\gamma {𝒲}_1(d_{\pi })(P^{\pi }(\cdot |x),P^{\pi }(\cdot |y)),x,y\in 𝒮.$ (1)

To achieve good generalization properties, we learn an embedding function $z_{\theta }:𝒮\to {ℝ}^k$ that reflects the information encoded in the policy similarity metric; this yields a policy similarity embedding (Section 4). We use contrastive methods (Hadsell et al., 2006; Oord et al., 2018) whose track record for representation learning is well-established. We adapt SimCLR (Chen et al., 2020), a popular contrastive method for learning embeddings of image inputs. Given two inputs $x$ and $y$, their embedding similarity is $s_{\theta }(x,y)=sim(z_{\theta }(x),z_{\theta }(y))$, where $sim(u,v)=\frac{u^{T}v}{\Vert u\Vert \Vert v\Vert }$ denotes the cosine similarity function. SimCLR aims to maximize similarity between augmented versions of an image (e.g., crops, colour changes) while minimizing its similarity to other images. The loss used by SimCLR for two versions $x,y$ of an image, and a set ${𝒳}^{\prime }$ containing other images is:

$${\mathrm{\ell }}_{\theta }(x,y;{𝒳}^{\prime })=-\log \frac{\exp (\lambda s_{\theta }(x,y))}{\exp (\lambda s_{\theta }(x,y))+{\sum }_{x^{\prime }\in {𝒳}^{\prime }\setminus \{x\}}\exp (\lambda s_{\theta }(x^{\prime },y))}$$ (2)

where $\lambda$ is an inverse temperature hyperparameter. The overall SimCLR loss is then the expected value of ${\mathrm{\ell }}_{\theta }(x,y;S)$, when $x$, $y$, and ${𝒳}^{\prime }$ are drawn from some augmented training distribution.

3 Policy Similarity Metric

Figure 2: Cyan edges represent actions with a positive reward, which are also the optimal actions. Zero rewards everywhere else. $x_0,y_0$ are the start states while $x_2,y_2$ are the terminal states.

A useful tool in learning a policy that generalizes is to understand which states result in similar behavior, and which do not. To be maximally effective, this similarity should go beyond the immediately chosen action and consider long-term behavior. In this regards, the $\pi$-bisimulation metrics are interesting as they are based on the full sequence of future rewards received from different states. However, considering rewards can be both too restrictive (when the policies are the same, but the obtained rewards are not; see Figure 2) or too permissive (when the policies are different, but the obtained rewards are not; see Figure 5(a)). In fact, $\pi$-bisimulation metrics actually lead to poor generalization in our experiments (Sections 5.1 and 5.2).

To address this issue, we instead consider the similarity between policies themselves. We replace the absolute reward difference by a probability pseudometric between policies, denoted Dist. Additionally, since we would like to perform well in unseen environments, we are interested in similarity in optimal behavior. We thus use ${\pi }^{\ast }$ as the grounding policy. This yields the policy similarity metric (PSM), for which states are close when the optimal policies in these states and future states are similar. For a given Dist, the PSM $d^{\ast }:𝒮\times 𝒮\to ℝ$ satisfies the recursive equation

$d^{\ast }(x,y)=\underset{\text{(A)}}{\underbrace{\text{Dist}({\pi }^{\ast }(x),{\pi }^{\ast }(y))}}+\underset{\text{(B)}}{\underbrace{\gamma {𝒲}_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |x),P^{{\pi }^{\ast }}(\cdot |y))}}.$ (3)

The Dist term captures the difference in local optimal behavior (A) while ${𝒲}_1$ captures long-term optimal behavior difference (B); the exact weights assigned to the two are given by the discount factor. Furthermore, when Dist is bounded, $d^{\ast }$ is guaranteed to be finite. While there are technically multiple PSMs (one for each Dist), we omit this distinction whenever clear from context. A proof of the uniqueness of $d^{\ast }$ is given in Proposition A.1.

Our main use of PSM will be to compare states across environments. In this context, we identify the terms in Equation 3 with specific environments for clarity and write (despite its technical inaccuracy)

$d^{\ast }(x,y)=\text{Dist}({\pi }_{𝒳}^{\ast }(x),{\pi }_{𝒴}^{\ast }(y))+\gamma {𝒲}_1(d^{\ast })(P_{𝒳}^{{\pi }^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y)).$

PSM is applicable to both discrete and continuous action spaces. In our experiments, Dist is the total variation distance ($TV$) when $𝒜$ is discrete, and we use the ${\mathrm{\ell }}_1$ distance between the mean actions of the two policies when $𝒜$ is continuous. PSM can be computed iteratively using dynamic programming (Ferns et al., 2011) (more details in Section C.1). Furthermore, when ${\pi }^{\ast }$ is unavailable on training environments, we replace it by an approximation ${\widehat{\pi }}^{\ast }$ to obtain an approximate PSM, which is close to the exact PSM depending on the suboptimality of ${\widehat{\pi }}^{\ast }$ (Proposition C.3).

Despite resembling $\pi$-bisimulation metrics in form, the PSM possesses different characteristics which are better suited to the problem of generalizing a learned policy. To illustrate this point, consider the following simple nearest-neighbour scheme: Given a state $y\in 𝒴$, denote its closest match in $𝒳$ by ${\tilde{x}}_y:=\arg {\min }_{x\in 𝒳}{d}^{\ast }(x,y)$. Suppose that we use this scheme to transfer ${\pi }_{𝒳}^{\ast }$ to ${ℳ}_{𝒴}$, in the sense that we behave according to the policy $\tilde{\pi }(y)={\pi }_{𝒳}^{\ast }({\tilde{x}}_y)$. We can then bound the difference between $\tilde{\pi }$ and ${\pi }_{𝒴}^{\ast }$, something that is not possible if $d^{\ast }$ is replaced by a $\pi$-bisimulation metric.

Theorem 1.

[Bound on policy transfer] For any $y\in 𝒴$, let $Y_y^t\sim P^{\tilde{\pi }}(\cdot |Y_y^{t-1})$ define the sequence of random states encountered starting in $Y_y^0=y$ and following policy $\tilde{\pi }$. We have:

$${𝔼}_{Y_y^t}\left[\sum _{t\ge 0}{\gamma }^{t}TV(\tilde{\pi }(Y_y^t),{\pi }^{\ast }(Y_y^t))\right]\le \frac{1+\gamma }{1-\gamma }{d}^{\ast }({\tilde{x}}_y,y).$$

The proof is in the Appendix (Section A). Theorem 1 is non-vacuous whenever . In particular, $d^{\ast }({\tilde{x}}_y,y)=0$ implies that the transferred policy is optimal. Although this scheme is not practical (computing $d^{\ast }$ requires knowledge of ${\pi }_{𝒴}^{\ast }$), it shows that meaningful policy generalization can be obtained if we can find a mapping that generalizes across $𝒮$. Put another way, PSM gives us a principled way of lifting generalization across inputs (a supervised learning problem) to generalization across environments. We now describe how to employ PSM for learning representations that put together states in which the agent’s long-term optimal behavior is similar.

4 Learning Contrastive Metric Embeddings

Figure 3: Architecture for learning CMEs. Given an input pair $(x,y)$, we first apply the (optional) data augmentation operator $\mathrm{\Psi }$ to produce the input augmentations ${\mathrm{\Psi }}_x:=\mathrm{\Psi }(x),{\mathrm{\Psi }}_y:=\mathrm{\Psi }(y)$. When not using data augmentation, $\mathrm{\Psi }$ is equal to the identity operator, that is, $\forall x\mathrm{\Psi }(x)=x$. The agent’s policy network then outputs the representations for these augmentations by applying the encoder $f_{\theta }$, that is, $f_x=f_{\theta }({\mathrm{\Psi }}_x),f_y=f_{\theta }({\mathrm{\Psi }}_y)$. These representations are projected using a non-linear projector $h_{\theta }$ to obtain the embedding $z_{\theta }$, that is, $z_{\theta }(x)=h_{\theta }(f_x),z_{\theta }(y)=h_{\theta }(f_y)$. These metric embeddings are trained using the contrastive loss defined in Equation (4). The policy ${\pi }_{\theta }$ is an affine function of the representation, that is, ${\pi }_{\theta }(\cdot |y)=W^T{f}_y+b$, where $W,b$ are learned weights and biases. The entire network is trained end-to-end jointly using the reinforcement learning (or imitation learning) loss in conjunction with the auxiliary contrastive loss.

To generalize a learned policy to new environments, we build on the success of contrastive representations (Section 2). Given a state similarity metric $d$, we develop a general procedure (Algorithm 1) to learn contrastive metric embeddings (CMEs) for $d$. We utilize the metric $d$ for defining the set of positive and negative pairs, as well as assigning importance weights to these pairs in the contrastive loss (Equation $4$).

Algorithm 1 Contrastive Metric Embeddings (CMEs)

1: Given: State embedding $z_{\theta }(\cdot )$, Metric $d(\cdot ,\cdot )$ Training environments ${\{{ℳ}_i\}}_{i=1}^N$. Hyperparameters: Temperature $1/\lambda$, Scale $\beta$, Total training steps $K$

2: for each step $k=1,\mathrm{\dots },K$ do

3: Sample a pair of training MDPs ${ℳ}_{𝒳},{ℳ}_{𝒴}$

4: Update $\theta$ to minimize ${ℒ}_{\mathrm{CME}}$ where ${ℒ}_{\mathrm{CME}}={𝔼}_{{ℳ}_{𝒳},{ℳ}_{𝒴}\sim \rho }\left[L_{\theta }({ℳ}_{𝒳},{ℳ}_{𝒴})\right]$

5: end for

We first apply a transformation to convert $d$ to a similarity measure $\mathrm{\Gamma }$, bounded in [0, 1] for “soft” similarities. In this work, we transform $d$ using the Gaussian kernel with a positive scale parameter $\beta$, that is, $\mathrm{\Gamma }(x,y)=\exp (-d(x,y)/\beta )$. $\beta$ controls the sensitivity of the similarity measure to $d$.

Second, we select the positive and negative pairs given a set of states ${𝒳}^{\prime }\subseteq 𝒳$ and $𝒴$ from MDPs ${ℳ}_{𝒳},{ℳ}_{𝒴}$ respectively. For each anchor state $y\in 𝒴$, we use its nearest neighbor in ${𝒳}^{\prime }$ based on the similarity measure $\mathrm{\Gamma }$ to define the positive pairs $\{({\tilde{x}}_y,y)\}$, where ${\tilde{x}}_y=\underset{x\in {𝒳}^{\prime }}{\mathrm{argmax}}\mathrm{\Gamma }(x,y)$. The remaining states in ${𝒳}^{\prime }$, paired with $y$, are used as negative pairs. This choice of pairs is motivated by Theorem 1, which shows that if we transfer the optimal policy in ${ℳ}_{𝒳}$ to the nearest neighbors defined using PSM, its performance in ${ℳ}_{𝒴}$ has suboptimality bounded by PSM.

Next, we define a soft version of the SimCLR contrastive loss (Equation 2) for learning the function $z_{\theta }$, which maps states (usually high-dimensional) to embeddings. Given a positive state pair $({\tilde{x}}_y,y)$, the set ${𝒳}^{\prime }$, and the similarity measure $\mathrm{\Gamma }$, the loss (pseudocode provided in Section I.2) is given by

${\mathrm{\ell }}_{\theta }({\tilde{x}}_y,y;{𝒳}^{\prime })=-\log {\displaystyle \frac{\mathrm{\Gamma }({\tilde{x}}_y,y)\exp (\lambda s_{\theta }({\tilde{x}}_y,y))}{\mathrm{\Gamma }({\tilde{x}}_y,y)\exp (\lambda s_{\theta }({\tilde{x}}_y,y))+{\sum }_{x^{\prime }\in {𝒳}^{\prime }\setminus \{{\tilde{x}}_y\}}(1-\mathrm{\Gamma }(x^{\prime },y))\exp (\lambda s_{\theta }(x^{\prime },y))}}(4),$

where we use the same notation as in Equation 2. Following SimCLR, we use a non-linear projection of the representation as $z_{\theta }$ (Figure 3). The agent’s policy is an affine function of the representation.

The total contrastive loss for ${ℳ}_{𝒳}$ and ${ℳ}_{𝒴}$ utilizes the optimal trajectories ${\tau }_{𝒳}^{\ast }={\{x_t\}}_{t=1}^N$ and ${\tau }_{𝒴}^{\ast }={\{y_t\}}_{t=1}^N$, where $x_{t+1}\sim P_{𝒳}^{{\pi }^{\ast }}(\cdot |x_t)$ and $y_{t+1}\sim P_{𝒴}^{{\pi }^{\ast }}(\cdot |y_t)$. We set ${𝒳}^{\prime }={\tau }_{𝒳}^{\ast }$ and define

$$L_{\theta }({ℳ}_{𝒳},{ℳ}_{𝒴})={𝔼}_{y\sim {\tau }_{𝒴}^{\ast }}\left[{\mathrm{\ell }}_{\theta }({\tilde{x}}_y,y;{\tau }_{𝒳}^{\ast })\right]\mathrm{where}{\tilde{x}}_y=\underset{x\in {\tau }_{𝒳}^{\ast }}{\mathrm{argmax}}\mathrm{\Gamma }(x,y).$$

We refer to CMEs learned with policy similarity metric as policy similarity embeddings (PSEs). PSEs can trivially be combined with data augmentation by using augmented states when computing state embeddings. We simultaneously learn PSEs with the RL agent by adding ${ℒ}_{\mathrm{CME}}$ (Algorithm 1) as an auxiliary objective during training. Next, we illustrate the benefits from this auxiliary objective.

5 Jumping Task from Pixels: A Case Study

Task Description. The jumping task (Tachet des Combes et al., 2018) (Figure 1) captures, using well-defined factors of variations, whether agents can learn the correct invariances required for generalization, directly from image inputs. The task consists of an agent trying to jump over an obstacle. The agent has access to two actions: right and jump. The agent needs to time the jump precisely, at a specific distance from the obstacle, otherwise it will eventually hit the obstacle. Different tasks consist in shifting the floor height and/or the obstacle position. To generalize, the agent needs to be invariant to the floor height while jump based on the obstacle position. The obstacle can be in 26 different locations while the floor has 11 different heights, totaling 286 tasks.

Problem Setup. We split the problem into 18 seen (training) and 268 unseen (test) tasks to stress test generalization using a few changes in the underlying factors of variations seen during training. The small number of positive examples results in a highly unbalanced classification problem with low amounts of data, making it challenging without additional inductive biases. Thus, we evaluate generalization in regimes with and without data augmentation. The different grids configurations (Figure 4) capture different types of generalization: the “wide” grid tests generalization via “interpolation”, the “narrow” grid tests out-of-distribution generalization via “extrapolation”, and the random grid instances evaluate generalization similar to supervised learning where train and test samples are drawn i.i.d. from the same distribution.

We used RandConv (Lee et al., 2020a), a state-of-the-art data augmentation for generalization. For hyperparameter selection, we evaluate all agents on a validation set containing 54 unseen tasks in the “wide” grid (Figure 4(a)) and pick the parameters with the best validation performance. We use these fixed parameters for all grid configurations to show the robustness of PSEs to hyperparameter tuning. We first compute the optimal trajectories in the training tasks. Using these trajectories, we compute PSM using dynamic programming (Section C.1). We train the agent by imitation learning, combined with an auxiliary loss for PSEs (Section 4). More details are in Section F.

Table 1: Percentage (%) of test tasks solved by different methods without and with data augmentation. The “wide”, “narrow”, and random grids are described in Figure 4. We report average performance across 100 runs with different random initializations, with standard deviation between parentheses.

Data Augmentation	Method	Grid Configuration (%)
		“Wide”	“Narrow”	Random
✗	Dropout and ${\mathrm{\ell }}_2$ reg.	17.8 (2.2)	10.2 (4.6)	9.3 (5.4)
	Bisimulation Transfer4	17.9 (0.0)	17.9 (0.0)	30.9 (4.2)
	PSEs	33.6 (10.0)	9.3 (5.3)	37.7 (10.4)
✓	RandConv	50.7 (24.2)	33.7 (11.8)	71.3 (15.6)
	RandConv + ${\pi }^{\ast }$-bisimulation	41.4 (17.6)	17.4 (6.7)	33.4 (15.6)
	RandConv + PSEs	87.0 (10.1)	52.4 (5.8)	83.4 (10.1)

(a) “Wide” grid

(b) “Narrow” grid

(c) Random grid

Figure 4: Jumping Task: Visualization of average performance of PSEs with data augmentation across different configurations. We plot the median performance across 100 runs. Each tile in the grid represents a different task (obstacle position/floor height combination). For each grid configuration, the height varies along the $y$-axis (11 heights) while the obstacle position varies along the $x$-axis (26 locations). The red letter $\top$ indicates the training tasks. Random grid depicts only one instance, each run consisted of a different test/train split. Beige tiles are tasks PSEs solved while black tiles are tasks PSEs did not solve when used with data augmentation. These results were chosen across all the 100 runs to demonstrate what the average reported performance looks like.

5.1 Evaluating Generalization on Jumping Task

We show the efficacy of PSEs compared to common generalization approaches such as regularization (e.g., Cobbe et al., 2019; Farebrother et al., 2018), and data augmentation (e.g., Lee et al., 2020a; Laskin et al., 2020a), which are quite effective on pixel-based RL tasks. We also contrast PSEs with bisimulation transfer (Castro & Precup, 2010), a tabular state-based transfer approach based on bisimulation metrics which does not do any learning and bisimulation preserving representations (Zhang et al., 2021), showing the advantages of PSM over a prevalent state similarity metric.

We first investigated how well PSEs generalize over existing methods without incorporating additional domain knowledge during training. Table 1 summarizes, in the setting without data augmentation, the performance of these methods in different train/test splits (c.f. Figure 4 for a detailed description). PSEs, with only 18 examples, already leads to better performance than standard regularization.

PSEs also outperform bisimulation transfer in the “wide” and random grids. Although bisimulation transfer is impractical when evaluating zero-shot generalization, we still performed this comparison, unfair to PSEs, to highlight their efficacy. PSEs perform better because, in contrast to bisimulation, PSM is reward agnostic (c.f. Proposition C.1) – the expected return of the jump action is quite different depending on the obstacle position (c.f. Figure F.2 for a visual juxtaposition of PSM and bisimulation). Overall, these results are promising because they place PSEs as an effective generalization method that does not rely on data augmentation.

Nevertheless, PSEs are complementary to data augmentation, which consistently improves generalization in deep RL. We compared RandConv combined with PSEs to simply using RandConv. Domain-specific augmentation also succeeds in the jumping task. Thus, it is not surprising that RandConv is so effective compared to techniques without augmentation. Table 1 ($2^{nd}$ row) shows that PSEs substantially improve the performance of RandConv across all grid configurations. Moreover, Table 1 ($2^{nd}$ row) illustrates that when combined with RandConv, bisimulation preserving representations (Zhang et al., 2021) diminish generalization by $30-50\%$ relative to PSEs.

Notably, Table 1 ($1^{st}$ row) indicates that learning-based methods are ineffective on the “narrow” grid without data augmentation. That said, PSEs do work quite well when combined with RandConv. However, even with data augmentation, generalization in “narrow” grid happens only around the vicinity of training tasks, exhibiting the challenge this grid poses for learning-based methods. We believe this is due to the poor extrapolation ability of neural networks (e.g., Haley & Soloway, 1992; Xu et al., 2020), which is more perceptible without prior inductive bias from data augmentation.

5.2 Understanding gains from PSEs: Ablations and Visualizations

(a) Optimal Trajectories

(b) PSEs

(c) PSM + ${\mathrm{\ell }}_2$ embeddings

(d) ${\pi }^{\ast }$-bisim. + CMEs

Figure 5: Embedding visualization. (a) Optimal trajectories on original jumping task (visualized as coloured blocks) with different obstacle positions. We visualize the hidden representations using UMAP, where the color of points indicate the tasks of the corresponding observations. Points with the same number label correspond to same distance of the agent from the obstacle, the underlying optimal invariant feature across tasks.

Table 2: Ablating PSEs. Percentage (%) of test tasks solved when we ablate the similarity metric and learning procedure for metric embeddings in the data augmentation setting on wide grid. PSEs, which combine CMEs with PSM, considerably outperform other embeddings. We report the average performance across 100 runs with standard deviation between parentheses. All ablations except PSEs deteriorate performance compared to just using data augmentation (RandConv), as reported in Table 1.

Metric / Embedding	${\mathrm{\ell }}_2$-embeddings	CMEs
${\pi }^{\ast }$-bisimulation	41.4 (17.6)	23.1 (7.6)
PSM	17.5 (8.4)	87.0 (10.1)

PSEs are contrastive metric embeddings (CMEs) learned with PSM. We investigate the gains from CMEs (Section 4) and PSM (Section 3) by ablating them. CMEs can be learned with any state similarity metric – we use ${\pi }^{\ast }$-bisimulation (Castro, 2020) as an alternative. Similarly, PSM can be used with any metric embedding – we use ${\mathrm{\ell }}_2$-embeddings (Section D) as an alternative, which Zhang et al. (2021) employed with ${\pi }^{\ast }$-bisimulation for learning representations in a single task RL setting. For a fair comparison, we tune hyperparameters for 128 trials for each ablation entry in Table 2.

Table 2 show that PSEs ($=$ PSM + CMEs) generalize significantly better than ${\pi }^{\ast }$-bisimulation with CMEs or ${\mathrm{\ell }}_2$-embeddings, both of which significantly degrade performance ($-60\%$ and $-45\%$, respectively). This is expected, as ${\pi }^{\ast }$-bisimulation imposes the incorrect invariances for the jumping task (c.f. Figures 2(a) and 2(d)). Additionally, looking at the rows of Table 2, CMEs are superior than ${\mathrm{\ell }}_2$-embeddings for PSM while inferior for ${\pi }^{\ast }$-bisimulation. This result is in line with the hypothesis that CMEs better enforce the invariances encoded by a similarity metric compared to ${\mathrm{\ell }}_2$-embeddings (c.f. Figures 5(b) and 5(c)).

Visualizing learned representations. We visualize the metric embeddings in the ablation above by projecting them to two dimensions with UMAP (McInnes et al., 2018), a popular visualization technique for high dimensional data which better preserves the data’s global structure compared to other methods such as t-SNE (Coenen & Pearce, 2019).

Figure 5 shows that PSEs partition the states into two sets: (1) states where a single suboptimal action leads to failure (all states before jump) and (2) states where actions do not affect the final outcome (states after jump). Additionally, PSEs align the labeled states in the first set, which have a PSM distance of zero. These aligned states have the same distance from the obstacle, the invariant feature that generalizes across tasks. On the other hand, ${\mathrm{\ell }}_2$-embeddings (Zhang et al., 2021) with PSM do not align states with zero PSM except the state with the jump action – poor generalization, as observed empirically, is likely due to states with the same optimal behavior ending up with distant embeddings. CMEs with ${\pi }^{\ast }$-bisimulation align states with ${\pi }^{\ast }$-bisimulation distance of zero – such states are equidistant from the start state and have different optimal behavior for any pair of tasks with different obstacle positions (Figure 2(c)).

5.3 Effect of Policy Suboptimality on PSEs

To understand the sensitivity of learning effective PSEs to the quality of the policies, we compute PSEs using $ϵ$-suboptimal policies on the jumping task, which take the optimal action with probability $1-ϵ$ and the subopotimal action with probability $ϵ$.

Figure 6: Percentage (%) of test tasks solved using PSEs computed using $ϵ$-suboptimal policies on the “wide” grid. We report the mean across 100 runs. Error bars show one standard deviation.

We evaluate the generalization performance of PSEs for increasingly suboptimal policies, ranging from the optimal policy ($ϵ=0$) to the uniform random policy ($ϵ=0.5$). To isolate the effect of suboptimality on PSEs, the agent still imitates the optimal actions during training for all $ϵ$.

Figure 6 shows that PSEs show near-optimal generalization with $ϵ\le 0.4$ while degrade generalization with an uniform random policy. This result is well-aligned with Proposition C.3, which shows that for policies with decreasing suboptimality, the PSM approximation becomes more accurate, resulting in improved PSEs. Overall, this study confirms that the utility of PSEs for generalization is robust to suboptimality. One reason for this robustness is that PSEs are likely to align states with similar long-term greedy optimal actions, resulting in good performance even with suboptimal policies that preserve these greedy actions.

5.4 Jumping Task with Colors: Where Task-dependent Invariance Matters

Figure 7: Percentage (%) of red obstacle test tasks solved when trained, jointly with red and green obstacles, on the “wide” grid. We report the mean across 100 runs. Error bars show 99% confidence interval for the mean.

The task-dependent invariances captured by PSEs are usually orthogonal to task-agnostic invariances from data augmentation. This difference is important because, for certain RL tasks, data augmentation can erroneously alias states with different optimal behavior. Domain knowledge is often required to select appropriate augmentations, otherwise augmentations can even hurt generalization. In contrast, PSEs do not require any domain knowledge but instead exploit the inherent structure of the RL tasks.

To demonstrate the difference between PSEs and data augmentation, we simply include colored obstacles in the jumping task (see Figure F.5). In this modified task, the optimal behavior of the agent depends on the obstacle color: the agent needs to jump over the red obstacle but strike the green obstacle to get a high return. The red obstacle task has the same difficulty as the original jumping task while the green obstacle task is easier. We jointly train the agent with 18 training tasks each, for both obstacle colors, on the “wide” grid and evaluate generalization on unseen red tasks.

Figure 7 shows the large performance gap between PSEs and data augmentation with RandConv. All methods solve the green obstacle tasks (Table F.1). As opposed to the original jumping task (c.f. Table 1), data augmentation inhibits generalization since RandConv forces the agent to ignore color, conflating the red and green tasks (Figure F.6). PSEs still outperform regularization and data augmentation. Furthermore, data augmentation performs better when combined with PSEs. Thus, PSEs are effective even when data augmentation hurts performance.

6 Additional Empirical Evaluation

In this section, we exhibit that PSM ignores spurious information for generalization using a LQR task (Song et al., 2019) with non-image inputs. Then, we demonstrate the scalability of PSEs without explicit access to optimal policies in an RL setting with continuous actions, using Distracting DM Control Suite (Stone et al., 2021).

6.1 LQR with Spurious Correlations

We show how representations learned using PSM, when faced with semantically equivalent environments, can learn the main factors of variation and ignore spurious correlations that hinder

Figure 8: LQR generalization: Absolute test error in LQR cost relative to the oracle (which has access to true state), of various methods trained with $n_d$ distractors on 2 training environments. Lower error is better. PSM achieves near-optimal performance. We report mean across 100 different seeds. Error bars (non-existent for most methods) show 99% confidence interval for mean.

generalization. We use LQR with distractors (Song et al., 2019; Sonar et al., 2020) to assess generalization in a feature-based RL setting with linear function approximation. The distractors are input features that are spuriously correlated with optimal actions and can be used for predicting these actions during training, but hurt generalization. The agent learns a linear policy using 2 environments with fixed distractors. This policy is evaluated on environments with unseen distractors.

We aggregate state pairs in training environments with near-zero PSM. We contrast this approach with (i) IPO (Sonar et al., 2020), a policy optimization method based on IRM (Arjovsky et al., 2019) for out-of-distribution generalization, (ii) overparametrization, which leads to better generalization via implicit regularization (Song et al., 2019), and (iii) weight sparsity using ${\mathrm{\ell }}_1$-regularization since the policy weights which generalize in this task are sparse.

All methods optimally solve the training environments; however, the baselines perform abysmally in terms of generalization compared to state aggregation with PSM (Figure 8), indicating their reliance on distractors. PSM obtains near-optimal generalization which we corroborate through this conjecture (Section G.1): Assuming zero state aggregation error with PSM, the policy learned using gradient descent is independent of the distractors. Refer to Section G for a detailed discussion.

6.2 Distracting DM Control Suite

Figure 9: Distracting Control Suite: Snapshots of training environments.

Finally, we demonstrate scalability of PSEs on the Distracting DM Control Suite (DCS) (Stone et al., 2021), which tests whether agents can ignore high-dimensional visual distractors irrelevant to the RL task. Since we do not have access to optimal training policies, we use learned policies as proxy for ${\pi }^{\ast }$ for computing PSM as well as collecting data for optimizing PSEs. Even with this approximation, PSEs outperform state-of-the-art data augmentation.

DCS extends DM Control (Tassa et al., 2020) with visual distractions. We use the dynamic background distractions (Stone et al., 2021; Zhang et al., 2018b) where a video is played in the background from a specific frame. The video and the frame are randomly sampled every new episode. We use 2 videos during training (Figure 9) and evaluate generalization on 30 unseen videos (Figure H.1).

All agents are built on top of SAC (Haarnoja et al., 2018) combined with DrQ (Kostrikov et al., 2020), an augmentation method with state-of-the-art performance on DM control. Without data augmentation on DM control, SAC performs poorly, even during training (Kostrikov et al., 2020). We augment DrQ with an auxiliary loss for learning PSEs and compare it with DrQ (Table 3). Orthogonal to DrQ, PSEs align representations of different states across environments based on PSM (c.f. Figure 3). All agents are trained for 500K environment steps with random crop augmentation. For computing PSM, we use policies learned by DrQ pretrained on training environments for 500K steps.

Table 3: Generalization performance with unseen distractions in the Distracting Suite at 500K steps. We report the average scores across 5 seeds $\pm$ standard error. All methods are added to SAC (Haarnoja et al., 2018). Pretrained initialization uses DrQ trained for 500K steps. Figures H.2 and H.3 show learning curves.

Initialization	Method	BiC-catch	C-swingup	C-run	F-spin	R-easy	W-walk
Random	DrQ	747$\pm$28	582$\pm$42	220$\pm$12	646$\pm$54	931$\pm$14	549$\pm$83
	DrQ + PSEs	821$\pm$17	749$\pm$19	308$\pm$12	779$\pm$49	955$\pm$10	789$\pm$28
Pretrained	DrQ	748$\pm$30	689$\pm$22	219$\pm$10	764$\pm$48	943$\pm$10	709$\pm$29
	DrQ + PSEs	805$\pm$25	753$\pm$13	282$\pm$8	803$\pm$19	962$\pm$11	829$\pm$21

First, we investigate how much better PSEs generalize relative to DrQ, assuming the agent is provided with PSM beforehand. The agent’s policy is randomly initialized so that additional gains over DrQ can be attributed to the auxiliary information from PSM. The substantial gains in Table 3 indicate that PSEs are more effective than DrQ for encoding invariance to distractors.

Since PSEs utilize PSM approximated using pretrained policies, we also compare to a DrQ agent where we initialize it using these pretrained policies. This comparison provides the same auxiliary information to DrQ as available to PSEs, thus, the generalization difference stems from how they utilize this information. Table 3 demonstrates that PSEs outperform DrQ with pretrained initialization, indicating that the additional pretraining steps are more judiciously utilized for computing PSM as opposed to just longer training with DrQ. More details, including learning curves, are in Section H.

7 Related Work

PSM (Section 3) is inspired by bisimulation metrics (Section B). However, different than traditional bisimulation (e.g., Larsen & Skou, 1991; Givan et al., 2003; Ferns et al., 2011), PSM is more tractable as it defined with respect to a single policy similar to the recently proposed ${\pi }^{\ast }$-bisimulation (Castro, 2020; Zhang et al., 2021). However, in contrast to PSM, bisimulation metrics rely on reward information and may not provide a meaningful notion of behavioral similarity in certain environments (Section 5). For example, states similar under PSM would have similar optimal policies, yet can have arbitrarily large ${\pi }^{\ast }$-bisimulation distance between them (Proposition C.1).

PSEs (Section 4) use contrastive learning to encode behavior similarity (Section 3) across MDPs. Previously, contrastive learning has been applied for imposing state self-consistency (Laskin et al., 2020b), capturing predictive information (Oord et al., 2018; Mazoure et al., 2020; Lee et al., 2020b) or encoding transition dynamics (van der Pol et al., 2020; Stooke et al., 2020; Schwarzer et al., 2020) within an MDP. These methods can be integrated with PSEs to encode additional invariances. Interestingly, in a similar spirit to PSEs, Pacchiano et al. (2020); Moskovitz et al. (2021) explore comparing behavioral similarity between policies to guide policy optimization within an MDP.

PSEs are complementary to data augmentation methods (Kostrikov et al., 2020; Lee et al., 2020a; Raileanu et al., 2020; Ye et al., 2020), which have recently been shown to significantly improve agents’ generalization capabilities. In fact, we combine PSEs to state-of-the-art augmentation methods including random convolutions (Lee et al., 2020a; Laskin et al., 2020a) in the jumping task and DrQ (Kostrikov et al., 2020) on Distracting Control Suite, leading to performance improvement. Furthermore, for certain RL tasks, it can be unclear what an optimality invariant augmentation would look like (Section 5.4). PSM can quantify the invariance of such augmentations (Proposition C.2).

8 Conclusion

This paper advances generalization in RL by two contributions: (1) the policy similarity metric (PSM) which provides a new notion of state similarity based on behavior proximity, and (2) contrastive metric embeddings, which harness the benefits of contrastive learning for representations based on a similarity metric. PSEs combine these two ideas to improve generalization. Overall, this paper shows the benefits of exploiting the inherent structure in RL for learning effective representations.

Appendix

Appendix A Proofs

We begin by defining some notation which will be used throughout these results:

• We denote ${𝔼}_{t\ge 0}[{\gamma }^{t}TV(\tilde{\pi }(Y_y^t),{\pi }^{\ast }(Y_y^t))]={𝔼}_{Y_y^t}\left[{\sum }_{t\ge 0}{\gamma }^{t}TV(\tilde{\pi }(Y_y^t),{\pi }^{\ast }(Y_y^t))\right]$

• For any $y\in Y$, let $Y_y^t\sim P^{\tilde{\pi }}(\cdot |Y_y^{t-1})$, where $Y_y^0=y$.

• .

We now proceed with some technical lemmas necessary for the main result.

Lemma 1.

Given any two pseudometrics $d,d^{\prime }\in 𝕄$ and probability distributions $P_{𝒳},P_{𝒴}$ where $𝒳,𝒴\subset 𝒮$, we have:

$${𝒲}^1(d)(P_{𝒳},P_{𝒴})\le \Vert d-d^{\prime }\Vert +{𝒲}^1(d^{\prime })(P_{𝒳},P_{𝒴})$$

Proof.

Note that the dual of the linear program for computing ${𝒲}^1(d)(P_{𝒳},P_{𝒴})$ is given by

$\underset{\mathrm{\Gamma }}{\min }{\displaystyle \sum _{x\in 𝒳,y\in 𝒴}}\mathrm{\Gamma }(x,y)d(x,y)$

$\mathrm{subject}\mathrm{to}{\displaystyle \sum _x}\mathrm{\Gamma }(x,y)=P_{𝒴}(y)\forall y,{\displaystyle \sum _y}\mathrm{\Gamma }(x,y)=P_{𝒳}(x)\forall x,\mathrm{\Gamma }(x,y)\ge 0\forall x,y$

Using the dual formulation subject to the constraints above, ${𝒲}^1(d)$ can be written as

${𝒲}^1(d)(P_{𝒳},P_{𝒴})\le \Vert d-d^{\prime }\Vert$ $={𝒲}^1(d-d^{\prime }+d^{\prime })(P_{𝒳},P_{𝒴})\le \Vert d-d^{\prime }\Vert$

$=\underset{\mathrm{\Gamma }}{\min }{\displaystyle \sum _{x\in 𝒳,y\in 𝒴}}\mathrm{\Gamma }(x,y)(d(x,y)-d^{\prime }(x,y)+d^{\prime }(x,y))$

$\le \underset{\mathrm{\Gamma }}{\min }{\displaystyle \sum _{x\in 𝒳,y\in 𝒴}}\mathrm{\Gamma }(x,y)(\Vert d-d^{\prime }\Vert +d^{\prime }(x,y))$

$=\Vert d-d^{\prime }\Vert +\underset{\mathrm{\Gamma }}{\min }{\displaystyle \sum _{x\in 𝒳,y\in 𝒴}}\mathrm{\Gamma }(x,y)d^{\prime }(x,y)$

$=\Vert d-d^{\prime }\Vert +{𝒲}^1(d^{\prime })(P_{𝒳},P_{𝒴})$

∎

Lemma 2.

Given any $y_0\in Y$, we have:

${\displaystyle \sum _{y_1\in Y}}\left(P^{\tilde{\pi }}(y_1|y_0)-P^{{\pi }^{\ast }}(y_1|y_0)\right)T{V}^n(Y_{y_1}^0)\le {\displaystyle \frac{2}{1-\gamma }}TV(\tilde{\pi }(y_0),{\pi }^{\ast }(y_0))$

Proof.

${\displaystyle \sum _{y_1\in Y}}\left(P^{\tilde{\pi }}(y_1|y_0)-P^{{\pi }^{\ast }}(y_1|y_0)\right)T{V}^n(Y_{y_1}^0)$ $\le \left|{\displaystyle \sum _{y_1\in Y}}(P^{\tilde{\pi }}(y_1|y_0)-P^{{\pi }^{\ast }}(y_1|y_0))T{V}^n(Y_{y_1}^0)\right|$

$\le {\displaystyle \sum _{y_1\in Y}}\left|{\displaystyle \sum _{a\in A}}P(y_1|y_0,a)(\tilde{\pi }(a|y_0)-{\pi }^{\ast }(a|y_0))\right|T{V}^n(Y_{y_1}^0)$

$\le {\displaystyle \frac{1}{1-\gamma }}{\displaystyle \sum _{y_1\in Y}}{\displaystyle \sum _{a\in A}}P(y_1|y_0,a)|\tilde{\pi }(a|y_0)-{\pi }^{\ast }(a|y_0)|$

$={\displaystyle \frac{1}{1-\gamma }}{\displaystyle \sum _{a\in A}}|\tilde{\pi }(a|y_0)-{\pi }^{\ast }(a|y_0)|{\displaystyle \sum _{y_1\in Y}}P(y_1|y_0,a)$

$={\displaystyle \frac{1}{1-\gamma }}{\displaystyle \sum _{a\in A}}|\tilde{\pi }(a|y_0)-{\pi }^{\ast }(a|y_0)|$

$={\displaystyle \frac{2}{1-\gamma }}TV(\tilde{\pi }(y_0),{\pi }^{\ast }(y_0))$

∎

Lemma 3.

Given any $y_0\in Y$, if $T{V}^n(Y_{y_1}^0)\le \frac{1+\gamma }{1-\gamma }{d}^{\ast }({\tilde{x}}_{y_1},y_1)$, we have:

$$\sum _{y_1\in Y}{P}^{{\pi }^{\ast }}(y_1|y_0)T{V}^n(Y_{y_1}^0)\le \frac{1+\gamma }{1-\gamma }{W}_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |{\tilde{x}}_{y_0}),P^{{\pi }^{\ast }}(\cdot |y_0))$$

Proof.

Note that we have the following equality, where $\mathrm{𝟎}$ is a vector of zeros:

$$\sum _{y_1\in Y}{P}^{{\pi }^{\ast }}(y_1|y_0)T{V}^n(Y_{y_1}^0)=\sum _{y_1\in Y}{P}^{{\pi }^{\ast }}(y_1|y_0)T{V}^n(Y_{y_1}^0)-\sum _{x\in X}{P}^{{\pi }^{\ast }}(x|{\tilde{x}}_{y_0})\mathrm{𝟎}$$

which is the same form as the primal LP for $W_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |y_0),P^{{\pi }^{\ast }}(\cdot |{\tilde{x}}_{y_0}))$. By assumption, we have that

$$T{V}^n(Y_{y_1}^0)\le \frac{1+\gamma }{1-\gamma }{d}^{\ast }({\tilde{x}}_{y_1},y_1)$$

This implies that $\frac{1-\gamma }{1+\gamma }T{V}^n(Y_{\cdot }^0)$ is a feasible solution to $W_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |y_0),P^{{\pi }^{\ast }}(\cdot |{\tilde{x}}_{y_0}))$:

$$\sum _{y_1\in Y}{P}^{{\pi }^{\ast }}(y_1|y_0)\frac{1-\gamma }{1+\gamma }T{V}^n(Y_{y_1}^0)\le W_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |{\tilde{x}}_{y_0}),P^{{\pi }^{\ast }}(\cdot |y_0))$$

and the result follows. ∎

Proposition A.1.

The operator $ℱ$ given by:

$$ℱ(d)(x,y)=\text{Dist}({\pi }^{\ast }(x),{\pi }^{\ast }(y))+\gamma {𝒲}_1(d)(P_{𝒳}^{{\pi }^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y))$$

is a contraction mapping and has a unique fixed point for a bounded $dist$.

Proof.

We first prove that $ℱ$ is contraction mapping. Then, a simple application of the Banach Fixed Point Theorem asserts that $ℱ$ has a unique fixed point. Note that for all pseudometrics $d,d^{\prime }\in 𝕄$, and for all states $x\in 𝒳$, $y\in 𝒴$,

$ℱ(d)(x,y)-ℱ(d^{\prime })(x,y)$

$=\gamma ({𝒲}_1(d)(P^{{\pi }_{𝒳}^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y))-{𝒲}_1(d^{\prime })(P^{{\pi }_{𝒳}^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y)))$

$\stackrel{\text{<a class="ltx\_ref ltx\_refmacro\_autoref ltx\_font\_italic" href="\#Thmlemma1" title="Lemma 1. ‣ Appendix A Proofs ‣ Appendix ‣ Contrastive Behavioral Similarity Embeddings for Generalization in Reinforcement Learning"><span class="ltx\_text ltx\_ref\_tag">1</span></a>}}{\le }\gamma (\parallel d-d^{\prime }\parallel +{𝒲}_1(d^{\prime })(P^{{\pi }_{𝒳}^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y))-{𝒲}_1(d^{\prime })(P^{{\pi }_{𝒳}^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y)))$

$=\gamma \Vert d-d^{\prime }\Vert$

Thus, $\Vert ℱ(d)-ℱ(d^{\prime })\Vert \le \gamma \Vert d-d^{\prime }\Vert$, so that $ℱ$ is a contracting mapping for and has an unique fixed point $d^{\ast }$. ∎

See 1

Proof.

We will prove this by induction. Assuming that the bound holds for $T{V}^n$, we prove the bound holds for $T{V}^{n+1}$. The base case for $n=1$ follows from $TV(\tilde{\pi }(y),{\pi }^{\ast }(y))=TV({\pi }^{\ast }({\tilde{x}}_y),{\pi }^{\ast }(y))\le d^{\ast }({\tilde{x}}_y,y)$. Note that $T{V}^n\le \frac{1-{\gamma }^{n+1}}{1-\gamma }$ since the $TV$ distance per time-step can be at most 1.

Let $P_t^{\pi }(y^{\prime }|y)$ denote the probability of ending in state $y^{\prime }\in Y$ after $t$ steps when following policy $\pi$ and starting from state $y$. We then have:

$T$ $V^{n+1}(Y_y^k)={\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y)TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))+\gamma T{V}^n(Y_y^{k+1})$

$={\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y)\left[TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))+\gamma {\displaystyle \sum _{y_{k+1}\in Y}}{P}^{\tilde{\pi }}(y_{k+1}|y_k)T{V}^n(Y_{y_{k+1}}^0)\right]$

$={\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y)[TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))+\gamma {\displaystyle \sum _{y_{k+1}\in Y}}(P^{\tilde{\pi }}(y_{k+1}|y_k)-P^{{\pi }^{\ast }}(y_{k+1}|y_k))T{V}^n(Y_{y_{k+1}}^0)$

$+\gamma {\displaystyle \sum _{y_{k+1}\in Y}}{P}^{{\pi }^{\ast }}(y_{k+1}|y_k)T{V}^n(Y_{y_{k+1}}^0)]$

$\stackrel{\text{<a class="ltx\_ref ltx\_refmacro\_autoref ltx\_font\_italic" href="\#Thmlemma2" title="Lemma 2. ‣ Appendix A Proofs ‣ Appendix ‣ Contrastive Behavioral Similarity Embeddings for Generalization in Reinforcement Learning"><span class="ltx\_text ltx\_ref\_tag">2</span></a>}}{\le }{\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y)\left[TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))+{\displaystyle \frac{2\gamma }{1-\gamma }}TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))+\gamma {\displaystyle \sum _{y_{k+1}\in Y}}{P}^{{\pi }^{\ast }}(y_{k+1}|y_k)T{V}^n(Y_{y_{k+1}}^0)\right]$

$\stackrel{\text{<a class="ltx\_ref ltx\_refmacro\_autoref ltx\_font\_italic" href="\#Thmlemma3" title="Lemma 3. ‣ Appendix A Proofs ‣ Appendix ‣ Contrastive Behavioral Similarity Embeddings for Generalization in Reinforcement Learning"><span class="ltx\_text ltx\_ref\_tag">3</span></a>}}{\le }{\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y)[TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))$

$+\gamma ({\displaystyle \frac{2}{1-\gamma }}TV(\tilde{\pi }(y_k),{\pi }^{\ast }(y_k))+{\displaystyle \frac{1+\gamma }{1-\gamma }}{W}_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |{\tilde{x}}_{y_k}),P^{{\pi }^{\ast }}(\cdot |y_k)))]$

$={\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y)\left[{\displaystyle \frac{1+\gamma }{1-\gamma }}(TV({\pi }^{\ast }({\tilde{x}}_{y_k}),{\pi }^{\ast }(y_k))+\gamma W_1(d^{\ast })(P^{{\pi }^{\ast }}(\cdot |{\tilde{x}}_{y_k}),P^{{\pi }^{\ast }}(\cdot |y_k)))\right]$

$\le {\displaystyle \sum _{y_k\in Y}}{P}_k^{\tilde{\pi }}(y_k|y){\displaystyle \frac{1+\gamma }{1-\gamma }}{d}^{\ast }({\tilde{x}}_{y_k},y_k)$

Thus, by induction, it follows that for all $n$:

$$T{V}^n(Y_y^0)\le \frac{1+\gamma }{1-\gamma }{d}^{\ast }({\tilde{x}}_y,y),$$

which completes the proof. ∎

Appendix B Bisimulation metrics

Notation. We use the notation as defined in Section 2.

Bisimulation metrics (Givan et al., 2003; Ferns et al., 2011) define a pseudometric $d_{\sim }:𝒮\times 𝒮\to ℝ$ where $d_{\sim }(x,y)$ is defined in terms of distance between immediate rewards and next state distributions. Define ${ℱ}_{\sim }^e:𝕄\to 𝕄$ by

${ℱ}_{\sim }^e(d)(x,y)=\underset{a\in 𝒜}{\max }|R(x,a)-R(y,a)|+\gamma {𝒲}_1(d)(P^a(\cdot |x),P^a(\cdot |y))$ (B.1)

then, ${ℱ}_{\sim }^e$ has a unique fixed point $d_{\sim }$ which is a bisimulation metric. ${ℱ}_{\sim }^e$ uses the 1-Wasserstein metric ${𝒲}_1:𝕄\to {𝕄}_p$. The 1-Wasserstein distance ${𝒲}_1(d)$ under the pseudometric $d$ can be computed using the dual linear program:

$$\underset{𝐮,𝐯}{\max }\sum _{x\in 𝒳}P(x)u_x-\sum _{y\in 𝒴}P(y)v_y\mathrm{subject}\mathrm{to}\forall x\in 𝒳,y\in 𝒴u_x-v_y\le d(x,y)$$

Since we are only interested in computing the coupling between the states in $𝒳$ and $𝒴$, the above formulation assumes that $P_{𝒳}(y)=0$ for all $y\in 𝒴$ and $P_{𝒴}(x)=0$ for all $x\in 𝒳$. The computation of $d_{\sim }$ is expensive and requires a tabular representation of the states, rendering it impractical for large state spaces. On-policy bisimulation (Castro, 2020) (e.g., ${\pi }^{\ast }$-bisimulation) is tied to specific behavior policies and is much easier to approximate than bisimulation.

Appendix C Policy Similarity Metric

C.1 Computing PSM

In general, PSM for a given Dist across MDPs ${ℳ}_{𝒳}$ and ${ℳ}_{𝒴}$ is given by

$d^{\ast }(x,y)=\text{Dist}({\pi }_{𝒳}^{\ast }(x),{\pi }_{𝒴}^{\ast }(y))+\gamma {𝒲}_1(d^{\ast })(P_{𝒳}^{{\pi }^{\ast }}(\cdot |x),P_{𝒴}^{{\pi }^{\ast }}(\cdot |y)).$ (C.1)

Since our main focus is showing the utility of PSM for generalization, we simply use environments where PSM can be computed using dynamic programming. Using a similar observation to Castro (2020), we assert that the recursion for $d^{\ast }$ takes the following form in deterministic environments:

$d^{\ast }(x,y)=\text{Dist}({\pi }_{𝒳}^{\ast }(x),{\pi }_{𝒴}^{\ast }(y))+\gamma d^{\ast }(x^{\prime },y^{\prime }).$ (C.2)

where $x^{\prime }=P_{𝒳}^{{\pi }^{\ast }}(x)$, $y^{\prime }=P_{𝒴}^{{\pi }^{\ast }}(y)$ are the next states from taking actions ${\pi }_{𝒳}^{\ast }(x)$, ${\pi }_{𝒳}^{\ast }(y)$ from states $x,y$ respectively. Furthermore, we assume that Dist between terminal states from ${ℳ}_{𝒳}$ and ${ℳ}_{𝒴}$ is zero. Note that the form of Equation C.2 closely resembles the update rule in Q-learning, and as such, can be efficiently computed with samples using approximate dynamic programming. Given access to optimal trajectories ${\tau }_{𝒳}^{\ast }={\{x_t\}}_{t=1}^N$ and ${\tau }_{𝒴}^{\ast }={\{y_t\}}_{t=1}^N$, where $x_{t+1}=P_{𝒳}^{{\pi }^{\ast }}(x_t)$ and $y_{t+1}=P_{𝒴}^{{\pi }^{\ast }}(y_t)$, Equation C.2 can be solved using exact dynamic programming; we provide pseudocode in Section I.1.

There are other ways to approximate the Wasserstein distance in bisimulation metrics (e.g., Ferns et al., 2006; 2011; Castro, 2020; Zhang et al., 2021). That said, approximating bisimulation (or PSM) for stochastic environments remains an exciting research direction (Castro, 2020). Investigating other distance metrics for long-term behavior difference in PSM is also interesting for future work.

C.2 PSM Connections to Data Augmentation and Bisimulation

Connection to bisimulation. Although bisimulation metrics have appealing properties such as bounding value function differences (e.g., (Ferns & Precup, 2014)), they rely on reward information and may not provide a meaningful notion of behavioral similarity in certain environments. Proposition C.1 implies that states similar under PSM would have similar optimal policies yet can have arbitrarily large bisimulation distance between them.

Proposition C.1.

There exists environments $M_{𝒳}$ and $M_{𝒴}$ such that $\forall (x,y)\in ℒ$ where $ℒ=\{(x,y)|x\in 𝒳,y\in 𝒴,d^{\ast }(x,y)=0\}$, $d_{\sim }^{\ast }(x,y)=\frac{|R_{\max }-R_{\min }|}{1-\gamma }-ϵ$ for any given $ϵ>0$.

For example, consider the two semantically equivalent environments in Figure 2 with ${\pi }_{𝒳}^{\ast }(x_0)={\pi }_{𝒴}^{\ast }(y_0)=a_0$ and ${\pi }_{𝒳}^{\ast }(x_1)={\pi }_{𝒴}^{\ast }(y_1)=a_1$ but different rewards $r_x,r_y$ respectively. Whenever $r_y>(1+1/\gamma )r_x$, bisimulation metrics incorrectly imply that $x_0$ is more behaviorally similar to $y_1$ than $y_0$.

For the MDPs shown in Figure 2, to determine which $y$ state is behaviorally equivalent to $x_0$, we look at the distances computed by bisimulation metric $d_{\sim }$ and ${\pi }^{\ast }$-bisimulation metric $d_{\sim }^{\ast }$:

$d_{\sim }(x_0,y_0)$ $=d_{\sim }^{\ast }(x_0,y_0)=(1+\gamma )|r_y-r_x|$

$d_{\sim }(x_0,y_1)$ $=\max ((1+\gamma )r_x,r_y),d_{\sim }^{\ast }(x_0,y_1)=|r_y-r_x|+\gamma r_x$

Thus, $r_y>(1+1/\gamma )r_x$ implies that as well as .

Connection to data augmentation. Data augmentation often assumes access to optimality invariant transformations, e.g., random crops or flips in image-based benchmarks (Laskin et al., 2020a; Kostrikov et al., 2020). However, for certain RL tasks, such augmentations can erroneously alias states with different optimal behavior and hurt generalization. For example, if the image observation is flipped in a goal reaching task with left and right actions, the optimal actions would also be flipped to take left actions instead of right and vice versa. Proposition C.2 states that PSMs can precisely quantify the invariance of such augmentations.

Proposition C.2.

For an MDP $M_{𝒳}$ and its transformed version $M_{\psi (𝒳)}$ for the data augmentation $\psi$, $d^{\ast }(x,\psi (x))$ indicates the optimality invariance of $\psi$ for any $x\in 𝒳$.

C.3 PSM with Approximately-Optimal Policies

Generalized Policy Similarity Metric for arbitrary policies. For a given Dist, we define a generalized PSM $d:(𝒮\times \mathrm{\Pi })\times (𝒮\times \mathrm{\Pi })\to ℝ$ where $\mathrm{\Pi }$ is the set of all policies over $𝒮$. $d$ satisfies the recursive equation:

$$d((x,{\pi }_1),(y,{\pi }_2))=\text{Dist}({\pi }_1(x),{\pi }_2(y))+\gamma {𝒲}_1(d)(P^{{\pi }_1}(\cdot |x),P^{{\pi }_2}(\cdot |y)).$$ (C.3)

Since Dist is assumed to be a pseudometric and ${𝒲}_1$ is a probability metric, it implies that $d$ is a pseudometric as (1) $d$ is non-negative, that is, $d((x,{\pi }_1),(y,{\pi }_2))\ge 0$, (2) $d$ is symmetric, that is, $d((x,{\pi }_1),(y,{\pi }_2))=d((x,{\pi }_1),(y,{\pi }_2))$, and $d$ satisfies the triangle inequality, that is, .

Using this notion of generalized PSM, we show that the approximation error in PSM from using a suboptimal policy is bounded by the policy’s suboptimality. Thus, for policies with decreasing suboptimality, the PSM approximation becomes more accurate, resulting in improved PSEs.

Proposition C.3.

[Approximation error in PSM] Let $\widehat{d}:𝒮\times 𝒮\to ℝ$ be the approximate PSM computed using a suboptimal policy $\widehat{\pi }$ defined over $𝒮$, that is, $\widehat{d}(x,y)=\text{Dist}(\widehat{\pi }(x),\widehat{\pi }(y))+\gamma {𝒲}_1(\widehat{d})(P^{\widehat{\pi }}(\cdot |x),P^{\widehat{\pi }}(\cdot |y))$ . We have:

Proof.

The PSM $d^{\ast }$ and approximate PSM $\widehat{d}$ are instantiations of the generalized PSM (Equation C.3) with both input policies as ${\pi }^{\ast }$ and $\widehat{\pi }$ respectively.

$d^{\ast }(x,y)=d((x,{\pi }^{\ast }),(y,{\pi }^{\ast }))$

$d^{\ast }(x,y)-\widehat{d}(x,y)$

$\text{Similarly},\widehat{d}(x,y)-d^{\ast }(x,y)$

∎

Appendix D L2 Metric Embeddings

Another common choice (Zhang et al., 2021) for learning metric embeddings is to use the squared loss (i.e., $l2$-loss) for matching the euclidean distance between the representations of a pair of states to the metric distance between those states. Concretely, for a given $d^{\ast }$ and representation $f_{\theta }$, the loss $ℒ(\theta )={𝔼}_{s_i,s_j}[{({\Vert f_{\theta }(s_i)-f_{\theta }(s_j)\Vert }_2-d^{\ast }(s_i,s_j))}^2]$ is minimized. However, it might be too restrictive to match the exact metric distances, which we demonstrate empirically by comparing $l2$ metric embeddings with CMEs (Section 5.2).

Appendix E Extended Related Work

Generalization across different tasks used to be described as transfer learning. In the past, most transfer learning approaches relied on fixed representations and tackled different problem formulations (e.g., assuming shared state space). Taylor & Stone (2009) present a comprehensive survey of the techniques at the time, before representation learning became so prevalent in RL. Recently, the problem of performing well in a different, but related task, started to be seen as a problem of generalization; with the community highlighting that deep RL agents tend to overfit to the environments they are trained on (Cobbe et al., 2019; Witty et al., 2018; Farebrother et al., 2018; Juliani et al., 2019; Kostrikov et al., 2020; Song et al., 2019; Justesen et al., 2018; Packer et al., 2018).

Prior generalization approaches are typically adapted from supervised learning, including regularization (Cobbe et al., 2019; Farebrother et al., 2018), stochasticity (Zhang et al., 2018c), noise injection (Igl et al., 2019; Zhang et al., 2018a), more diverse training conditions (Rajeswaran et al., 2017; Witty et al., 2018) and self-attention architectures (Tang et al., 2020). In contrast, PSEs exploits behavior similarity (Section 3), a property related to the sequential aspect of RL.

Meta-learning is also related to generalization. Meta-learning methods try to find a parametrization that requires a small number of gradient steps to achieve good performance on a new task (Finn et al., 2017). In this context, various meta-learning approaches capable of zero-shot generalization have been proposed (Li et al., 2018; Agarwal et al., 2019; Balaji et al., 2018). These approaches typically consist in minimizing the loss in the environments the agent is while adding an auxiliary loss for ensuring improvement in the other (validation) environments available to the agent. Nevertheless, Tachet des Combes et al. (2018) has shown that meta-learning approaches fail in the jumping task which we also observed empirically. Others have also reported similar findings (e.g., Farebrother et al., 2018).

There are several other approaches for tackling zero-shot generalization in RL, but they often rely on domain-specific information. Some examples include knowledge about equivalences between entities in the environment (Oh et al., 2017) and about what is under the agent’s control (Ye et al., 2020). Causality-based methods are a different way of tackling generalization, but current solutions do not scale to high-dimensional observation spaces (e.g., Killian et al., 2017; Perez et al., 2020; Zhang et al., 2020).

Appendix F Jumping Task with Pixels

Figure F.1: Optimal trajectories on the jumping tasks for two different environments. Note that the optimal trajectory is a sequence of right actions, followed by a single jump at a certain distance from the obstacle, followed by right actions.

Detailed Task Description. The jumping task consists of an agent trying to jump over an obstacle on a floor. The environment is deterministic with the agent observing a reward of $+1$ at each time step. If the agent successfully reaches the rightmost side of the screen, it receives a bonus reward of $+100$; if the agent touches the obstacle, the episode terminates. The observation space is the pixel representation of the environment, as depicted in Figure 1. The agent has access to two actions: right and jump. The jump action moves the agent vertically and horizontally to the right.

(a) Bisimulation metric

(b) Policy Similarity metric

(c) ${\pi }^{\ast }$-Bisimulation metric

(d) True State Similarity

Figure F.2: PSM vs. Bisimulation. Visualizing PSM and bisimulation metrics with discount factor $\gamma =0.99$. $x_i$ and $y_i$ correspond to the states visited by optimal policies in the two environments with obstacle at positions 25 and 45, respectively. For each grid, the ${(i,j)}^{th}$ location shows the distance assigned by the metric to the states $x_i$ and $y_j$. Lower distances (lighter shades) imply higher similarity between two states. The bottom right figure shows which states are equivalent to each other. The ranges are different for each metric as bisimulation metrics utilize reward differences while PSM uses the total variation ($TV$) distance between policies. Note that the large bisimulation distances are due to the fact that the reward at the terminal state is set to 100.

Architecture. The neural network used for Jumping Task experiment is adapted from the Nature DQN architecture. Specifically, the network consists of 3 convolutional layers of sizes 32, 64, 64 with filter sizes $8\times 8$, $4\times 4$ and $3\times 3$ and strides 4, 2, and 1, respectively. The output of the convnet is fed into a single fully connected layer of size 256 followed by ’ReLU’ non-linearity. Finally, this FC layer output is fed into a linear layer which computes the policy which outputs the probability of the jump and right actions.

Contrastive Embedding. For all our experiments, we use a single ReLU layer with $k=64$ units for the non-linear projection to obtain the embedding $z_{\theta }$ (Figure 3). We compute the embedding using the penultimate layer in the jumping task network. Hyperparameters are reported in Table F.2.

Total Loss. For jumpy world, the total loss is given by ${ℒ}_{\text{IL}}+\alpha {ℒ}_{\mathrm{CME}}$ where ${ℒ}_{\text{IL}}$ is the imitation learning loss, ${ℒ}_{\mathrm{CME}}$ is the auxiliary loss for learning PSEs with the coefficient $\alpha$.

Figure F.3: Test performance curves in the setting without data augmentation on the “wide”, “narrow”, and random grids described in Figure 4. We plot the median performance across 100 runs. Shaded regions show 25 and 75 percentiles.

Figure F.4: Test performance curves in the setting with data augmentation on the “wide”, “narrow”, and random grids described in Figure 4. We plot the median performance across 100 runs. Shaded regions show 25 and 75 percentiles.

F.1 Jumping Task from Colors

Figure F.5: Jumping task with colored obstacles. The agent needs to jump over the red obstacle but strike the green obstacle.

Table F.1: Percentage (%) of test tasks solved when trained on the “wide” grid with both red and green obstacles. The numbers we report are averaged across 100 runs. Standard error is reported between parentheses.

Method	Red (%)		Green (%)
RandConv	6.2	(0.4)	99.6	(0.2)
Dropout and $l_2$ reg.	19.5	(0.2)	100.0	(0.0)
RandConv + PSEs	29.8	(1.3)	99.6	(0.2)
PSEs	37.9	(1.9)	100.0	(0.0)

Figure F.6: Randconv enforces color invariance. The first column shows the original observation where the top row corresponds to the observation from task with green obstacle jumping task while the bottom row corresponds to the red obstacle jumping task. Columns 2-6 in each row show 5 augmentations by applying RandConv. The augmentations show that RandConv tries to encode invariance with respect to the obstacle color.

F.2 Hyperparameters

For hyperparameter selection, we evaluate all agents on a validation set containing 54 unseen tasks in the “wide” grid and pick the parameters with the best validation performance. The validation set (Figure F.7) was selected by using the environments nearby to the training environments whose floor height differ by 1 or whose obstacle position differ by 1.

Table F.2: Common hyperparameters across all methods for all jumping task experiments.

Hyperparameter	Value
Learning rate decay	0.999
Training epochs	2000
Optimizer	Adam
Batch size (Imitation)	256
Num training tasks	18
$\mathrm{\Gamma }$-scale Parameter ($\beta$)	0.01
Embedding size ($k$)	64
Batch Size (${ℒ}_{\mathrm{CME}}$)	57
$|{\tau }_{𝒳}^{\ast }|$	57

Figure F.7: Unseen 56 validation tasks, labeled with V, used for hyperparameter selection.

Table F.3: Optimal hyperparameters for reporting results in Table 1. These hyperparameters are selected using the “wide” grid by maximizing final performance on a validation set containing 56 unseen tasks. All grid configurations in Table 1 use these hyperparameters.

Hyperparameter	Dropout and ${\mathrm{\ell }}_2$-reg.	PSEs	RandConv	RandConv + PSEs
Learning Rate	$4\text{×}{10}^{-3}$	$3.2\text{×}{10}^{-3}$	$7\text{×}{10}^{-3}$	$2.6\text{×}{10}^{-3}$
${\mathrm{\ell }}_2$-reg. coefficient	$4.3\text{×}{10}^{-4}$	$1\text{×}{10}^{-5}$	–	–
Dropout coefficient	$3\text{×}{10}^{-1}$	–	–	–
Contrastive Temperature ($1/\lambda$)	–	$1.0$	–	$5\text{×}{10}^{-1}$
Auxiliary loss coefficient ($\alpha$)	–	$1\text{×}{10}^1$	–	$5.0$

Table F.4: Optimal hyperparameters for reporting results in Figure 5.4. These hyperparameters are selected using the “wide” grid by maximizing final performance on a validation set containing 56 unseen tasks.

Hyperparameter	Dropout and ${\mathrm{\ell }}_2$-reg.	PSEs	RandConv	RandConv + PSEs
Learning Rate	$4\text{×}{10}^{-3}$	$6\text{×}{10}^{-3}$	$5\text{×}{10}^{-3}$	$2.6\text{×}{10}^{-3}$
${\mathrm{\ell }}_2$-reg. coefficient	$4.3\text{×}{10}^{-4}$	$7\text{×}{10}^{-5}$	–	–
Dropout coefficient	$3\text{×}{10}^{-1}$	–	–	–
Contrastive Temperature ($1/\lambda$)	–	$5\text{×}{10}^{-1}$	–	$5\text{×}{10}^{-1}$
Auxiliary loss coefficient ($\alpha$)	–	$5.0$	–	$5.0$

Table F.5: Optimal hyperparameters for reporting ablation results in Table 2. These hyperparameters are selected using the “wide” grid by maximizing final performance on a validation set containing 56 unseen tasks.

Hyperparameter	PSM		${\pi }^{\ast }$-bisimulation
	CMEs	${\mathrm{\ell }}_2$-embeddings	CMEs	${\mathrm{\ell }}_2$-embeddings
Learning Rate	$4\text{×}{10}^{-3}$	$5\text{×}{10}^{-4}$	$4.7\text{×}{10}^{-4}$	$1\text{×}{10}^{-4}$
Contrastive Temperature ($1/\lambda$)	$1.0$	–	$5\text{×}{10}^{-1}$	–
Auxiliary loss coefficient ($\alpha$)	$5.0$	$1\text{×}{10}^{-1}$	$1\text{×}{10}^{-1}$	$1\text{×}{10}^{-6}$

Please note that Table F.3 and Table F.4 correspond to two different tasks: one uses the standard jumping task with white obstacles, while the other uses colored obstacles where the optimal policies depend on color. For fair comparison, we tune hyperparameters for all the methods using Bayesian optimization (Golovin et al., 2017). We use the best parameters among these tuned hyperparameters and the ones found in Table F.3, leading to different parameters for both PSEs as well as RandConv. Evaluating PSEs with the jumping task hyperparameters from Table F.3 instead of the ones in Table F.4 leads to a small drop (-4%) on the jumping task with colors (Section 5.4). Nevertheless, PSEs still outperform other methods in Section 5.4.

Appendix G LQR: Additional Details

Table G.1: LQR generalization performance: Absolute error in LQR cost, w.r.t. the oracle solver (which has access to true state), of various methods trained with $n_d$ distractors on $N=2$ environments. The reported mean and standard deviations are across 100 different seeds. Lower error is better.

Method	Number of Distractors ($n_d$)
	500	1000	10000
Overparametrization (Song et al., 2019)	25.8 (1.5)	24.9 (1.1)	24.9 (0.4)
IPO (Sonar et al., 2020) (IRM + Policy opt.)	32.6 (5.0)	27.3 (2.8)	24.8 (0.4)
Weight Sparsity (${\mathrm{\ell }}_1$-reg.)	28.2 (0.0)	28.2 (0.0)	28.2 (0.0)
PSM (State aggregation)	0.03 (0.0)	0.03 (0.0)	0.02 (0.0)

Optimal control with linear dynamics and quadratic cost, commonly known as LQR, has been increasingly used as a simplified surrogate for deep RL problems (Recht, 2019). Following Song et al. (2019); Sonar et al. (2020), we analyze the following LQR problem for assessing generalization:

(G.1)

where $𝒟$ is the initial state distribution, $s_t\in {ℝ}^{n_s}$ is the (hidden) true state at time t, $a_t\in {ℝ}^{n_a}$ is the control action, and $K$ is the linear policy matrix. The agent receives the input observation $o_t$, which is a linear transformation of the state $s_t$. $W_c$ and $W_d$ are semi-orthogonal matrices to prevent information loss for predicting optimal actions. An environment corresponds to a particular choice of $W_d$; all other system parameters ($A,B,Q,R,W_c$) are fixed matrices which are shared across environments and unknown to the agent. The agent learns the policy matrix $K$ using $N$ training environments based on Equation G.1. At test time, the learned policy is evaluated on environments with unseen $W_d$.

The generalization challenge in this setup is to ignore the distractors: $W_c{s}_t\in {ℝ}^{n_s}$ represents the state features invariant across environments while $W_d{s}_t\in {ℝ}^{n_d}$ is a high-dimensional distractor of size $n_s,n_d$, respectively, such that . Furthermore, the policy matrix which generalizes across all environments is $K_{\star }={\left[\begin{array}{c}10W_{c}P_{\star }^{}{}_{}{}^T\\ 0\end{array}\right]}^T$, where $P_{\star }$ corresponds to the optimal LQR solution with access to state $s_t$. However, for a single environment with distractor $W_d$, multiple solutions exist, for instance, $K_{\star }^{}{}_{}{}^{\prime }={\left[\begin{array}{c}10\alpha W_{c}P_{\star }^{}{}_{}{}^T\\ (1-\alpha )W_{d}P_{\star }^{}{}_{}{}^T\end{array}\right]}^T\forall \alpha \in [0,1]$. Note that the distractors are an order of magnitude larger than the invariant features in $o_t$ and dependence on them is likely to cause the agent to act erratically on inputs with unseen distractors, resulting in poor generalization.

We use overparametrized policies with two linear layers, i.e., $K=K_1{K}_2$, where $K_1(o)$ is the learned representation for observation $o$. We learn $K$ using gradient descent using the combined cost on 2 training environments with varying number of distractors. We aggregate observation pairs with near-zero PSM by matching their representations using a squared loss. We use the open-source code released by Sonar et al. (2020) for our experiments.

Table G.2: An overview of hyper-parameters for LQR.

Parameter	Setting
A	Orthogonal matrix, scaled 0.8
B	$I_{20\times 20}$
$n_x$	20
$n_a$	20
Q	$I_{20\times 20}$
R	$I_{20\times 20}$
$K_i\forall i$	Orthogonal Initialization, scaled 0.001
$W_d$	Random semi-orthogonal matrix

The reliance on distractors for IPO also highlights a limitation of IRM: if a model can achieve a solution with zero training error, then any such solution is acceptable by IRM regardless of its generalization ability – a common scenario with overparametrized deep neural nets (Jin et al., 2020).

G.1 Near-optimality of PSM Aggregation

Conjecture 1.

Assuming zero state aggregation error with policy similarity metric (PSM), the policy matrix $K$ learned using gradient descent is independent of the distractors.

Proof.

For LQR domains $x,y$, an observation pair ($o_t^{}{}_{}{}^x$, $o_t^{}{}_{}{}^y$) has zero PSM iff the underlying state $s_t$ is same for both the observations in the pair. This is true, as (a) both domains has the same transition dynamics, as specified by Equation G.1, and (b) the optimal policy is deterministic and is completely determined by the current state $s_t$ at any time $t$.

Assume $o_t^{}{}_{}{}^x=\left[\begin{array}{c}0.1W_c\\ W_{d^x}\end{array}\right]s_t$ and $o_t^{}{}_{}{}^y=\left[\begin{array}{c}0.1W_c\\ W_{d^y}\end{array}\right]s_t$ for distractor semi-orthogonal matrices $W_{d^{𝒳}}$ and $W_{d^{𝒴}}$, respectively. Furthermore, the representation is given by $K_1(o_t^{}{}_{}{}^x)$ and $K_1(o_t^{}{}_{}{}^y)$ respectively. Assume that where $K_s\in {ℝ}^{h\times n_s}$ and $K_d\in {ℝ}^{h\times n_d}$ and $K_1\in {ℝ}^{h\times (n_s+n_d)}$.

Zero state-aggregation error with squared loss implies that for pair ($o_t^{}{}_{}{}^x$, $o_t^{}{}_{}{}^y$) corresponding to $s_t$,

$$K_1(o_t^{}{}_{}{}^x-o_t^{}{}_{}{}^y)=K_1\left[\begin{array}{c}0\\ W_{d^x}-W_{d^y}\end{array}\right]s_t=0⟹K_d(W_{d^x}-W_{d^y})s_t=0$$ (G.2)

As Equation G.2 holds for all states visited by the optimal policy in an infinite horizon LQR, it follows that $K_d(W_{d^x}-W_{d^y})=0$.

Furthermore, it is well-known that gradient descent tends to find low-rank solutions due to implicit regularization (Arora et al., 2019; Gunasekar et al., 2017), e.g.,with small enough step sizes and initialization close enough to the origin, gradient descent on matrix factorization converges to the minimum nuclear norm solution for 2 layer linear networks (Gunasekar et al., 2017). Based on this, we conjecture that $K_d=0$ which we found to be true in practice. ∎

Appendix H Distracting Control Suite

Figure H.1: DCS Test Environments: Snapshots of test environments used for evaluating generalization on Distracting Control Suite. Random backgrounds are projected from scenes of the first 30 videos of DAVIS 2017 validation dataset.

We use the same setup as Kostrikov et al. (2020); Stone et al. (2021) for implementation details and training protocol. For completeness, we describe the details below.

Dynamic Background Distractions. In Distracting Control Suite (Stone et al., 2021), random backgrounds are projected from scenes of the DAVIS 2017 dataset (Pont-Tuset et al., 2017) onto the skybox of the scene. To make these backgrounds visible for all tasks and views, the floor grid is semi-transparent with transparency coefficient 0.3. We take the first 2 videos in the DAVIS 2017 training set and randomly sample a scene and a frame from those at the start of every episode. In the dynamic setting, the video plays forwards or backwards until the last or first frame is reached at which point the video is played backwards. This way, the background motion is always smooth and without “cuts”.

Soft Actor-Critic. Soft Actor-Critic (SAC) (Haarnoja et al., 2018) learns a state-action value function $Q_{\theta }$, a stochastic policy ${\pi }_{\theta }$ and a temperature $\alpha$ to find an optimal policy for an MDP $(𝒮,𝒜,p,r,\gamma )$ by optimizing a $\gamma$-discounted maximum-entropy objective. $\theta$ is used generically to denote the parameters updated through training in each part of the model. The actor policy ${\pi }_{\theta }(a_t|s_t)$ is a parametric $\tanh$-Gaussian that given $s_t$ samples $a_t=\tanh ({\mu }_{\theta }(s_t)+{\sigma }_{\theta }(s_t)ϵ)$, where $ϵ\sim 𝒩(0,1)$ and ${\mu }_{\theta }$ and ${\sigma }_{\theta }$ are parametric mean and standard deviation.

The policy evaluation step learns the critic $Q_{\theta }(s_t,a_t)$ network by optimizing a single-step of the soft Bellman residual

$J_Q(𝒟)$ $=E_{\begin{array}{c}(s_t,a_t,s_t^{\prime })\sim 𝒟\\ a_t^{\prime }\sim \pi (\cdot |s_t^{\prime })\end{array}}[{(Q_{\theta }(s_t,a_t)-y_t)}^2]$

$y_t$ $=r(s_t,a_t)+\gamma [Q_{{\theta }^{\prime }}(s_t^{\prime },a_t^{\prime })-\alpha \log {\pi }_{\theta }(a_t^{\prime }|s_t^{\prime })],$

where $𝒟$ is a replay buffer of transitions, ${\theta }^{\prime }$ is an exponential moving average of the weights. SAC uses clipped double-Q learning, which we omit for simplicity but employ in practice.

The policy improvement step then fits the actor policy ${\pi }_{\theta }(a_t|s_t)$ network by optimizing the objective

$J_{\pi }(𝒟)$ $=E_{s_t\sim 𝒟}[D_{\mathrm{KL}}({\pi }_{\theta }(\cdot |s_t)||\exp \{{\displaystyle \frac{1}{\alpha }}{Q}_{\theta }(s_t,\cdot )\})].$

Finally, the temperature $\alpha$ is learned with the loss

$J_{\alpha }(𝒟)$ $=E_{\begin{array}{c}{s}_t\sim 𝒟\\ a_t\sim {\pi }_{\theta }(\cdot |s_t)\end{array}}[-\alpha \log {\pi }_{\theta }(a_t|s_t)-\alpha \overline{\mathscr{H}}],$

where $\overline{\mathscr{H}}\in ℝ$ is the target entropy hyper-parameter that the policy tries to match, which in practice is usually set to $\overline{\mathscr{H}}=-|𝒜|$.

H.1 Actor and Critic Networks

Following Kostrikov et al. (2020), we use clipped double Q-learning for the critic, where each $Q$-function is parametrized as a 3-layer MLP with ReLU activations after each layer except of the last. The actor is also a 3-layer MLP with ReLUs that outputs mean and covariance for the diagonal Gaussian that represents the policy. The hidden dimension is set to $1024$ for both the critic and actor.

H.2 Encoder Network

We employ the encoder architecture from Kostrikov et al. (2020). This encoder consists of four convolution layers with $3\times 3$ kernels and $32$ channels. The ReLU activation is applied after each convolutional layer. We use stride to $1$ everywhere, except of the first convolutional layer, which has stride $2$. The output of the convnet is feed into a single fully-connected layer normalized by LayerNorm. Finally, we apply tanh nonlinearity to the $50$ dimensional output of the fully-connected layer. We initialize the weight matrix of fully-connected and convolutional layers with the orthogonal initialization and set the bias to be zero. The actor and critic networks both have separate encoders, although we share the weights of the conv layers between them. Furthermore, only the critic optimizer is allowed to update these weights (i.e.,we stop the gradients from the actor before they propagate to the shared convolutional layers).

H.3 Contrastive Metric Embedding Loss

For all our experiments, we use a single ReLU layer with $k=256$ units for the non-linear projection to obtain the embedding $z_{\theta }$ (Figure 3). We compute the embedding using the penultimate layer in the actor network. For picking the hyperparameters, we used 3 temperatures $[0.1,0.01,1.0]$ and 3 auxiliary ${ℒ}_{\mathrm{CME}}$ loss coefficients $[1,3,10]$ using “Ball In Cup Catch” as the validation environment. All other hyperparameters are the same as prior work (see Table H.2).

We approximate optimal policies with the policies obtained after training a DrQ agent for 500K environment steps. Since a given action sequence from this approximate policy has the same performance across different training environments, we compute the PSM across training environments, via dynamic programming (see Section I.1 for pseudo-code), using such action sequences.

Total Loss. The total loss is given by ${ℒ}_{\text{RL}}+\alpha {ℒ}_{\mathrm{CME}}$ where ${ℒ}_{\text{RL}}$ is the reinforcement learning loss which combines $J_{\pi }(𝒟)$, $J_{\pi }(𝒟)$, and $J_{\alpha }(𝒟)$), while ${ℒ}_{\mathrm{CME}}$ is the auxiliary loss for learning PSEs with the coefficient $\alpha$.

H.4 Training and Evaluation Setup

For evaluation, we use the first 30 videos from the DAVIS 2017 validation dataset (see Figure H.1). Each checkpoint is evaluated by computing the average episode return over 100 episodes from the unseen environments. All experiments are performed with five random seeds per task used to compute means and standard deviations/errors of their evaluations. We use $K=2,M=2$ as prescribed by Kostrikov et al. (2020) for DrQ. Following Kostrikov et al. (2020) and Stone et al. (2021), we use a different action repeat hyper-parameter for each task, which we summarize in Table H.3. We construct an observational input as a $3$-stack of consecutive frames (Kostrikov et al., 2020), where each frame is an RGB rendering of size $84\times 84$ from the $0$th camera. We then divide each pixel by $255$ to scale it down to $[0,1]$ range. For data augmentation, we maintain temporal consistency by using the same crop augmentation across consecutive frames.

Table H.1: Optimal hyperparameters for PSE auxiliary loss for reporting results in Table 3.

Hyperparameter	Setting
Contrastive temperature ($1/\lambda$)	0.1
Auxiliary loss coefficient ($\alpha$)	1.0
$\mathrm{\Gamma }$-scale parameter ($\beta$)	0.1
Batch Size (${ℒ}_{\mathrm{CME}}$)	128
$|{\tau }_{𝒳}^{\ast }|$	1000 // Action Repeat

Table H.2: Hyper-parameters taken from Kostrikov et al. (2020) in the Distracting Control Suite experiments.

Parameter	Setting
Replay buffer capacity	$100,000$
Seed steps	$1,000$
Batch size (DrQ)	$512$
Discount $\gamma$	$0.99$
Optimizer	Adam
Learning rate	${10}^{-3}$
Critic target update frequency	$2$
Critic Q-function soft-update rate $\tau$	$0.01$
Actor update frequency	$2$
Actor log stddev bounds	$[-10,2]$
Init temperature	$0.1$

Table H.3: The action repeat hyper-parameter used for each task in the Distracting Control Suite benchmark.

Task name	Action repeat
Cartpole Swingup	$8$
Reacher Easy	$4$
Cheetah Run	$4$
Finger Spin	$2$
Ball In Cup Catch	$4$
Walker Walk	$2$

H.5 Generalization Curves

Figure H.2: Random Initialization. Generalization performance on unseen environments over the course of training. The agent is initialized randomly. DrQ augmented with PSEs outperforms DrQ. We plot the average episode return across 5 seeds and the shaded region shows the standard deviation. Each checkpoint is evaluated using 100 episodes with unseen backgrounds.

Figure H.3: Pretrained Initialization. Generalization performance on unseen environments over the course of training. The agent is initialized using a pretrained DrQ agent. DrQ augmented with PSEs outperforms DrQ on most of the environments. We plot the average return across 5 seeds and the shaded region shows the standard deviation. Each checkpoint is evaluated using 100 episodes with unseen backgrounds.

Appendix I Pseudo code

I.1 Dynamic Programming for Computing PSM

⬇

1def metric_fixed_point(cost_matrix, gamma=0.99, eps=1e-7):

2 """DP for calculating PSM in environments with deterministic dynamics.

3

4 Args:

5 cost_matrix: DIST matrix where entries at index (i, j) is DIST(x_i, y_j)

6 gamma: Metric discount factor.

7 eps: Threshold for stopping the fixed point iteration.

8 """

9 d = np.zeros_like(cost_matrix)

10 def operator(d_cur):

11 d_new = 1 * cost_matrix

12 discounted_d_cur = gamma * d_cur

13 d_new[:-1, :-1] += discounted_d_cur[1:, 1:]

14 d_new[:-1, -1] += discounted_d_cur[1:, -1]

15 d_new[-1, :-1] += discounted_d_cur[-1, 1:]

16 return d_new

17

18 while True:

19 d_new = operator(d)

20 if np.sum(np.abs(d - d_new)) < eps:

21 break

22 else:

23 d = d_new[:]

24 return d

I.2 Contrastive Loss

⬇

1def contrastive_loss(similarity_matrix,

2 metric_values,

3 temperature,

4 beta=1.0):

5 """Contrative Loss with embedding similarity ."""

6 metric_shape = tf.shape(metric_values)

7 ## z_\theta(X): embedding_1 = nn_model.representation(X)

8 ## z_\theta(Y): embedding_2 = nn_model.representation(Y)

9 ## similarity_matrix = cosine_similarity(embedding_1, embedding_2

10 ## metric_values = PSM(X, Y)

11 similarity_matrix /= temperature

12 neg_logits1 = similarity_matrix

13

14 col_indices = tf.cast(tf.argmin(metric_values, axis=1), dtype=tf.int32)

15 pos_indices1 = tf.stack(

16 (tf.range(metric_shape[0], dtype=tf.int32), col_indices), axis=1)

17 pos_logits1 = tf.gather_nd(similarity_matrix, pos_indices1)

18

19 metric_values /= beta

20 similarity_measure = tf.exp(-metric_values)

21 pos_weights1 = -tf.gather_nd(metric_values, pos_indices1)

22 pos_logits1 += pos_weights1

23 negative_weights = tf.math.log((1.0 - similarity_measure) + 1e-8)

24 neg_logits1 += tf.tensor_scatter_nd_update(

25 negative_weights, pos_indices1, pos_weights1)

26

27 neg_logits1 = tf.math.reduce_logsumexp(neg_logits1, axis=1)

28 return tf.reduce_mean(neg_logits1 - pos_logits1) # Equation 4

a