Title: LEARNING IN IMPLICIT GENERATIVE MODELS

Abstract: Generative adversarial networks (GANs) provide an algorithmic framework for constructing generative models with several appealing properties: they do not require a likelihood function to be specified, only a generating procedure; they provide samples that are sharp and compelling; and they allow us to harness our knowledge of building highly accurate neural network classifiers. Here, we develop our understanding of GANs with the aim of forming a rich view of this growing area of machine learning-to build connections to the diverse set of statistical thinking on this topic, of which much can be gained by a mutual exchange of ideas. We frame GANs within the wider landscape of algorithms for learning in implicit generative models-models that only specify a stochastic procedure with which to generate data-and relate these ideas to modelling problems in related fields, such as econometrics and approximate Bayesian computation. We develop likelihood-free inference methods and highlight hypothesis testing as a principle for learning in implicit generative models, using which we are able to derive the objective function used by GANs, and many other related objectives. The testing viewpoint directs our focus to the general problem of density ratio estimation. There are four approaches for density ratio estimation, one of which is a solution using classifiers to distinguish real from generated data. Other approaches such as divergence minimisation and moment matching have also been explored in the GAN literature, and we synthesise these views to form an understanding in terms of the relationships between them and the wider literature, highlighting avenues for future exploration and cross-pollination.

1. IMPLICIT GENERATIVE MODELS:

It is useful to make a distinction between two types of probabilistic models: prescribed and implicit models (Diggle and Gratton, 1984) . Prescribed probabilistic models are those that provide an explicit parametric specification of the distribution of an observed random variable x, specifying a log-likelihood function log q θ (x) with parameters θ. Most models in machine learning and statistics are of this form, whether they be state-of-the-art classifiers for object recognition, complex sequence models for machine translation, or fine-grained spatio-temporal models tracking the spread of disease. Alternatively, we can specify implicit probabilistic models that define a stochastic procedure that directly generates data. Such models are the natural approach for problems in climate and weather, population genetics, and ecology, since the mechanistic understanding of such systems can be used to directly create a data simulator, and hence the model. It is exactly because implicit models are more natural for many problems that they are of interest and importance.

Implicit generative models use a latent variable z and transform it using a deterministic function G θ that maps from R m → R d using parameters θ. Such models are amongst the most fundamental of models, e.g., many of the basic methods for generating non-uniform random variates are based on simple implicit models and one-line transformations (Devroye, 2006) . In general, implicit generative models specify a valid density on the output space that forms an effective likelihood function:

EQUATION

EQUATION

where q(z) is a latent variable that provides the external source of randomness and equation ( 2) is the definition of the transformed density as the derivative of the cumulative distribution function. When the function G is well-defined, such as when the function is invertible, or has dimensions m = d with easily characterised roots, we recover the familiar rule for transformations of probability distributions.

We are interested in developing more general and flexible implicit generative models where the function G is a non-linear function with d > m, specified by deep networks. The integral (2) is intractable in this case: we will be unable to determine the set {G θ (z) ≤ x}, the integral will often be unknown even when the integration regions are known and, the derivative is high-dimensional and difficult to compute. Intractability is also a challenge for prescribed models, but the lack of a likelihood term significantly reduces the tools available for learning. In implicit models, this difficulty motivates the need for methods that side-step the intractability of the likelihood (2), or are likelihood-free.

Generative adversarial networks (GANs) (Goodfellow et al., 2014) provide a solution for exactly this type of problem. This connection makes it possible to export GAN-type solutions to other areas, and to import new approaches from areas addressing the same research problem. GANs specify an algorithmic framework for learning in implicit generative models, also referred to as generator networks or generative neural samplers or differentiable simulators in the literature. The framework relies on a learning principle based on discriminating real from generated data, which we shall show instantiates a core principle of likelihood-free inference, that of hypothesis and two-sample testing.

Note on notation. We denote data by the random variable x, the (unknown) true data density by p * (x), our (intractable) model density by q θ (x). q(z) is a density over latent variables z. Parameters of the model are θ, and parameters of the ratio and discriminator functions are φ.

2. HYPOTHESIS TESTING AND DENSITY RATIOS:



2.1. LIKELIHOOD-FREE INFERENCE:

Without a likelihood function, many of the widely-used tools for inference and parameter learning become unavailable. But there are tools that remain, including the method-of-moments (Hall, 2005) , the empirical likelihood (Owen, 1988) , Monte Carlo sampling (Marin et al., 2012) , and mean-shift estimation (Fukunaga and Hostetler, 1975) . Since we can easily draw samples from the model, we can use any method that compares two sets of samples-one from the true data distribution and one from the model distribution-to drive learning. This is a process of density estimation-by-comparison that tests the hypothesis that the true data distribution p * (x) and our model distribution q(x) are equal, using the density ratio function r(x) = p * (x) /q(x). The density ratio provides information about the departure of our model distribution from the true distribution, and our aim is to see this ratio be one.

The density ratio r(x) is the core quantity for hypothesis testing, motivated by either the Neyman-Pearson lemma or the Bayesian posterior evidence, appearing as the likelihood ratio or the Bayes factor (Kass and Raftery, 1995) , respectively. This is the focus of our approach for likelihood-free inference: estimating density ratios and using them as the driving principle for learning in implicit generative models.

The direct approach of computing the density ratio by first computing the individual marginals is not possible with implicit models. By directly estimating the ratio and exploiting knowledge of the probabilities involved, it will turn out that computing the ratio can be a much easier problem than computing the marginal likelihoods, and is what will make this approach appealing. There are four general approaches to consider (Sugiyama et al., 2012a) : 1) class-probability estimation, 2) divergence minimisation, 3) ratio matching, and 4) moment matching. These are highly developed research areas in their own right, but their role in providing a learning principle for density estimation is under-appreciated and opens up an exciting range of approaches for learning in implicit generative models. Figure 1 summarises these approaches by showing pathways available for learning, which follow from the choice of inference driven by hypothesis testing and comparison.

2.2. CLASS PROBABILITY ESTIMATION:

The density ratio can be computed by building a classifier to distinguish observed data from that generated by the model. This is the most popular approach for density ratio estimation and the first port of call for learning in implicit models. Hastie et al. (2013) call this approach unsupervised-assupervised learning, Qin (1998) explore this for analysis of case-control in observational studies, both Neal (2008) and Cranmer et al. (2015) explore this approach for high-energy physics applications, Gutmann and Hyvärinen (2012) exploit it for learning un-normalised models, Lopez-Paz and Oquab

Density Estimation by Comparison:

Density Difference

r = p ⇤ q ✓ Density Ratio r = p ⇤ q ✓ f-Divergence Class Probability Estimation Bregman Divergence Moment Matching Bf [r ⇤ kr] f (u) = u log u (u + 1) log(u + 1)

Mixtures with identical moments L(✓, )

Max Mean Discrepency:

H 0 : p ⇤ = q ✓ vs. p ⇤ 6 = q ✓ Figure 1: Summary of approaches for learning in implicit models. We define a joint function L(φ, θ) and alternate between optimising the loss w.r.t. φ and θ.

(2016) for causal discovery, and Goodfellow et al. (2014) for learning in implicit generative models specified by neural networks.

We denote the domain of our data by X ⊂ R d . The true data distribution has a density p * (x) and our model has density q θ (x), both defined on X . We also have access to a set of n samples

X p = {x (p) 1 , . . . , x (p)

n } from the true data distribution, and a set of n samples X q = {x (q) 1 , . . . , x (q) n } from our model. We introduce a random variable y, and assign a label y = 1 to all samples in X p and y = 0 to all samples in X q . We can now represent p * (x) = p(x|y = 1) and q θ (x) = p(x|y = 0). By application of Bayes' rule, we can compute the ratio r(x) as:

EQUATION

which indicates that the problem of density ratio estimation is equivalent to that of class probability estimation, since the problem is reduced to computing the probability p(y = 1|x). We assume that the marginal probability over classes is p(y = 1) = π, which allows the relative proportion of data from the two classes to be adjusted if they are imbalanced; in most formulations π = ½ for the balanced case, and in imbalanced cases 1-π π ≈ n /n. Our task is now to specify a scoring function, or discriminator, D(x; φ) = p(y = 1|x): a function bounded in [0,1] with parameters φ that computes the probability of data belonging to the positive (real data) class. This discriminator is related to the density ratio through the mapping D = r /(r + 1); r = D /(1 -D). Conveniently, we can use our knowledge of building classifiers and specify these functions using deep neural networks. Given this scoring function, we must specify a proper scoring rule (Gneiting and Raftery, 2007; Buja et al., 2005) for binary discrimination to allow for parameter learning, such as those in Table 1 . A natural choice is to use the Bernoulli (logarithmic) loss:

L(φ, θ) = E p(x|y)p(y) [-y log D(x; φ) -(1 -y) log(1 -D(x; φ))] (4) = πE p * (x) [-log D(x; φ)] + (1 -π)E q θ (x) [-log(1 -D(x; φ))].

(5)

Since we know the underlying generative process for q θ (x), using a change of variables, we can express the loss in term of an expectation over the latent variable z and the generative model G(z; θ): 

EQUATION

D := D(x; φ)) Bernoulli loss πE p * (x) [-log D] + (1 -π)E q θ (x) [-log(1 -D)] Brier score πE p * (x) [(1 -D) 2 ] + (1 -π)E q θ (x) [D 2 ] Exponential loss πE p * (x) 1-D D 1 2 + (1 -π)E q θ (x) D 1-D 1 2 Misclassification πE p * (x) [I[D ≤ 0.5]] + (1 -π)E q θ (x) [I[D > 0.5]] Hinge loss πE p * (x) max 0, 1 -log D 1-D + (1 -π)E q θ (x) max 0, 1 + log D 1-D Spherical πE p * (x) [-αD] + (1 -π)E q θ (x) [-α(1 -D)] ; α = (1 -2D + 2D 2 ) -1 /2

The final form of this objective ( 6) is exactly that used in generative adversarial networks (GANs) (Goodfellow et al., 2014) . In practice, the expectations are computed by Monte Carlo integration using samples from p * and q θ . Equation ( 6) allows us to specify a bi-level optimisation (Colson et al., 2007) by forming a ratio loss and a generative loss, using which we perform an alternating optimisation. Our convention throughout the paper will be to always form the ratio loss by extracting all terms in L related to the ratio function parameters φ, and minimise the resulting objective. For the generative loss, we will similarly extract all terms related to the model parameters θ, flip the sign, and minimise the resulting objective. For equation ( 6), the bi-level optimisation is:

EQUATION

Generative loss:

EQUATION

The ratio loss is minimised since it acts as a surrogate negative log-likelihood; the generative loss is minimised since we wish to minimise the probability of the negative (generated-data) class. We explicitly write out these two stages to emphasise that the objectives used are separable. While we can derive the generative loss from the ratio loss as we have done, any generative loss that drives q θ to p * , such as minimising the widely-used E q(z) [-log D(G(z; θ))] (Goodfellow et al., 2014; Nowozin et al., 2016) or (Sønderby et al., 2016) , is possible. Any scoring rule from Table 1 can be used to give a loss function for optimisation. These rules are amenable to stochastic approximation and alternating optimisation, as described by Goodfellow et al. (2014) , along with many of the insights for optimisation that have been developed since (Salimans et al., 2016; Radford et al., 2015; Zhao et al., 2016; Sønderby et al., 2016) . The Bernoulli loss can be criticised in a number of ways and makes other scoring rules interesting to explore. The Brier loss provides a similar decision rule, and its use in calibrated regression makes it appealing; the motivations behind many of these scoring rules are discussed in (Gneiting and Raftery, 2007) . Finally, while we have focussed on the two sample hypothesis test, we believe it could be advantageous to extend this reasoning to the case of multiple testing, where we simultaneously test several sets of data (Bickel et al., 2008) .

E q(z) [-log D(G(z;θ)) 1-D(G(z;θ)) ] = E q(z) [-log r(G(z; θ))

The advantage of using a proper scoring rule is that the global optimum is achieved iff q θ = p * (cf. the proof in (Goodfellow et al., 2014) for the Bernoulli loss); however there are no convergence guarantees since the optimisation is non-convex. Goodfellow et al. (2014) discuss the relationship to maximum likelihood estimation, which minimises the divergence KL[p * q], and show that the GAN objective with Bernoulli loss is instead related to the Jensen Shannon divergence JS [p * ||q]. In the objective (6), π denotes the marginal probability of the positive class; however several authors have proposed choosing π depending on the problem. In particular, Huszár (2015) showed that varying π is related to optimizing a generalised Jensen-Shannon divergence JS π [p * ||q]. Creswell and Bharath (2016) presented results showing that different values of π are desirable, depending on whether we wish to fit one of the modes (a 'high precision, low recall' task such as generation) or explain all of the modes (a 'high recall, low precision' task such as retrieval).

2.3. DIVERGENCE MINIMISATION:

A second approach to testing is to use the divergence between the true density p * and our model q, and use this as an objective to drive learning of the generative model. A natural class of divergences to use are the f -divergences (or Ali-Silvey (Ali and Silvey, 1966) or Csiszar's φ-divergence (Csiszàr, 1967) ) since they are fundamentally linked to the problem of two-sample hypothesis testing (Liese and Vajda, 2008) : f -divergences represent an integrated Bayes risk since they are an expectation of the density ratio. Nowozin et al. (2016) develop f -GANs using this view. The f -divergences contain the KL divergence as a special case and are equipped with an exploitable variational formulation:

D f [p * (x) q θ (x)] = q θ (x)f p * (x) q θ (x) dx = E q θ (x) [f (r(x))] ≥ sup t E p * (x) [t(x)]-E q θ (x) [f † (t(x))] (9)

where f is a convex function with derivative f and Fenchel conjugate f † ; this divergence class instantiates many familiar divergences, such as the KL and Jensen-Shannon divergence. The variational formulation introduces the functions t(x) whose optimum is related to the density ratio since t * (x) = f (r(x)). Substituting t * in (9), we transform the objective (10) into supremum over r φ (which is attained when r φ = r * = p * /q θ ). For self-consistency, we flip the sign to make it a minimisation problem in r φ , leading to the bi-level optimisation:

EQUATION

Ratio loss:

EQUATION

Generative loss:

EQUATION

where we derived equation ( 12) by extracting all the terms involving q θ (x) in equation ( 10), used the change of variables to express it in terms of the underlying generative model and flipping the sign to obtain a minimisation. There is no discriminator in this formulation, and this role is taken by the ratio function. We minimise the ratio loss, since we wish to minimise the negative of the variational lower bound; we minimise the generative loss since we wish to drive the ratio to one. By using the function f (u) = u log u, we recover an objective using the KL divergence; when f (u) = u log u -(u + 1) log(u + 1), we recover the objective function in equation ( 6) from the previous section using the Bernoulli loss, and hence the objective for GANs.

The density ratio implies that p * (x) ≈ p = r φ (x)q θ (x), since it is the amount by which we must correct our model q θ (x) to match the true distribution. This led us to a divergence formulation that evaluates the divergence between the distributions p * and p, using the KL divergence:

EQUATION

where the first equation is the KL for un-normalised densities (Minka, 2005) . This leads to a convenient and valid ratio loss since all terms independent of r can be ignored. But we are unable to derive a useful generative loss from this expression ( 14) , since the third term with log q cannot be ignored, and is unavailable. Since the generative loss and ratio losses need not be coupled, any other generative loss can be used, e.g., equation (8). But this is not ideal, since we would prefer to derive valid ratio losses from the same principle to make reasoning about correctness and optimality easier. We include this discussion to point out that while the formulation of equation ( 9) is generally applicable, the formulation (14), while useful for ratio estimation, is not useful for generative learning. The difficulty of specifying stable and correct generative losses is a common theme of work in GANs; we will see a similar difficulty in the next section.

The equivalence between divergence minimisation and class probability estimation is also widely discussed, and most notably developed by Reid and Williamson (2011) using the tools of weighted integral measures, and more recently by Menon and Ong (2016) . This divergence minimisation viewpoint (9) was used in f -GANs and explored in depth by Nowozin et al. (2016) who provide a detailed description, and explore many of the objectives that become available and practical guidance.

2.4. RATIO MATCHING:

A third approach is to directly minimise the error between the true density ratio and an estimate of it. If we denote the true density ratio as r * (x) = p * (x) /q θ (x) and its approximation as r φ (x), we can define a loss using the squared error:

EQUATION

where the final objective is obtained by ignoring terms independent of r φ (x) and is used to derive ratio and generative losses. When used to learn the ratio function, Sugiyama et al. (2012b) refer to this objective as KL importance estimation procedure (KLIEP). Concurrently with this work, Uehara et al. (2016) recognised the centrality of the density ratio, the approach for learning by ratio matching, and its connections to GANs (Goodfellow et al., 2014) and f -GANs (Nowozin et al., 2016) , and provide useful guidance on practical use of ratio matching.

We can generalise ( 16) to loss functions beyond the squared error using the Bregman divergence for density ratio estimation (Sugiyama et al., 2012a; Uehara et al., 2016) , and is the unifying tool exploited in previous work (Reid and Williamson, 2011; Sriperumbudur et al., 2009; Sugiyama et al., 2012a; Menon and Ong, 2016) . This leads to a minimisation of the Bregman divergence B f between ratios:

EQUATION

where we have used p * = r * q θ , and D f is the f -divergence defined in equation ( 9). We can derive a ratio loss from ( 19) by extracting all the terms in r φ , leading to the minimisation of L B (r φ ). The role of the discriminator in GANs is again taken by the ratio r that provides information about the relationship between the two distributions. This ratio loss, as pointed out by Uehara et al. (2016) , is equivalent to the ratio loss we derived using divergence minimisation (10), since:

EQUATION

)

The equivalence of the second terms in ( 20) and ( 21) can be derived by using the definition of the dual function:

f † (f (x)) = max r rf (x) -f (r). The maximum is attained at x = r, leading to the identity f † (f (r φ (x))) = r φ (x)f (r φ (x)) -f (r φ (x)).

If we follow the strategy we used in previous sections to obtain a generative loss, by collecting the terms in equation ( 19) dependent on q θ , we obtain:

EQUATION

)

But this does not lead to a useful generative loss since there are terms involving q θ whose density is always unknown, similar to the difficulty we encountered with divergence minimisation in equation ( 14). Since we can use any generative loss that drives the ratio to one, we can employ other losses, like the alternatives described in section 2.2 . One approximation suggested by Uehara et al. (2016) is to assume p * ≈ r φ q θ , i.e. assume a near-optimal ratio, which reduces the f -divergence to:

D f [p * (x) q θ (x)] = E q θ (x) f p * q θ (x) ≈ E q θ (x) f q θ (x)r φ (x) q θ (x) = E q θ (x) [f (r φ (x))] (23)

As a result, the last two terms in equation ( 22) cancel, leaving a generative loss that can be used for learning. Using ratio matching under the Bregman divergence we obtain the bi-level optimisation:

Ratio loss:

EQUATION

Generative loss:

EQUATION

2.5. MOMENT MATCHING:

A final approach for testing is to evaluate whether the moments of the distributions p * and q are the same, i.e. by moment matching. We compare the moments of the two distributions by minimising their distance, using test statistics s(x) that provides the moments of interest:

EQUATION

The choice of test statistics s(x) is critical, since ideally, we wish to match all moments of the two distributions. When the functions s(x) are defined within a reproducing kernel Hilbert space, we obtain kernel-based forms of these objectives, which are highly flexible and allow easy handling of data such as images, strings and graphs. The objective (26) can then be re-expressed in closed-form in terms of kernel functions, leading to the maximum mean discrepancy criterion (MMD). The role of the MMD for learning in implicit generative models defined by neural networks has been explored by both Li et al. (2015) and Dziugaite et al. (2015) . While typically expensive to compute, there are now highly efficient approaches for computing the MMD (Chwialkowski et al., 2015) . The objective using feature matching by Salimans et al. ( 2016) is a form of moment matching that defines the test statistics using intermediate layers of the discriminator function. And a further set of objective functions is possible by framing this problem as one of density difference (rather than ratio) estimation (Sugiyama et al., 2013) .

The other significant body of research that focusses on learning in implicit generative models using the lens of moment matching is approximate Bayesian computation (ABC) (Marin et al., 2012) . The models in ABC are often related to problems in population genetics and ecology, where knowledge of these systems leads to ABC being easily exploited to learn simulator parameters. The ABC literature has widely-explored the use of test statistics, including fixed functions and kernels, and Markov chain Monte Carlo methods to learn the posterior distribution of the parameters. There is a great deal of opportunity for exchange between GANs, ABC and ratio estimation in aspects of scalability, applications, and theoretical understanding. A further insight obtained from a moment-matching approach is that other empirical measures such as the Dudley and Wasserstein distances can easily be considered, establishing even further connections to other growing areas, such as optimal transport (Frogner et al., 2015; Sriperumbudur et al., 2009) .

The moment matching approach can be generalised by representing them as an integral probability metric (Sriperumbudur et al., 2009) , which is what makes it possible to describe the relationships and ways of transforming from density ratio estimation using using moment-matching, f -divergences, and class-probability estimation. Sriperumbudur et al. (2009) showed that f -divergences and integral probability metrics (MMD, Wasserstein) intersect only at the total variation distance. As we described previously, by cleverly formulating our problem as a Bregman divergence minimisation, we show that all the methods we described for density ratio estimation are very closely related (Sugiyama et al., 2012a; Reid and Williamson, 2011; Sriperumbudur et al., 2009; Sugiyama et al., 2012b; Uehara et al., 2016) .

3. DISCUSSION:

By using an inferential principle driven by hypothesis testing, we have been able to develop a number of indirect methods for learning the parameters of generative models. These methods do not compute the probability of the data or posterior distributions over latent variables, but instead only involve relative statements of probability by comparing populations of data from the generative model to observed data. This view allows us to better understand how algorithms such as generative adversarial networks, approximate Bayesian computation, noise-contrastive estimation, and density ratio estimation are related. Ultimately, these techniques make it possible for us to make contributions to applications in climate and weather, economics, population genetics, and epidemiology, all areas whose principal tools are implicit generative models.

Distinction between implicit and prescribed models. The distinction between implicit and prescribed models is useful to keep in mind for at least two reasons: the choice of model has direct implications on the types of learning and inferential principles that can be called upon; and it makes explicit that there are many different ways in which to specify a model that captures our beliefs about data generating processes. Any implicit model can be easily turned into a prescribed model by adding a simple likelihood function (noise model) on the generated outputs, so the distinction is not essential. And models with likelihood functions also regularly face the problem of intractable marginal likelihoods. But the specification of a likelihood function provides knowledge of p * that leads to different algorithms by exploiting this knowledge, e.g., NCE resulting from class-probability based testing in un-normalised models (Gutmann and Hyvärinen, 2012) , or variational lower bounds for directed graphical models. We strive to maintain a clear distinction between the choice of model, choice of inference, and the resulting algorithm, since it is through such a structured view that we can best recognise the connections between research areas that rely on the same sets of tools.

Model misspecification and non-maximum likelihood methods. Once we have made the choice of an implicit generative model, we cannot use likelihood-based techniques, which then makes testing and estimation-by-comparison appealing. What is striking, is that this leads us to principles for parameter learning that do not require inference of any underlying latent variables, side-stepping one of the major challenges in statistical practice. This piques our interest in more general approaches for non-maximum likelihood and likelihood-free estimation methods, of which there is much work (Lyu, 2011; Gutmann and Hyvärinen, 2012; Hall, 2005; Marin et al., 2012; Frogner et al., 2015) . We often deal with misspecified models where q θ cannot represent p * , and non maximum likelihood methods could be a more robust choice depending on the task (see figure 1 in (Huszár, 2015) for an illustrative example).

Choice of loss function. This density ratio viewpoint has led us to derive multiple different objective functions and an understanding of how they are related. But it has left us unsatisfied since we have not gained the insight needed to choose between them. Sugiyama et al. (2012a) make statements on this choice in the simpler setting where only the density ratio is to be estimated. But when we wish to learn implicit generative models, this choice, at present, remains unclear.

Perceptual losses. Several authors have also proposed using pre-trained discriminative networks to define the test functions since the difference in activations (of say a pre-trained VGG classifier) can better capture perceptual similarity than the reconstruction error in pixel space. This provides a strong motivation for further research into joint models of images and labels. However, it is not completely unsupervised as the pre-trained discriminative network contains information about labels and invariances. This makes evaluation difficult since we lack good metrics and can only fairly compare to other joint models that use both label and image information.

Bayesian inference. We have mainly discussed point estimation methods for parameter learning.

It is also desirable to perform Bayesian inference in implicit models, where we learn the posterior distribution over the model parameters p(θ|x), allowing knowledge of parameter uncertainty to be used in risk minimisation and other decision-making tasks. This is the aim of approximate Bayesian computation (ABC) (Marin et al., 2012) . The most common approach for thinking about ABC is through moment-matching, but as we explored, there are other approaches available. An approach through class-probability estimation is appealing and leads to classifier ABC (Gutmann et al., 2014) .

We have highly diverse approaches for Bayesian reasoning in prescribed models, and it is desirable to develop a similar breadth of choice for implicit models.

Evaluation. Our view suggests that value of the density ratio, a quantity we can always compute using the approaches described, can be a useful metric to report. But density ratio estimation has not yet illuminated a method for consistent evaluation of the generative model that is learnt. Without such a measure, it is hard to make fair comparisons; but this is a problem faced by all related fields, and establishing such measures is an important contribution to be made. Furthermore, we have also not yet gained the tools to make theoretical statements that allow us to assess the correctness of the model learning framework, although theoretical developments in the literature on approximate Bayesian computation may help in this regard, e.g., Frazier et al. (2016) .

Non-differentiable models. We have restricted our development to implicit models that are differentiable. In many practical applications, the implicit model (or simulator) will be non-differentiable, discrete or defined in other ways, such as through a stochastic differential equation . The stochastic optimisation problem we are generally faced with (for differentiable and non-differentiable models), is to compute ∆ = ∇ θ E q θ (x) [f (x)], the gradient of the expectation of a function. As our exposition followed, when the implicit model is differentiable, the pathwise derivative estimator can be used, i.e. ∆ = E q(z) [∇ θ f (G θ (z))] by rewriting the expectation in terms of the known and easy to sample distribution q(z). It is commonly assumed that when we encounter non-differentiable functions that the score function estimator (or likelihood ratio or reinforce estimator) can be used; the score-function estimator is ∆ = E q θ (x) [f (x)∇ θ log q θ (x)]. For implicit models, we do not have knowledge of the density q(x) whose log derivative we require, making this estimator inapplicable. This leads to the first of three tools available for non-differentiable models: the use of the weak derivative and related stochastic finite difference estimators, which require forward-simulation only and compute gradients by perturbation of the parameters (Glasserman, 2003; Fu, 2005) .

The two other approaches are: moment matching and ABC-MCMC (Marjoram et al., 2003) , which has been successful for many problems with moderate dimension; and the natural choice of gradientfree optimisation methods, which include familiar tools such as Bayesian optimisation (Gutmann and Corander, 2016) , evolutionary search like CMA-ES, and the Nelder-Mead method, amongst others (Conn et al., 2009) . For all three approaches, new insights will be needed to help scale to high-dimensional data with complex dependency structures.

Ultimately, these concerns serve to highlight the many opportunities that remain for advancing our understanding of inference and parameter learning in implicit generative models.

:



Figures and Tables:

Proper scoring rules that can be minimised in class probability-based learning of implicit generative models.

Table content: <table><tr><td>Loss</td><td>Objective Function (</td></tr></table>

Peer review (iclr2017):
Thank you for an interesting read. Given the huge interest in generative modelling nowadays, this paper is very timely and does provide very clear connections between methods that don*t use maximum likelihood for training. It made a very useful observation that the generative and the discriminative loss do **not** need to be coupled with each other. I think this paper in summary provides some very useful insights to the practitioners on how to select the objective function to train the implicit generative model. The only reason that I decided to hold back my strong acceptance recommendation is that I don*t understand the acceptance criteria of ICLR. First this paper has the style very similar to the Sugiyama et al. papers that are cited (e.g. presenting in different perspectives that were all covered in those papers but in a different context), making me unsure about how to evaluate the novelty. Second this paper has no experiment nor mathematical theorem, and I*m not exactly sure what kinds of contributions the ICLR committee is looking for.