Title: THE CRAMÉR DISTANCE AS A SOLUTION TO BIASED WASSERSTEIN GRADIENTS

Abstract: The Wasserstein probability metric has received much attention from the machine learning community. Unlike the Kullback-Leibler divergence, which strictly measures change in probability, the Wasserstein metric reflects the underlying geometry between outcomes. The value of being sensitive to this geometry has been demonstrated, among others, in ordinal regression and generative modelling, and most recently in reinforcement learning. In this paper we describe three natural properties of probability divergences that we believe reflect requirements from machine learning: sum invariance, scale sensitivity, and unbiased sample gradients. The Wasserstein metric possesses the first two properties but, unlike the Kullback-Leibler divergence, does not possess the third. We provide empirical evidence suggesting this is a serious issue in practice. Leveraging insights from probabilistic forecasting we propose an alternative to the Wasserstein metric, the Cramér distance. We show that the Cramér distance possesses all three desired properties, combining the best of the Wasserstein and Kullback-Leibler divergences. We give empirical results on a number of domains comparing these three divergences. To illustrate the practical relevance of the Cramér distance we design a new algorithm, the Cramér Generative Adversarial Network (GAN), and show that it has a number of desirable properties over the related Wasserstein GAN.

1. INTRODUCTION:

In machine learning, the Kullback-Leibler (KL) divergence is perhaps the most common way of assessing how well a probabilistic model explains observed data. Among the reasons for its popularity is that it is directly related to maximum likelihood estimation and is easily optimized. However, the KL divergence suffers from a significant limitation: it does not take into account how close two outcomes might be, but only their relative probability. This closeness can matter a great deal: in image modelling, for example, perceptual similarity is key (Rubner et al., 2000; Gao & Kleywegt, 2016) . Put another way, the KL divergence cannot reward a model that "gets it almost right".

To address this limitation, researchers have turned to the Wasserstein metric, which does incorporate the underlying geometry between outcomes. The Wasserstein metric can be applied to distributions with non-overlapping supports, and has good out-of-sample performance (Esfahani & Kuhn, 2015) . Yet, practical applications of the Wasserstein distance, especially in deep learning, remain tentative. In this paper we provide a clue as to why that might be: estimating the Wasserstein metric from samples yields biased gradients, and may actually lead to the wrong minimum. This precludes using stochastic gradient descent (SGD) and SGD-like methods, whose fundamental mode of operation is sample-based, when optimizing for this metric.

As a replacement we propose the Cramér distance (Székely, 2002; Rizzo & Székely, 2016) , also known as the continuous ranked probability score in the probabilistic forecasting literature (Gneiting & Raftery, 2007) . The Cramér distance, like the Wasserstein metric, respects the underlying geometry but also has unbiased sample gradients. To underscore our theoretical findings, we demonstrate a significant quantitative difference between the two metrics when employed in typical machine learning scenarios: categorical distribution estimation, regression, and finally image generation. In the latter case, we use a multivariate generalization of the Cramér distance, the energy distance (Székely, 2002) , itself an instantiation of the MMD family of metrics (Gretton et al., 2012) .

2. PROBABILITY DIVERGENCES AND METRICS:

In this section we provide the notation to mathematically distinguish the Wasserstein metric (and later, the Cramér distance) from the Kullback-Leibler divergence and probability distances such as the total variation.

Let P be a probability distribution over R. When P is continuous, we will assume it has density µ P . The expectation of a function f : R → R with respect to P is

E x∼P f (x) := ∞ -∞

f (x)P (dx) = f (x)µ P (x)dx if P is continuous, and f (x)P (x) if P is discrete.

We will suppose all expectations and integrals under consideration are finite. We will often associate P to a random variable X, such that for a subset of the reals A ⊆ R, we have Pr{X ∈ A} = P (A).

The (cumulative) distribution function of P is then

F P (x) := Pr{X ≤ x} = x -∞ P (dx).

Finally, the inverse distribution function of P , defined over the interval (0, 1], is F -1 P (u) := inf{x : F P (x) = u}.

2.1. DIVERGENCES AND METRICS:

Consider two probability distributions P and Q over R. A divergence d is a mapping (P, Q) → R + with d(P, Q) = 0 if and only if P = Q almost everywhere. A popular choice is the Kullback-Leibler (KL) divergence

KL(P Q) := ∞ -∞ log P (dx) Q(dx) P (dx),

with KL(P Q) = ∞ if P is not absolutely continuous w.r.t. Q. The KL divergence, also called relative entropy, measures the amount of information needed to encode the change in probability from Q to P (Cover & Thomas, 1991) .

A probability metric is a divergence which is also symmetric (d(P, Q) = d(Q, P )) and respects the triangle inequality: for any distribution R, d(P, Q) ≤ d(P, R) + d(R, Q). We will use the term probability distance to mean a symmetric divergence satisfying the relaxed triangle inequality d(P, Q) ≤ c [d(P, R) + d(R, Q)] for some c ≥ 1.

We will first study the p-Wasserstein metrics w p (Dudley, 2002) . For 1 ≤ p < ∞, a practical definition is through the inverse distribution functions of P and Q:

EQUATION

We will sometimes find it convenient to deal with the p th power of the metric, which we will denote by w p p ; note that w p p is not a metric proper, but is a probability distance. We will be chiefly concerned with the 1-Wasserstein metric, which is most commonly used in practice. The 1-Wasserstein metric has a dual form which is theoretically convenient and which we mention here for completeness. Define F ∞ to be the class of 1-Lipschitz functions. Then

EQUATION

This is a special case of the celebrated Monge-Kantorovich duality (Rachev et al., 2013) , and is the integral probability metric (IPM) with function class F ∞ (Müller, 1997) . We invite the curious reader to consult these two sources as a starting point on this rich topic.

2.2. PROPERTIES OF A DIVERGENCE:

As noted in the introduction, the fundamental difference between the KL divergence and the Wasserstein metric is that the latter is sensitive not only to change in probability but also to the geometry of possible outcomes. To capture this notion we now introduce the concept of an ideal divergence.

Consider a divergence d, and for two random variables X, Y with distributions P, Q write d(X, Y ) := d(P, Q). We say that d is scale sensitive (of order β), i.e. it has property (S), if there exists a β > 0 such that for all X, Y , and a real value c > 0,

d(cX, cY ) ≤ |c| β d(X, Y ). (S)

A divergence d has property (I), i.e. it is sum invariant, if whenever A is independent from X, Y

d(A + X, A + Y ) ≤ d(X, Y ). (I)

Following Zolotarev (1976) , an ideal divergence d is one that possesses both (S) and (I).1 

We can illustrate the sensitivity of ideal divergences to the value of outcomes by considering Dirac functions δ x at different values of x. If d is scale sensitive of order β = 1 then the divergence d(δ 0 , δ 1/2 ) can be no more than half the divergence d(δ 0 , δ 1 ). If d is sum invariant, then the divergence of δ 0 to δ 1 is equal to the divergence of the same distributions shifted by a constant c, i.e. of δ c to δ 1+c . As a concrete example of the importance of these properties, Bellemare et al. (2017) recently demonstrated the importance of ideal metrics in reinforcement learning, specifically their role in providing the contraction property of the distributional Bellman operator. In particular, the contraction modulus is γ β , where γ ∈ [0, 1) is a discount factor and β is the scale sensitivity order.

In machine learning we often view the divergence d as a loss function. Specifically, let Q θ be some distribution parametrized by θ, and consider the loss θ → d(P, Q θ ). We are interested in minimizing this loss, that is finding θ * := arg min θ d(P, Q θ ). We now describe a third property based on this loss, which we call unbiased sample gradients.

Let X m := X 1 , X 2 , . . . , X m be independent samples from P and define the empirical distribution Pm := Pm (X m ) := 1 m m i=1 δ Xi (note that Pm is a random quantity). From this, define the sample loss θ → d( Pm , Q θ ). We say that d has unbiased sample gradients when the expected gradient of the sample loss equals the gradient of the true loss for all P and m:

E Xm∼P ∇ θ d( Pm , Q θ ) = ∇ θ d(P, Q θ ). (U)

The notion of unbiased sample gradients is ubiquitous in machine learning and in particular in deep learning. Specifically, if a divergence d does not possess (U) then minimizing it with stochastic gradient descent may not converge, or it may converge to the wrong minimum. Conversely, if d possesses (U) then we can guarantee that the distribution which minimizes the expected sample loss is Q = P . In the probabilistic forecasting literature, this makes d a proper scoring rule (Gneiting & Raftery, 2007) .

We now characterize the KL divergence and the Wasserstein metric in terms of these properties. As it turns out, neither simultaneously possesses both (U) and (S). Proposition 1. The KL divergence has unbiased sample gradients (U), but is not scale sensitive (S).

Proposition 2. The Wasserstein metric is ideal (I, S), but does not have unbiased sample gradients.

We will provide a proof of the bias in the sample Wasserstein gradients just below; the proof of the rest and later results are provided in the appendix.

3. BIAS IN THE SAMPLE GRADIENT ESTIMATES OF THE WASSERSTEIN DISTANCE:

In this section we give theoretical evidence of serious issues with gradients of the sample Wasserstein loss. We will consider a simple Bernoulli distribution P with parameter θ * ∈ (0, 1), which we would like to estimate from samples. Our model is Q θ , a Bernoulli distribution with parameter θ. We study the behaviour of stochastic gradient descent w.r.t. θ over the sample Wasserstein loss, specifically using the p th power of the metric (as is commonly done to avoid fractional exponents). Our results build on the example given by Bellemare et al. ( 2017), whose result is for θ * = 1 2 and m = 1.

Consider the estimate ∇ θ w p p ( Pm , Q θ ) of the gradient ∇ θ w p p (P, Q θ ). We now show that even in this simplest of settings, this estimate is biased, and we exhibit a lower bound on the bias for any value of m. Hence the Wasserstein metric does not have property (U). More worrisome still, we show that the minimum of the expected empirical Wasserstein loss θ → E Xm w p p ( Pm , Q θ ) is not the minimum of the Wasserstein loss θ → w p p (P, Q θ ). We then conclude that minimizing the sample Wasserstein loss by stochastic gradient descent may in general fail to converge to the minimum of the true loss. Theorem 1. Let Pm = 1 m m i=1 δ Xi be the empirical distribution derived from m independent samples X m = X 1 , . . . , X m drawn from a Bernoulli distribution P . Then for all 1 ≤ p < ∞,

• Non-vanishing minimax bias of the sample gradient. For any m ≥ 1 there exists a pair of Bernoulli distributions P , Q θ for which

E Xm∼P ∇ θ w p p ( Pm , Q θ ) -∇ θ w p p (P, Q θ ) ≥ 2e -2 ;

• Wrong minimum of the sample Wasserstein loss. The minimum of the expected sample loss θ = arg min θ E Xm w p p ( Pm , Q θ ) is in general different from the minimum of the true Wasserstein loss θ * = arg min θ w p p (P, Q θ ). • Deterministic solutions to stochastic problems. For any m ≥ 1, there exists a distribution P with nonzero entropy whose sample loss is minimized by a distribution Q θ with zero entropy.

Taken as a whole, Theorem 1 states that we cannot in general minimize the Wasserstein loss using naive stochastic gradient descent methods. Although our result does not imply the lack of a stochastic optimization procedure for this loss,2 we believe our result to be cause for concern. We leave as an open question whether an unbiased optimization procedure exists and is practical.

WASSERSTEIN BIAS IN THE LITERATURE:

Our result is surprising given the prevalence of the Wasserstein metric in empirical studies. We hypothesize that this bias exists in published results and is an underlying cause of learning instability and poor convergence often remedied to by heuristic means. For example, Frogner et al. (2015) and Montavon et al. (2016) reported the need for a mixed KL-Wasserstein loss to obtain good empirical results, with the latter explicitly discussing the issue of wrong minima when using Wasserstein gradients.

We remark that our result also applies to the dual (2), since the losses are the same. This dual was recently considered by Arjovsky et al. (2017) as an alternative loss to the primal (1). The adversarial procedure proposed by the authors is a two time-scale process which first maximizes (2) w.r.t f ∈ F ∞ using m samples, then takes a single stochastic gradient step w.r.t. θ. Interestingly, this approach does seem to provide unbiased gradients as m → ∞. However, the cost of a single gradient is now significantly higher, and for a fixed m we conjecture that the minimax bias remains.

4. THE CRAMÉR DISTANCE:

We are now ready to describe an alternative to the Wasserstein metric, the Cramér distance (Székely, 2002; Rizzo & Székely, 2016) . As we shall see, the Cramér distance has the same appealing properties as the Wasserstein metric, but also provides us with unbiased sample gradients. As a result, we believe this underappreciated distance is an appealing alternative to the Wasserstein metric for many machine learning applications.

4.1. DEFINITION AND ANALYSIS:

Recall that for two distributions P and Q over R, their (cumulative) distribution functions are respectively F P and F Q . The (squared) Cramér distance between P and Q is 10) is significantly more distant than the two others (0, 1). Rest. Distributions minimizing the divergences discussed in this paper, under the constraint Q(1) = Q(10). Both Wasserstein metric and Cramér distance underemphasize Q(0) to better match the cumulative distribution function. The sample Wasserstein loss result is for m = 1.

l 2 2 (P, Q) := ∞ -∞ (F P (x) -F Q (x)) 2 dx.

The Cramér distance is a Bregman divergence, and is a member of the l p family of divergences

l p (P, Q) := ∞ -∞ |F P (x) -F Q (x)| p dx 1/p .

The l p and Wasserstein metrics are identical at p = 1, but are otherwise distinct. As the following theorem shows, the Cramér distance possesses unique properties. Theorem 2. Consider two random variables X, Y , a random variable A independent of X, Y , and a real value c > 0. Then for 1 ≤ p ≤ ∞,

(I) l p (A + X, A + Y ) ≤ l p (X, Y ) (S) l p (cX, cY ) ≤ |c| 1/p l p (X, Y ).

Furthermore, the Cramér distance has unbiased sample gradients. That is, given

X m := X 1 , . . . , X m drawn from a distribution P , the empirical distribution Pm := 1 m m i=1 δ Xi , and a distribution Q θ , E Xm∼P ∇ θ l 2 2 ( Pm , Q θ ) = ∇ θ l 2 2 (P, Q θ ),

and of all the l p distances, only the Cramér (p = 2) has this property.

We conclude that the Cramér distance enjoys both the benefits of the Wasserstein metric and the SGD-friendliness of the KL divergence. Given the close similarity of the Wasserstein and l p metrics, it is truly remarkable that only the Cramér distance has unbiased sample gradients.

4.2. COMPARISON TO THE 1-WASSERSTEIN METRIC:

To illustrate how the Cramér distance compares to the 1-Wasserstein metric, we consider modelling the discrete distribution P depicted in Figure 1 (left). Since the trade-offs between metrics are only apparent when using an approximate model, we use an underparametrized discrete distribution Q θ which assigns the same probability to x = 1 and x = 10. That is,

Q θ (0) := Q θ {x = 0} = 1 1 + 2e θ Q θ (1) = Q θ (10) = e θ 1 + 2e θ .

Figure 1 depicts the distributions minimizing the various divergences under this parametrization. In particular, the Cramér solution is relatively close to the 1-Wasserstein solution. Furthermore, the minimizer of the sample Wasserstein loss (m = 1) clearly provides a bad solution (most of the mass is on 0). Note that, as implied by Theorem 1, the bias shown here would arise even if the distribution could be exactly represented.

To further show the impact of the Wasserstein bias we used gradient descent to minimize either the true or sample losses with a fixed step-size (α = 0.001). In the stochastic setting, at each step we construct the empirical distribution Pm from m samples (a Dirac when m = 1), and take a gradient step. We measure the performance of each method in terms of the true 1-Wasserstein loss. small sample sizes stochastic gradient descent fails to find reasonable solutions, and for m = 1 even converges to a solution worse than the KL minimizer. This small experiment highlights the cost incurred from minimizing the sample Wasserstein loss, and shows that increasing the sample size may not be sufficient to guarantee good behaviour.

ORDINAL REGRESSION:

We next trained a neural network in an ordinal regression task using either of the three divergences. The task we consider is the Year Prediction MSD dataset (Lichman, 2013) . In this task, the model must predict the year a song was written (from 1922 to 2011) given a 90-dimensional feature representation. In our setting, this prediction takes the form of a probability distribution. We measure each method's performance on the test set (Figure 2 ) in two ways: root mean squared error (RMSE) -the metric minimized by Hernández-Lobato & Adams (2015) -and the sample Wasserstein loss. Full details on the experiment may be found in the appendix.

The results show that minimizing the sample Wasserstein loss results in significantly worse performance. By contrast, minimizing the Cramér distance yields the lowest RMSE and Wasserstein loss, confirming the practical importance of having unbiased sample gradients. Naturally, minimizing for one loss trades off performance with respect to the others, and minimizing the Cramér distance results in slightly higher negative log likelihood than when minimizing the KL divergence (Figure 7 in appendix). We conclude that, in the context of ordinal regression where outcome similarity plays an important role, the Cramér distance should be preferred over either KL or the Wasserstein metric.

5. MULTIVARIATE DISTRIBUTIONS:

The energy distance (Székely, 2002) is a natural extension of the Cramér distance to the multivariate case. Let P, Q be probability distributions over R d and let X, X and Y, Y be independent random variables distributed according to P and Q, respectively. The energy distance (sometimes called the squared energy distance, see e.g. Rizzo & Székely, 2016) is

E(P, Q) := E(X, Y ) := 2 E X -Y 2 -E X -X 2 -E Y -Y 2 .

(3)

Székely showed that, in the univariate case, l 2 2 (P, Q) = 1 2 E(P, Q). Interestingly enough, the energy distance can also be written in terms of a difference of expectations. For

EQUATION

The energy distance is closely related to the distances known as maximum mean discrepancies (MMDs; Gretton et al., 2012) ; in particular, Sejdinovic et al. (2013) showed that the energy distance is equivalent to the squared MMD with kernel k(x, y) = x 2 + y 2 -xy 2 . Finally, we remark that E also possesses properties (I), (S), and (U) (proof in the appendix). (Liu et al., 2015) .

Algorithm 1: Cramér GAN Losses.

Parameter. Gradient penalty coefficient λ. Sample x r ∼ P , x g , x g ∼ Q, ∼ Uniform(0, 1).

Interpolate real and generated samples:

x = x r + (1 -)x g

Sample generator loss ( 12):

Lg = h(x r ) -h(x g ) 2 + h(x r ) -h(x g ) 2 -h(x g ) -h(x g

) 2 Sample surrogate generator loss (13) and critic loss:

Ls (u, v) = h(x r ) -h(v) 2 -h(x r ) 2 -h(u) -h(v) 2 + h(u) 2 L s = 1 2 Ls (x g , x g ) + Ls (x g , x g ) f (x) = h(x) -h(x g ) 2 -h(x) -h(x r ) 2 L critic = -L s + λ( ∇ xf (x) 2 -1) 2 Generator Steps Wasserstein Distance WGAN-GP WGAN-GP Cramer GAN Cramer GAN Nu = 1 Nu = 5 Nu = 1 Nu = 5

Figure 4 : Approximate Wasserstein distances between CelebA test set and the generators. N u is the number critic updates per generator update.

5.1. CRAMÉR GAN:

We now consider the Generative Adversarial Networks (GAN) framework (Goodfellow et al., 2014) , in particular issues arising in the Wasserstein GAN (Arjovsky et al., 2017) , and propose a better GAN based on the Cramér distance. A GAN is composed of a generative model Q (in our experiments, over images), called the generator, a target source P , and a trainable loss function called a discriminator or critic. GANs are particularly interesting because we can establish a direct comparison between the two distances. Our choice of name reflects this fact, and we prefer Cramér GAN to the perhaps more technically correct, but less palatable Energy Distance GAN. In theory, the Wasserstein GAN algorithm requires training the critic until convergence, but this is rarely achievable: we would require a critic that is a very powerful network to approximate the Wasserstein distance well (Arora et al., 2017) . Simultaneously, training this critic to convergence would overfit the empirical distribution of the training set, which is undesirable.

Our proposed loss function allows for useful learning with imperfect critics by combining the energy distance with a transformation function h : R d → R k , where d is the input dimensionality and k = 256 in our experiments. The generator then seeks to minimize the energy distance of the transformed variables E(h(X), h(Y )), where X is a real sample and Y is a generated sample. The critic itself seeks to maximize this same distance by changing the parameters of h, subject to a soft constraint (the gradient penalty used by Gulrajani et al., 2017) . Specifically, the critic maximizes a surrogate loss whose gradient can be estimated from a single real sample. The Cramér GAN losses are summarized in Algorithm 1, with additional design choices detailed in Appendix C.

We note that MMDs such as the energy distance have in the last year become an appealing tool for training GANs. Among others, the squared MMD is used within Generative Moment Matching Networks (Li et al., 2015; Dziugaite et al., 2015) ; Bouchacourt et al. (2016) trained a model to minimize the energy distance for hand pose estimation. Our use of the tranformation h(x) reflects our anecdotal finding that the direct minimization of the energy distance over raw images does not work well (see Figure 10 in appendix). Similar findings can be found in the work of Mroueh et al. (2017) and the independently developed MMD GAN (Li et al., 2017) , which additionally uses an auto-encoder loss to make the transformation injective.

The Cramér GAN we present here complements our comparison of the Wasserstein and Cramér distance from previous sections. At the same time, our experiments also provide novel GAN-related contributions, including the ability to perform conditional modelling using a surrogate generator loss, which lets us train the critic even when only one independent sample from P is available. We note also that in our experiments, xy 2 distances were more stable than distances generated by Gaussian or Laplacian kernels.

5.2. CRAMÉR GAN EXPERIMENTS:

We now show that, compared to the improved Wasserstein GAN (WGAN-GP) of Gulrajani et al. (2017) , the Cramér GAN leads to more stable learning and increased diversity in the generated samples. In both cases we train generative models that predict the right half of an image given the left half; samples from unconditional models are provided in the appendix (Figure 10 ). The dataset we use here is the CelebA 64 × 64 dataset (Liu et al., 2015) of celebrity faces.

Increased diversity. In our first experiment, we compare the qualitative diversity of completed faces by showing three sample completions generated by either model given the left half of a validation set image (Figure 3 ). We observe that the completions produced by WGAN-GP are almost deterministic. Our findings echo those of Isola et al. (2016) , who observed that "the generator simply learned to ignore the noise." By contrast, the completions produced by Cramér GAN are fairly diverse, including different hairstyles, accessories, and backgrounds. We view this lack of diversity in WGAN-GP as undesirable given that the main requirement of a generative model is that it should provide a variety of outputs.

Theorem 1 provides a clue as to what may be happening here. We know that minimizing the sample Wasserstein loss will find the wrong minimum. In particular, when the target distribution has low entropy, the sample Wasserstein minimizer may actually be a deterministic distribution. But a good generative model of images must lie in this "almost deterministic" regime, since the space of natural images makes up but a fraction of all possible pixel combinations and hence there is little perpixel entropy. We hypothesize that the increased diversity in the Cramér GAN comes exactly from learning these almost deterministic predictions.

More stable learning. In a second experiment, we varied the number of critic updates (N u ) per generator update. To compare performance between the two architectures, we measured the loss computed by an independent WGAN-GP critic trained on the validation set, following a similar evaluation previously done by Danihelka et al. (2017) . Figure 4 shows the independent Wasserstein critic distance between each generator and the test set during the course of training. Echoing our results with the toy experiment and ordinal regression, the plot shows that when a single critic update is used, WGAN-GP performs particularly poorly. We note that additional critic updates also improve Cramér GAN. This indicates that it is helpful to keep adapting the h(x) transformation.

6. CONCLUSION:

There are many situations in which the KL divergence, which is commonly used as a loss function in machine learning, is not suitable. The desirable alternatives, as we have explored, are the divergences that are ideal and allow for unbiased estimators: they allow geometric information to be incorporated into the optimization problem; because they are scale-sensitive and sum-invariant, they possess the convergence properties we require for efficient learning; and the correctness of their sample gradients means we can deploy them in large-scale optimization problems. Among open questions, we mention deriving an unbiased estimator that minimizes the Wasserstein distance, and variance analysis and reduction of the Cramér distance gradient estimate.

A PROOFS A.1 PROPERTIES OF A DIVERGENCE:

Proof (Proposition 1 and 2). The statement regarding (U) for the KL divergence is well-known, and forms the basis of most stochastic gradient algorithms for classification. Chung & Sobel (1987) have shown that the total variation does not have property (S); by Pinsker's inequality, it follows that the same holds for the KL divergence. A proof of (I) and (S) for the Wasserstein metric is given by Bickel & Freedman (1981) , while the lack of (U) is shown in the proof of Theorem 1.

A.2 BIASED ESTIMATE:

Proof (Theorem 1). Minimax bias: Consider P = B(θ * ), a Bernoulli distribution of parameter θ * and Q θ = B(θ) a Bernoulli of parameter θ. The empirical distribution Pm is a Bernoulli with parameter θ := 1 m m i=1 X i . Note that with P and Q θ both Bernoulli distributions, the p th powers of the p-Wasserstein metrics are equal, i.e. w 1 (P, Q θ ) = w p p (P, Q θ ). This gives us an easy way to prove the stronger result that all p-Wasserstein metrics have biased sample gradients. The gradient of the loss w p p (P,

Q θ ) is, for θ = θ * , g := ∇w p p (P, Q θ ) = ∇ θ * -θ = sgn(θ -θ * ),

and similarly, the gradient of the sample loss is, for θ = θ,

EQUATION

which is different from g for any θ * ∈ (0, 1). In particular for m = 1, E P ĝ = 1 -2θ * does not depend on θ, thus a gradient descent using a one-sample gradient estimate has no chance of minimizing the Wasserstein loss as it will converge to either 1 or 0 instead of θ * . Now observe that for m ≥ 2, and any θ > m-1 m , Pr{ θ < θ} = Pr{∃i s.t.

X i = 0} = 1 -(θ * ) m ,

and therefore

E ĝ = 1 -2(θ * ) m . Taking θ * = m-1 m , we find that g -E ĝ = 1 -[1 -2(θ * ) m ] = 2 1 - 1 m m ≥ 2e -2 .

Thus for any m, there exists P = B(θ * ) and Q θ = B(θ) with θ * = m-1 m < θ < 1 such that the bias g -E ĝ is lower-bounded by a numerical constant. Thus the minimax bias does not vanish with the number of samples m.

Notice that a similar argument holds for θ * and θ being close to 0. In both situations where θ * is close to 0 or 1, the bias is non vanishing when |θ * -θ| is of order 1 m . However this is even worse when θ * is away from the boundaries. For example chosing θ * = 1 2 , we can prove that the bias is non vanishing even when |θ * -θ| is (only) of order 1 √ m .

Indeed, using the anti-concentration result of Veraar (2010) (Proposition 2), we have that for a sequence Y 1 , . . . , Y m of Rademacher random variables (i.e. +/ -1 with equal probability),

Pr 1 n m i=1 Y i ≥ ≥ (1 -m 2 ) 2 /3.

This means that for samples X 1 , . . . , X m drawn from a Bernoulli B(θ * = 1 2 ) (i.e., Y i = 2X i -1 are Rademacher), we have thus for 1/2 = θ * < θ < θ * + 1/ √ 8m we have the following lower bound on the bias:

Pr θ ≥ θ * + /2 ≥ (1 -m 2 ) 2 /3,

g -E ĝ = 2 Pr θ ≥ θ ≥ 1/6.

Thus the bias is lower-bounded by a constant (independent of m) when θ * = 1 2 and |θ * -θ| = O(1/ √ m).

Wrong minimum: From (5), we deduce that a stochastic gradient descent algorithm based on the sample Wasserstein gradient will converge to a θ such that Pr{ θ < θ} = 1 2 , i.e., θ is the median of the distribution over θ, whereas θ * is the mean of that distribution. Since θ follows a (normalized) binomial distribution with parameters m and θ * , we know that the median θ and the mean θ * do not necessarily coincide, and can actually be as far as 1 2m -away from each other. For example for any odd m and any θ * ∈ 1 2 , 1 2 -1 2m the median is θ * -1 2m . It follows that the minimum of the expected sample Wasserstein loss (the fixed point of the stochastic gradient descent using the sample Wasserstein gradient) is different from the minimum of the true Wasserstein loss:

arg min θ E[w p p ( Pm , Q θ )] = arg min θ [w p p (P, Q θ )].

This is illustrated in Figure 5 .

Notice that the fact that the minima of these losses differ is worrisome as it means that minimizing the sample Wasserstein loss using (finite) samples will not converge to the correct solution.

Deterministic solutions: Consider the specific case where (1/2) 1/n < θ * < 1 (illustrated in the right plot of Figure 5 ). Then the expected sample gradient ∇ E[w p p ( Pm , Q θ * )] = Eĝ = 1-2(θ * ) n < 0 for any θ, so a gradient descent algorithm will converge to 1 instead of θ * . Notice that a symmetric argument applies for θ * close to 0.

In this simple example, minimizing the sample Wasserstein loss may lead to degenerate solutions (i.e., deterministic) when our target distributions have low (but not zero) entropy.

A.3 CONSISTENCY OF THE SAMPLE 1-WASSERSTEIN GRADIENT:

We provide an additional result here showing that the sample 1-Wasserstein gradient converges to the true gradient as m → ∞.

Theorem 3. Let P and Q θ be probability distributions, with Q θ parametrized by θ. Assume that the set x ∈ X, such that F P (x) = F Q θ (x) has measure zero, and that for any x ∈ X, the map θ → F Qθ (x) is differentiable in a neighborhood V(θ) of θ with a uniformly bounded derivative (for θ ∈ V(θ) and x ∈ X). Let Pm = 1 m m i=1 δ Xi be the empirical distribution derived from m independent samples X 1 , . . . , X m drawn from P . Then

lim m→∞ ∇w 1 ( Pm , Q θ ) = ∇w 1 (P, Q θ ), almost surely.

We note that the measure requirement is strictly to keep the proof simple, and does not subtract from the generality of the result.

Proof. Let ∇ := ∇ θ . Since p = 1 the Wasserstein distance w 1 (P, Q) measures the area between the curves defined by the distribution function of P and Q, thus w 1 (P,

Q) = l 1 (P, Q) = F P (x) - F Q (x) dx and ∇w 1 (P, Q θ ) = lim ∆→0 w 1 (P, Q θ+∆ ) -w 1 (P, Q θ ) ∆ = lim ∆→0 1 ∆ F P (x) -F Q θ+∆ (x) -F P (x) -F Q θ (x) dx.

Now since we have assumed that for any x ∈ X, the map θ → F Q θ (x) is differentiable in a neighborhood V(θ) of θ and its derivative uniformly (over V(θ) and x) bounded by M , we have

1 ∆ F P (x) -F Q θ+∆ (x) -F P (x) -F Q θ (x) ≤ 1 ∆ F Q θ+∆ (x) -F Q θ (x) ≤ M.

Thus the dominated convergence theorem applies and

∇w 1 (P, Q θ ) = lim ∆→0 1 ∆ F P (x) -F Q θ+∆ (x) -F P (x) -F Q θ (x) dx = ∇ F P (x) -F Q θ (x) dx = sgn F P (x) -F Q θ (x) ∇F Q θ (x)dx,

since we have assumed that the set of x ∈ X such that F P (x) = F Q θ (x) has measure zero. Now, using the same argument for w 1 ( Pm , Q θ ) we deduce that

∇w 1 ( Pm , Q θ ) = lim ∆→0 1 ∆ F Pm (x) -F Q θ+∆ (x) -F Pm (x) -F Q θ (x) A(x)

dx.

Let us decompose this integral over X as the sum of two integrals, one over X \ Ω m and the other one over Ω m , where

Ω m = x ∈ X, F Pm (x) = F Q θ (x) . We have X\Ωm A(x)dx = X\Ωm sgn F Pm (x) -F Q θ (x) ∇F Q θ (x)dx,

and

Ωm A(x)dx ≤ Ωm lim ∆→0 1 ∆ F Q θ+∆ (x) -F Q θ (x) dx ≤ M |Ω m |.

Now from the strong law of large numbers, we have that for any x, the empirical cumulative distribution function F Pm (x) converges to the cumulative distribution F P (x) almost surely. We deduce that Ω m converges to the set x, F P (x) = F Q θ (x) which has measure zero, thus |Ω m | → 0 and

lim m→∞ ∇w 1 ( Pm , Q θ ) = lim m→∞ X sgn F Pm (x) -F Q θ (x) ∇F Q θ (x)dx. Now, since |∇F Q θ (x)| ≤ M

, we can use once more the dominated convergence theorem to deduce that

lim m→∞ ∇w 1 ( Pm , Q θ ) = X lim m→∞ sgn F Pm (x) -F Q θ (x) ∇F Q θ (x)dx = X sgn F P (x) -F Q θ (x) ∇F Q θ (x)dx = ∇w 1 (P, Q θ ).

The following lemma will be useful in proving that the Cramér distance has property (U).

Lemma 1. Let X m := X 1 , . . . , X m be independent samples from P , and let Pm :=

1 m i δ Xi . Then E Xm∼P F Pm (x) = F P (x).

Proof. Because the X i 's are independent,

F Pm (x) = x -∞ Pm (dx) = 1 m m i=1 I [X i ≤ x] .

Now, taking the expectation w.r.t. X m ,

E Xm∼P F Pm (x) = E Xm∼P 1 m m i=1 I [X i ≤ x] = 1 m m i=1 E Xi∼P I [X i ≤ x] = 1 m m i=1 Pr{X i ≤ x} = F P (x),

since the X i are identically distributed according to P .

Proof (Theorem 2). Like the Wasserstein metrics, the l p metrics have dual forms as integral probability metrics (see Dedecker & Merlevède, 2007 , for a proof):

EQUATION

where F q := {f : f is absolutely continuous, df dx q ≤ 1} and q is the conjugate exponent of p, i.e. p -1 + q -1 = 1.3 We will use this dual form below.

We will prove that l p has properties (I) and (S) for p ∈ [1, ∞); the case p = ∞ follows by a similar argument. Begin by observing that

F cX (x) = P r{cX ≤ x} = P r X ≤ x c = F X x c .

Then we may rewrite l p p (cX, cY ) as

l p p (cX, cY ) = ∞ -∞ F X x c -F Y x c p dx (a) = c ∞ -∞ F X (z) -F Y (z) p dz,

where (a) uses a change of variables z = x/c. Taking both sides to the power 1/p proves that the l p metric possesses property (S) of order 1/p. For (I), we use the IPM formulation (6):

l p (A + X, A + Y ) = sup f ∈Fq E A+X f (x) -E A+Y f (y) (a) = sup f ∈Fq E A E X f (x + a) -E A E Y f (y + a) = sup f ∈Fq E A E X f (x + a) -E Y f (y + a) (b) ≤ E A sup f ∈Fq E X f (x + a) -E Y f (y + a) ,

where (a) is by independence of A and X, Y , and (b) is by Jensen's inequality. Next, recall that F q is the set of absolutely continuous functions whose derivative has bounded L q norm. Hence if f ∈ F q , then also for all a the translate g a (x) := f (x + a) is also in F q . Therefore,

l p (A + X, A + Y ) ≤ E A sup f ∈Fq E X f (x + a) -E Y f (y + a) = E A sup g∈Fq E X g(x) -E Y g(y) = sup g∈Fq E X g(x) -E Y g(y) = l p (X, Y ).

Now, to prove (U).:

Here we make use of the introductory requirement that "all expectations under consideration are finite." Specifically, we require that the mean under P , E x∼P [x], is well-defined and finite, and similarly for Q θ . In this case,

EQUATION

This mild requirement guarantees that the tails of the distribution function F P are light enough to avoid infinite Cramér distances and expected gradients (a similar condition was set by Dedecker & Merlevède (2007) ). Now, by definition,

∇ θ l 2 2 (P, Q θ ) = ∇ θ ∞ -∞ (F Q θ (x) -F P (x)) 2 dx (a) = ∞ -∞ ∇ θ (F Q θ (x) -F P (x)) 2 dx = ∞ -∞ 2 (F Q θ (x) -F P (x)) ∇ θ F Q θ (x)dx (b) = ∞ -∞ 2 F Q θ (x) -E Xm F Pm (x) ∇ θ F Q θ (x)dx = ∞ -∞ 2 E Xm F Q θ (x) -F Pm (x) ∇ θ F Q θ (x)dx (c) = E Xm ∞ -∞ 2 F Q θ (x) -F Pm (x) ∇ θ F Q θ (x)dx = E Xm ∇ θ l 2 2 ( Pm , Q θ )

, where (a) follows from the hypothesis (7) (the convergence of the squares follows from the convergence of the ordinary values), (b) follows from Lemma 1 and (c) follows from Fubini's theorem, again invoking (7).

Finally, we prove that of all the l p p distances (1 ≤ p ≤ ∞) only the Cramér distance, l 2 2 , has the (U) property.

Without loss of generality, let us suppose P is not a Dirac, and further suppose that for any X m ∼ P , F Q θ (x) ≥ F Pm (x) everywhere. For example, when Q θ has bounded support we can take P to be a sufficiently translated version of Q θ , such that the two distributions' supports do not overlap.

We have already established that the 1-Wasserstein does not have the (U) property, and is equivalent to l p p for p = 1. We will thus assume that p > 1, and also that p < ∞, the latter being recovered through standard limit arguments. Begin with the gradient for l p p (P, Q θ ),

∇ θ l p p (P, Q θ ) = ∇ θ ∞ -∞ F Q θ (x) -F P (x) p dx (a) = p ∞ -∞ F Q θ (x) -F P (x) p-1 ∇ θ F Q θ (x)dx = p ∞ -∞ φ p (F Q θ (x) -F P (x))∇ θ F Q θ (x)dx = p ∞ -∞ φ p E Xm F Q θ (x) -F Pm (x) ∇ θ F Q θ (x)dx,

for φ p (z) = z p-1 ; in (a) we used the same argument as in Theorem 3. Now, φ p is convex on [0, ∞) when p ≥ 2 and concave on the same interval when 1 < p < 2. From Jensen's inequality we know that for a convex (concave) function φ and a random variable Z, E φ(Z) is greater than (less than) or equal to φ(E Z), with equality if and only if φ is linear or Z is deterministic. By our first assumption we have ruled out the latter. By our second assumption

F Q θ (x) ≥ F Pm (x)

, we can apply Jensen's inequality at every x to deduce that

E Xm ∇ θ l p p ( Pm , Q θ ) < ∇ θ l p p (P, Q θ ), if 1 < p < 2, E Xm ∇ θ l p p ( Pm , Q θ ) > ∇ θ l p p (P, Q θ ), if p > 2, E Xm ∇ θ l p p ( Pm , Q θ ) = ∇ θ l p p (P, Q θ ), if p = 2.

We conclude that of the l p p distances, only the Cramér distance has unbiased sample gradients.

Proposition 3. The energy distance E(P, Q) has properties (I), (S), and (U).

Proof. As before, write E(X, Y ) := E(P, Q). Recall that

E(X, Y ) = 2 E X -Y 2 -E X -X 2 -E Y -Y 2 .

Consider a random variable A independent of X and Y . First, we want to prove property (I):

E(A + X, A + Y ) ≤ E(X, Y ).

We will use Proposition 2 from Székely & Rizzo (2013) to express the energy distance in terms of characteristic functions φ X , φ Y of d-dimensional random variables X and Y :

E(X, Y ) = 1 c d R d |φ X (t) -φ Y (t)| 2 |t| d+1 dt

where

c d = π (d+1)/2 Γ( d+1 2 )

.

The proof then uses properties of characteristic functions (|φ A (t)| ≤ 1 and φ A+X (t) = φ A (t)φ X (t) for independent variables A and X) to show:

E(A + X, A + Y ) = 1 c d R d |φ A+X (t) -φ A+Y (t)| 2 |t| d+1 dt = 1 c d R d |φ A (t)φ X (t) -φ A (t)φ Y (t)| 2 |t| d+1 dt = 1 c d R d |φ X (t) -φ Y (t)| 2 |t| d+1 |φ A (t)| 2 dt ≤ 1 c d R d |φ X (t) -φ Y (t)| 2 |t| d+1 dt = E(X, Y ).

This proves (I). Next, consider a real value c > 0. We have

E(cX, cY ) = 2 E cX -cY 2 -E cX -cX 2 -E cY -cY 2 = 2c E X -Y 2 -c E X -X 2 -c E Y -Y 2 = cE(X, Y ).

This proves (S). Finally, suppose that Y is distributed according to Q θ parametrized by θ. Let X m = X 1 , . . . , X m be drawn from P , and let Pm := 1 m m i=1 δ Xi . Let X be the random variable distributed according to Pm , and X an independent copy of X. Then

E( Pm , Q θ ) = E( X, Y ) = 2 E X -Y 2 -E X -X 2 -E Y -Y 2 .

The gradient of the true loss w.r.t. θ is

EQUATION

)

Now, taking the gradient of the sample loss w.r.t. θ,

EQUATION

)

Since the second terms of the gradients match, all we need to show is that the first terms are equal, in expectation. Assuming that ∇ θ and the expectation over X m commute, we write

E Xm ∇ θ E X -Y 2 = ∇ θ E Xm E X -Y 2 = ∇ θ E Xm E E x∼ Pm x -Y 2 ,

by independence of X and Y . But now we know that the expected empirical distribution is P , that is

E Xm E x∼ Pm E x -Y 2 = E x∼P E x -Y 2 = E X -Y 2 .

It follows that the first terms of ( 8) and ( 9) are also equal, in expectation w.r.t. X m . Hence we conclude that the energy distance has property (U), that is

E Xm∼P ∇ θ E( Pm , Q θ ) = ∇ θ E(P, Q θ ).

B COMPARISON WITH THE WASSERSTEIN DISTANCE:

Figure 2 (left) provides learning curves for the toy experiment described in Section 4.2. 

B.1 ORDINAL REGRESSION:

We compare the different losses on an ordinal regression task using the Year Prediction MSD dataset from (Lichman, 2013) . The task is to predict the year of a song (taking on values from 1922 to 2011), from 90-dimensional feature representation of the song. 4 Previous work has used this dataset for benchmarking regression performance (Hernández-Lobato & Adams, 2015) , treating the target as a continuous value. Following Hernández-Lobato & Adams (2015) , we train a network with a single hidden layer with 100 units and ReLU non-linearity, using SGD with 40 passes through the training data, using the standard train-test split for this dataset (Lichman, 2013) . Unlike (Hernández-Lobato & Adams, 2015) , the network outputs a probability distribution over the years (90 possible years from 1922-2011).

We train models using either the 1-Wasserstein loss, the Cramér loss, or the KL loss, the latter of which reduces the ordinal regression problem to a classification problem. In all cases, we compare performance for three different minibatch sizes, i.e. the number of input-target pairs per gradient step. Note that the minibatch size only affects the gradient estimation, but has otherwise no direct relation to the number of samples m previously discussed, since each sample corresponds to a different input vector. We report results as a function of number of passes over the training data so that our results are comparable with previous work, but note that smaller batch sizes get more updates.

The results are shown in Figure 2 . Training using the Cramér loss results in the lowest root mean squared error (RMSE) and the final RMSE value of 8.89 is comparable to regression (Hernández-Lobato & Adams, 2015) which directly optimizes for MSE. We further observe that minimizing the Wasserstein loss trains relatively slowly and leads to significantly higher KL loss. Interestingly, larger minibatch sizes do seem to improve the performance of the Wasserstein-based method somewhat, suggesting that there might be some beneficial bias reduction from combining similar inputs. By contrast, using with the Cramér loss trains significantly faster and is more robust to choice of minibatch size. 

B.2 IMAGE MODELLING WITH PIXELCNN:

As additional supporting material, we provide here the results of experiments on learning a probabilistic generative model on images using either the 1-Wasserstein, Cramér, or KL loss. We trained a PixelCNN model (Van den Oord et al., 2016) on the CelebA 32x32 dataset (Liu et al., 2015) , which is constituted of 202,599 images of celebrity faces. At a high level, probabilistic image modelling involves defining a joint probability Q θ over the space of images. PixelCNN forms this joint probability autoregressively, by predicting each pixel using a histogram distribution conditional on a probability-respecting subset of its neighbours. This kind of modelling task is a perfect setting to study Wasserstein-type losses, as there is a natural ordering on pixel intensities. This is also a setting in which full distributions are almost never available, because each prediction is conditioned on very different context; and hence we require a loss that can be optimized from single samples.

Here the true losses are not available. Instead we report the sample Wasserstein loss, which is an upper bounds on the true loss Bellemare et al. (proof is provided by 2017). For the KL divergence we report the cross-entropy loss, as is typically done; the KL divergence itself corresponds to the expected cross-entropy loss minus the real distribution's (unknown) entropy.

Figure 8 shows, as in the toy example, that minimizing the Wasserstein distance by means of stochastic gradient fails. The Cramér distance, on the other hand, is as easily minimized as the KL and in fact achieves lower Wasserstein and Cramér loss. We note that the resulting KL loss is higher than when directly minimizing the KL, reflecting the very real trade-off of using one loss over another. We conclude that in the context of learning an autoregressive image model, the Cramér should be preferred to the Wasserstein metric.

C CRAMÉR GAN C.1 LOSS FUNCTION DETAILS:

Our critic has a special form:

f (x) = E Y ∼Q h(x) -h(Y ) 2 -E X ∼P h(x) -h(X ) 2

where Q is the generator and P is the target distribution. The critic has trainable parameters only inside the deep network used for the transformation h. From (4), we define the generator loss to be

EQUATION

as in Wasserstein GAN, except that no max f operator is present and we can obtain unbiased sample gradients. At the same time, to provide helpful gradients for the generator, we train the transformation h to maximize the generator loss. Concretely, the critic seeks to maximize the generator loss while minimizing a gradient penalty:

EQUATION

)

where GP is the gradient penalty from the original WGAN-GP algorithm (Gulrajani et al., 2017) (the penalty is given in Algorithm 1). The gradient penalty bounds the critic's outputs without using a saturating function. We chose λ = 10 from a short parameter sweep. Our training is otherwise similar to the improved training of Wasserstein GAN (Gulrajani et al., 2017) .

In the next two sections, we describe how to practically compute gradients of these losses with respect to the generator and transformation parameters, respectively.

C.2 GRADIENT ESTIMATES FOR THE GENERATOR:

Recall that the energy distance is:

E(X, Y ) = 2 E X∼P Y ∼Q X -Y 2 -E X∼P X ∼P X -X 2 -E Y ∼Q Y ∼Q Y -Y 2

If Y is generated from the standard normal noise Z ∼ N (0, 1) by a differentiable generator Y = G(Z) and the generator has an integrable gradient, we can use the reparametrization trick (Kingma & Welling, 2014) to compute the gradient with respect to the generator parameters:

∇ θ G E(X, Y ) = 2 E X∼P Z∼N (0,1) ∇ θ G X -G(Z) 2 - E Z∼N (0,1) Z ∼N (0,1) ∇ θ G G(Z) -G(Z) 2 .

We see that we only need one real sample X to estimate the gradient, because the X -X term does not depend on the generator parameters. This allows us to define a generator loss usable for situations with only one real sample (e.g., for conditional modeling):

EQUATION

C.3 GRADIENT ESTIMATES FOR THE TRANSFORMATION:

As shown in the previous section, we can obtain an unbiased gradient estimate of the generator loss (12) from three samples: two from the generator, and one from the target distribution. However, to estimate the gradient of the Cramér GAN loss with respect to the transformation parameters we need four independent samples: two from the generator and two from the target distribution. In many circumstances, for example when learning conditional densities, we do not have access to two independent target samples. We will instead define a surrogate objective for the critic. The surrogate critic will have the following form: Both generators had the same DCGAN architecture (Radford et al., 2015) .

f s (x) = E Y ∼Q h(x) -h(Y ) 2 -h(x) 2

which we use to define a surrogate loss L s (X, Y ) similar to (10):

EQUATION

The surrogate loss emulates an integral probability metric (IPM) (Müller, 1997) and can be used to train the critic. The maximization of this loss will force E h(X)-h(Y ) 2 and E h(Y )-h(Y ) 2 to be informative about the underlying distributions.

The generator can be then trained to minimize the energy distance Lg (12) of the transformed variables. It is also possible to obtain training more similar to Wasserstein GAN by training the generator to minimize the surrogate loss (13). We recommend trying both possibilities, because they were both stable and produced diverse conditional samples. The whole training procedure is summarized as Algorithm 1.

Finally, when estimating the losses in Algorithm 1, we use two independent samples x g , x g from the generator. However, in constructing the surrogate loss Ls , an asymmetry arises. We reduce variance by averaging the two losses Ls (x g , x g ) and Ls (x g , x g ).

C.4 GENERATOR ARCHITECTURE:

The generator architecture is the U-Net (Ronneberger et al., 2015) previously used for Image-to-Image translation (Isola et al., 2016) . We used no batch normalization and no dropout in the generator and in the critic. The network conditioned on the left half of the image and on extra 12 channels with Gaussian noise. We generated two independent samples for a given image to compute the Cramér GAN loss. To be computationally fair to WGAN-GP, we trained WGAN-GP with twice the minibatch size (i.e., the Cramér GAN minibatch size was 64, while the WGAN-GP minibatch size was 128).

C.5 CRITIC ARCHITECTURE:

Our h(x) transformation is a deep network with 256 outputs (more is better). The network has the traditional deep convolutional architecture (Radford et al., 2015) . We do not use batch normalization, as it would conflict with the gradient penalty.

C.6 PERFORMANCE EVALUATION:

We report the Inception score (Salimans et al., 2016) and the Fréchet Inception Distance (FID) (Heusel et al., 2017) in Figure 11 (left), which are commonly used measures of evaluation for GANs.

:

Properties (S) and (I) are called regularity and homogeneity by Zolotarev; we believe our choice of terms is more machine learning-friendly.

For example, if P has finite support keeping track of the empirical distribution suffices.

This relationship is the reason for the notation F∞ in the definition the dual of the 1-Wasserstein (2).

We refer to https://archive.ics.uci.edu/ml/datasets/YearPredictionMSD for the details.

http://download.tensorflow.org/models/image/imagenet/inception-2015-12-05.tgz

Figures and Tables:

Figure 1: Leftmost. Target distribution. One outcome (10) is significantly more distant than the two others (0, 1). Rest. Distributions minimizing the divergences discussed in this paper, under the constraint Q(1) = Q(10). Both Wasserstein metric and Cramér distance underemphasize Q(0) to better match the cumulative distribution function. The sample Wasserstein loss result is for m = 1.

Figure2(left) plots the resulting training curves in the 1-Wasserstein regime, with the KL and Cramér solutions indicated for reference. We first note that, compared to the KL solution, the Cramér solution has significantly smaller Wasserstein distance to the target distribution. Second, for

Figure 3: Generated right halves of the faces for WGAN-GP (left) and Cramér GAN (right). The given left halves are from CelebA 64x64 validation set(Liu et al., 2015).

Figure 5: Wasserstein loss (black curve) θ → |θ * -θ| versus expected sample Wasserstein loss (red curve) θ → E[| θ -θ|], for different values of m and θ * and p = 1. Left: m = 1, θ * = 0.6. A stochastic gradient using a one-sample Wasserstein gradient estimate will converge to 1 instead of θ * . Middle: m = 6, θ * = 0.6. The minimum of the expected sample Wasserstein loss is the median of θ which is here θ = 2 3 = θ * = 0.6. Right: m = 5, p = 0.9. The minimum of the expected sample Wasserstein is θ = 1 and not θ * = 0.9.

Figure 6: Wasserstein while training to minimize different loss functions (Wasserstein, KL, Cramér). Averaged over 10 random initializations. Error-bands indicate one standard deviation. Note the different y-axes.

Figure 8: Left, middle. Sample Wasserstein and cross-entropy loss curves on the CelebA validation data set. Right. Test loss at the end of training, in function of loss minimized (see text for details).

Figure 9: Generated right halves for WGAN-GP (left) and Cramér GAN (right) for left halves from the validation set of Downsampled ImageNet 64x64 (Van den Oord et al., 2016). The low diversity in WGAN-GP samples is consistent with the observations of Isola et al. (2016): "the generator simply learned to ignore the noise."

Figure 10: Left. Generated images from a generator trained to minimize the energy distance of raw images, E(X, Y ). Right. Generated images if minimizing the Cramér GAN loss, E(h(X), h(Y )).Both generators had the same DCGAN architecture(Radford et al., 2015).

Peer review (iclr2018):
 The manuscript proposes to use the Cramer distance as a measure between distributions (acting as a loss) when optimizing an objective function using stochastic gradient descent (SGD). Cramer distance is a Bregman divergence and is a member of the Lp family of divergences. Here a "distance" means a symmetric divergence measure that satisfies the relaxed triangle inequality. The motivation for using the Cramer distance is that it has unbiased sample gradients while still enjoying some other properties such as scale sensitivity and sum invariant. The authors also proof that for the Bernoulli distribution, there is a lower bound independent of the sample size for the deviation between the gradient of the Cramer distance, and the expectation of the estimated gradient of the Cramer distance. Then, the multivariate case of the Cramer distance, called the energy distance, is also briefly presented. The paper closes with some experiments on ordinal regression using neural networks and training GANs using the Cramer distance. In general, the manuscript is well written and the ideas are smoothly presented. While the manuscript gives some interesting insights, I find that the contribution could have been explained in a more broader sense, with a stronger compelling message. Some remarks and questions: 1. The KL divergence considered here is sum invariant but not scale sensitive, and has unbiased sample gradients. The authors are considering here the standard (asymmetric) KL divergence (sec. 2.1). Is it the case that losing scale sensitivity make the KL divergence insensitive to the geometry of the outcomes? or is it due to the fact the KL divergence is not symmetric? or ? 2. The main argument for the paper is that the simple sample-based estimate for the gradient using the Wasserstein metric is a biased estimate for the true gradient of the Wasserstein distance, and hence it is not favored with SGD-type algorithms. Are there any other estimators in the literature for the gradient of the Wasserstein distance? Was this issue overlooked in the literature? 3. I am not sure if a biased estimate for the gradient will lead to a ``wrong minimum** in an energy space that has infinitely many local minima. Of course one should use an unbiased estimate for the gradient whenever this is possible. However, even when this is possible, there is no guarantee that this will consistently lead to deeper and ``better** minima, and there is no guarantee as well that these deep local minima reflect meaningful results. 4. To what extent can one generalize theorem 1 to other probability distributions (continuous and discrete) and to the multivariate cases as well? 5. I also don*t think that the example given in sec. 4.2 and depicted in Fig. 1 is the best and simplest way to illustrate the benefit of Cramer distance over Wasserstein. Similarly, the experiments for the multivariate case using GANs and Neural Networks do not really deliver tangible, concrete and conclusive results. Partly, these results are very qualitative, which can be understood within the context of GANs. However, the authors could have used other models/algorithms where they can obtain concrete quantitative results (for this type of contribution). In addition, such sophisticated models (with various hyper-parameters) can mask the true benefit for the Cramer distance, and can also mask the extent of how good/poor is the sample estimate for the Wasserstein gradient.