extfieldsofstudy: Economics Mathematics

s2fieldsofstudy: Economics Mathematics

text:
Robust estimation of superhedging prices

We consider statistical estimation of superhedging prices using historical stock returns in a frictionless market with d traded assets. We introduce a plugin estimator based on empirical measures and show it is consistent but lacks suitable robustness. To address this we propose novel estimators which use a larger set of martingale measures defined through a tradeoff between the radius of Wasserstein balls around the empirical measure and the allowed norm of martingale densities. We establish consistency and robustness of these estimators and argue that they offer a superior performance relative to the plugin estimator. We generalise the results by replacing the superhedging criterion with acceptance relative to a risk measure. We further extend our study, in part, to the case of markets with traded options, to a multiperiod setting and to settings with model uncertainty. We also study convergence rates of estimators and convergence of superhedging strategies.

Introduction
Computation of risk associated to a given financial position is one of the fundamental operations market participants have to perform. For institutional players, like banks, it is regulated by the Basel Committee [56] which dictates rules and requirements for such risk assessments. A golden standard has long been given by Value-at-Risk (VaR), however more recently this is being replaced by convex risk measures like Average VaR (Expected Shortfall) or more sophisticated approaches which include market modelling. Consequently, there is an abundant literature on VaR estimation and some more recent works related to statistical estimation of law-invariant risk measures, see [13,44,45,43,58,19,20]. All these works consider a static situation with no trading involved.
In contrast, in this paper we consider estimation of risk for an agent who can trade in the market to offset her risk exposure. To put in evidence the novelty and relevance of our setting, we concentrate on one, simple but canonical, way to assess risk: the superhedging price. Consider a one-period frictionless market with prices (S t , S t+1 ) denominated in units of a fixed numeraire. The current stock prices S t are known and the future prices S t+1 are modelled as random variables, say with return r := S t+1 /S t drawn from a distribution P on R d + . For a payoff g : R d + → R, its superhedging price is given by: (1.1) π P (g) := inf{x ∈ R | ∃H ∈ R d s.t.
In this simple setting, an arbitrage strategy is H ∈ R d such that P(H(r − 1) ≥ 0) = 1 and P(H(r − 1) > 0) > 0 and if no such strategy exists we say that noarbitrage NA(P) holds. By the Fundamental Theorem of Asset Pricing, absence of arbitrage is equivalent to existence of a probability measure Q, equivalent to P, under which S is a martingale, i.e., E Q [r] = 1. There might be more than one such measure and they can all be used for pricing. Taking the supremum over E Q [g] enables to compute the maximal feasible price for g and this, by the fundamental pricing-hedging duality, is the same as the superhedging price of g: for all Borel g, cf. [22,Thm. 1.31]. Despite its theoretical importance and practical relevance, to the best of our knowledge, there has been no attempt to study statistical estimation of the superhedging price. Our paper fills this important gap. Instead of postulating a measure P, we build estimators of π P (g) directly from historical observations of returns r 1 , . . . , r N and study their properties. Furthermore, we extend the estimators to take into account also the option price data. This is practically relevant and methodologically novel in that it allows a coherent and simultaneous use of historical time-series data with current option price data or, in mathematical finance jargon, the physical measure data and the risk neutral measure data. In contrast, in existing approaches historical returns are, if at all, only used indirectly to compute π P (g). In classical mathematical finance, one first postulates a family of plausible models {P θ : θ ∈ Θ}. Such choice may be influenced by stylised features of historical returns, see [26] for a recent example, as well as by other considerations, e.g., of computational tractability. Thereon, historical returns are not used and only the "future facing" options price data is exploited to select a candidate pricing measure Q θ . More recently, pioneered by Mykland [52,54,53] in a continuous-time setting and pursued within the so-called robust approach to pricing and hedging, it was suggested to use historical returns to select a prediction set, i.e., the set of paths on which the superhedging property is required, and then to compute the resulting cheapest superhedge which trades in stocks and options, see [33,6]. Our approach inherits from that perspective but takes a statistical viewpoint and evolves it into a dynamic and asymptotically consistent methodology.
To describe our approach, suppose we observe d-dimensional historical returns r 1 = S 1 /S 0 , . . . , r N = S N /S N −1 and for simplicity assume that these are nonnegative i.i.d. realisations of a distribution P which satisfies the no-arbitrage condition. We can equivalently represent the observations through their associated empirical measuresP which are well known to converge weakly to P as N → ∞, see [66,Theorem 19.1,p. 266]. This suggests a very natural way to approximate the superhedging price by simply usingP N in place of P. We show in Theorem 2.1 below that the resulting plugin estimatorπ N (g) := πP N (g) is asymptotically consistent: lim N →∞π N (g) = π P (g), P ∞ -a.s., where P ∞ denotes the law of the process (r N ) N ≥1 . However, we also show that π N has serious shortcomings. First, it is not (statistically) robust: small perturbations of P can lead to large changes in the distribution ofπ N . We argue that the Lévy-Prokhorov metric used in the classical definition of statistical robustness, Definition 2.6, is not appropriate when looking at the financial context of derivatives pricing. We propose and study alternative metrics and ensuing notions of statistical robustness in Section 4. Second, the plugin estimator also lacks robustness from the financial point of view of risk management. In fact,π N is monotone in N and converges from below so it is always a lower estimate of the risk:π N ≤ π P . In Theorem 2.11, and in more detail in [55,Section 9], we study the convergence rates for the plugin estimators. This, in the one-dimensional case d = 1, could be exploited to build conservative estimates for the superhedging price π P .
A first intuition to improve the plugin estimator could be to turn to estimators of the support of P. Indeed, the superhedging price π P (g), say for a continuous g, only depends on P via its support. We could thus replace theP N -a.s. inequality in the plugin estimator by an inequality on an estimator of the support of P. Such estimators are well studied in statistics, going back to [27,11,16,28], see also [14,41,50,31,59,65,15,9]. Unfortunately, this approach does not seem to hold any ground. First, convergence of support estimators usually imposes strong conditions on P, e.g., compactness and convexity of the support and/or existence of a density. Second, for the convergence of the superhedging prices we would also need to impose some uniform continuity assumptions on g. However, under such conditions on P and g, we could directly improve the plugin estimator and consider a suitableπ N + a N , see Section 2.4.
Instead, to address the shortcomings of the plugin estimator, we propose novel estimators, which we introduce in Section 3. They exploit the dual formulation of the superhedging price in (1.2). In order to achieve financial robustness and to increase our point estimates we need to consider a larger class of martingale measures. Thus we consider π Q N (g) = sup where B N (P N ) is some "ball" in the space of probability measures around the empirical measureP N . We show that this can lead to a consistent estimator if we use a sufficiently strong metric, e.g., the Wasserstein infinity metric W ∞ . In general however such Q N is too large. Instead, our main insight is to consider a tradeoff between the radius of the balls and the behaviour of martingale densities: Q N := {Q ∈ M | dQ/dP ∞ ≤ k N for someP ∈ B p ε N (P N )}, where B p ε N (P N ) denotes the p-Wasserstein ball of radius ε N aroundP N and ε N → 0 as well as k N → ∞. With a suitable choice of ε N , k N , we establish consistency of πQ N (g) for a regular g, see Theorem 3.6, and also their financial robustness, see Corollary 3.7. This also allows us to study the cases when the estimator naturally extends to the setting of superhedging under model uncertainty about P, see Corollary 3.8. The statistical robustness of πQ N (g) is shown in Section 4, see Theorem 4.2. In Section 5 we extend our analysis to the case when risk is assessed not using the superhedging capital but rather via a generic risk measure ρ admitting a Kusuoka representation ( [46], see (5.1) for a definition). We stress that this is substantially different to all the works recalled at the beginning of this Introduction since we consider an agent who can trade and optimises her position to offset the risk. We propose an estimator, inspired by πQ N (g), and show its consistency. Finally, we also propose another estimator: which is inspired by penalty methods used in risk measures and their representations as non-linear expectations. Asymptotic consistency of this estimator is shown in Thereom 3.12 and holds for an arbitrary measurable bounded g. The rest of the paper is organised as follows. In Section 2 we study the plugin estimatorπ N : its consistency, convergence rates and robustness, both statistical and financial. In Section 3 we propose improved estimators and establish consistency for all of them, under different sets of assumptions. Subsequently, in Section 4, we discuss statistical robustness of all the estimators. We show in particular that no estimator can be robust in the classical sense of Tukey-Huber-Hampel, and suggest ways to amend the classical definition to make it more appropriate to the superhedging price estimation. We then detail further applications of the main results. In particular we discuss convergence of superhedging strategies in Section 6 and partially extend the results to a multiperiod setting in Section 7. Section 8 discusses estimators π Q N for generic sets of martingale measures Q N and derives necessary and sufficient conditions for asymptotic consistency of estimators. In particular, it motivates the estimators studied in Section 3. Finally, Section 9 studies in more detail convergence rates of the plugin estimatorπ N when d = 1 and is auxiliary to Section 2.4. Lastly, Appendix A contains the remaining proofs, along with auxiliary results. Notation. We write P(A) for the set of probability measures on A ⊂ R d . P n ⇒ P denotes weak convergence of measures. P ∈ P(R d + ) is a generic distribution for returns r so that E P [r] = R d + xP(dx). We let M = {Q ∈ P(R d + ) : E Q [r] = 1} denote the set of martingale measures for the stock prices (S t , S t+1 ), where S t > 0 is fixed and S t+1 = rS t . We write M A for the set of martingale measures supported on A. We say that P ∈ P(R d + ) does not admit arbitrage, or that NA(P) holds, if {Q ∈ M : Q ∼ P} = ∅. Above, and throughout, H is a row vector, r is a column vector and 1 denotes either a scalar or a column vector (1, . . . , 1) T .

The plugin estimator
Recall that we want to build an estimator for the superhedging price π P (g). The easiest and possibly most natural way to do this is simply to replace the measure P with the empirical measuresP N . This yields the plugin estimator : In this section we develop the necessary tools to show asymptotic consistency of this estimator and understand its properties. The proofs are reported in Appendix A.
2.1. Consistency. We now state the main result of this section: Theorem 2.1. Let P 1 , P ∈ P(R d + ) and g : R d + → R be Borel measurable. Assume that r 1 , r 2 , . . . are realisations of a time-homogeneous ergodic Markov chain with initial distribution P 1 and unique invariant distribution P such that P 1 P. Then lim N →∞π N (g) = π P (g) P ∞ -a.s., (2.2) where P ∞ denotes the law of the Markov process started from P 1 .
Remark 2.2. The assumptions in the above theorem are standard in econometric theory and cover a variety of models frequently used for modelling of financial returns data. We refer to Corollary A.4 in [55] for sufficient conditions for stationarity (with exponential decay rates) for various random coefficient autoregressive models, e.g., linear and power GARCH and stochastic autoregressive volatility models, which are frequently used for option pricing. Nevertheless we remark that this assumption rules out deterministic trends, structural breaks and seasonalities, which need to be treated separately.
The proof for a general g follows by Lusin's theorem from the case of a continuous claim g which in turn depends on the characterisation of the superhedging price using concave envelopes, which we now recall.
We define the P-a.s. concave envelope aŝ It is well known that in the definition of concave envelopes above we could take infimum over affine functions instead of concave functions. It follows from the definition of the superhedging price in (1.1) that we have . Properties and computational methods for concave envelopes, or more generally for convex hulls of a set of discrete points, have been studied in many applied sciences and there are a number of efficient numerical routines available for their calculation. Naturally computational complexity increases with higher dimensions. Nevertheless there exist algorithms determining approximative convex hulls, whose complexity is independent of the dimension, see for instance [61].
To establish a dual formulation for the plugin estimator, assume now that P 1 P as well as no-P-arbitrage, NA(P), holds and recall this implies the pricing-hedging duality, cf. (1.2). It turns out that since supp(P N ) ⊆ supp(P) this already implies that NA(P N ) holds for N large enough. More generally we have: Proposition 2.4. Let P ∈ P(R d + ) and (P N ) N ∈N be a sequence of probability measures on R d + such that P N ⇒ P and supp(P N ) ⊆ supp(P). Then NA(P) In particular, if NA(P) holds then in the setup of Theorem 2.1 we also have lim N →∞π We close this section considering an extended setup where in addition to the traded assets S, whose historical prices we observe, there also exist options in the market, which can be used for hedging g. If the market enlarged with those options does not allow for an arbitrage, the superhedging price of g in this market is again approximated by the plugin estimator, which now also allows for trading in the options. More precisely, we have the following: Corollary 2.5. Let P ∈ P(R d + ) and g : R d + → R be Borel-measurable. In addition to the assets S, assume that there ared traded options with continuous payoffs f 1 (r) and prices f 0 in the market. Define the evaluation map e(r) = r 1 , . . . , r d , f 1 1 (r)/f 1 0 , . . . , fd 1 (r)/fd 0 T andP := P • e −1 . Finally assume no arbitrage, NA(P), holds. Then, under the assumptions of Theorem 2.1, we have P ∞ -a.s.
It is worth stressing that in the classical approach to pricing and hedging, the historical returns are seen as physical measure inputs and might be used, e.g., for extracting stylised features which models should exhibit. In contrast, option prices f 0 are risk-neutral measure inputs and would be used to calibrate the pricing measures. To the best of our knowledge consistent use of both in one estimator has not been achieved before.

Statistical robustness.
Robustness of estimators is concerned with their sensitivity to perturbation of the sampling measure P. To formalise this, suppose we have a sequence of estimators T N which can be expressed as a fixed functional T : P(R d + ) → R evaluated on the sequence of empirical measures, i.e., T N = T (P N ). This is clearly the case with the plugin estimator of the superhedging price in (2.1). Hampel [30] proposed the following definition of statistical robustness: where d L is the Lévy-Prokhorov metric We sometimes say that T N is robust with respect to d L to stress the dependency on the particular choice of the metric. A classical result of Hampel, see [36,Thm. 2.21], states that if T is asymptotically consistent, i.e.
T N = T (P N ) −→ T (P), for all P ∈ P(R d + ) then T N is robust at P if and only if T (·) is continuous at P. The following theorem characterises weak continuity of the superhedging price and hence also robustness of its estimators. In particular, it implies that even for i.i.d. returnsπ N is robust only for special combinations of g and P.
Theorem 2.7. Let g be continuous and P ∈ P(R d + ). Then the functionalP → πP(g) is lower semicontinuous at P. It is continuous if and only if In consequence, any asymptotically consistent estimator T N is robust at P only if the above equality holds true.
In particular we see that, in general, the plugin estimatorπ N (g) is not robust w.r.t. d L . The fact that this holds for any asymptotically consistent estimator suggests strongly that the classical definition of robustness is not adequate in the present context. The superhedging price π P (g) is concerned with the support of P in the sense that for P 1 , P 2 ∈ P(R d + ) with equal supports, and for a continuous g, we have π P1 (g) = π P2 (g). In contrast, any δ-perturbation in the Lévy-Prokhorov sense allows for arbitrary changes to the support, see Lemma A.1. In particular, even if P satisfies no-arbitrage, measures in its neighbourhood may not and one may not employ (1.2) for these. To control the support, we can consider d H (supp(P), supp(P)), for P,P ∈ P(R d + ) and where d H denote the Hausdorff metric on closed subsets of R d + . Proposition 2.8. Let g : R d + → R be uniformly continuous and let P ∈ P(R d + ) such that NA(P) holds. Then the functional P(R d + ) P → πP(g) is continuous w.r.t. the pseudo-metric d H (supp(P), supp(P)).
Alas, this does not allow us to recover statistical robustness of the plugin estimator as the pseudometric above does not admit control over the tails of P. Instead, in Section 4.2, we consider a stronger W ∞ metric which allows to obtain an analogue to Hampel's robustness result.
2.3. Financial robustness. The plugin estimatorπ N not only lacks statistical robustness, as seen above, but is also not a financially robust estimate of risk. In fact, if P 1 P, it converges to the superhedging price from below, i.e.,π N π P . From a risk-management perspective one would like to find a consistent estimator for the P-a.s. superhedging price converging from above. However, as we now show, this is not possible in general. As a direct consequence of the discontinuity of the superhedging functional with respect to the Lévy-Prokhorov metric d L , the convergence from above at some confidence level cannot be achieved in practical applications.
Proposition 2.9. Let P ∈ P(R d + ) satisfy NA(P) and g be bounded and Lipschitz continuous. Then, there exists no consistent estimator T N of π P (g) such that for a confidence level α ∈ [0, 1] there exists N 0 ∈ N and for all P ∈ P(R d + ). Thus, in order to achieve the above property (2.7) it is necessary to make additional regularity assumptions on P and g. We show that this is possible for suitably conservative estimators, see Section 3.2 below. In the case of the plugin estimator, we can never achieve convergence from above but we can develop an understanding of the order of magnitude of the difference π P (g) −π N (g). We first do this by studying the convergence rates, see also [55,Section 9]. Secondly, we achieve this via notions of statistical robustness suited for the plugin estimator, see Section 4.2.
2.4. Convergence rates. We now investigate the convergence rate in (2.4). While motivated by financial considerations, the question is of independent interest. We focus on the one-dimensional case. We let F P be the cumulative distribution function of P ∈ P(R + ) and d N = sup r∈R+ |FP N (r) − F P (r)| denote the Kolmogorov-Smirnov distance betweenP N and P. Definition 2.10. For N ∈ N and k = 1, . . . , 1/(3d N ) we define the interquantile distance if P has bounded support, 0 otherwise. and let κ N = sup k∈{0,..., 1/(3d N ) } κ N k .
We can now establish the speed of convergence for the plugin estimator.
Remark 2.12. When the support of P is bounded, the above result holds for all continuous g. Furthermore κ N tends to 0 as N → ∞.
3. Improved estimators for the P-a.s. superhedging price In the last section we have seen that the plugin estimator is asymptotically consistent but has important shortcomings from a statistical and financial point of view. To address these, we propose now new estimators and investigate their asymptotic behaviour as well as their robustness. To construct these, we consider "balls" around the empirical measureP N and we rely on recent convergence rate results ofP N to P for the choice of the radii.
We start by considering balls in the Wasserstein−∞ metric, which offers a very good control over the support but where we need to make strong assumptions on P to control the rate of convergence forP N . Subsequently, in Section 3.2, we consider Wasserstein−p metrics, p ≥ 1. The use of weaker metrics allows us to treat all measures admitting suitable finite moments but requires a penalisation over the dual (pricing) measures. In fact, our estimators rely on a suitable combination of results on convergence of empirical measures with insights into pricing and control over martingale densities. Similarly to the spirit of Corollary 2.5 above, we combine the physical measure-and the risk neutral measure-arguments, see (3.4). We note that using Wasserstein metrics, as opposed to weaker metrics, allows us to control the first moment which is important for no-arbitrage reasons, see [55,Section 8]. Finally, in Section 3.3, we consider much larger balls, indeed all of M, and let penalisation select the appropriate measures. Short proofs are given here, the proofs are reported in Appendix A.2.
3.1. Wasserstein W ∞ balls. When considering robustness of the plugin estimator we saw that to consider measures in a ball aroundP N we have to consider a notion of distance which, unlike the Lévy-Prokhorov metric, controls the supports. This is achieved by the Wasserstein-∞ distance where Π(P,P) denotes the set of all probability measures γ with marginals P and P and where the equality between the definition and the second representation is a consequence of the Skorokhod representation theorem. A direct comparison of (3.1) with (2.5) reveals that W ∞ controls support in the way that d L does not. However, one immediate issue with considering W ∞ is that if P has unbounded support then W ∞ (P,P N ) = ∞ for all N ∈ N sinceP N are finitely supported. For this reason, and also to obtain appropriate confidence intervals, in order to build a good estimator using W ∞ -balls we have to impose relatively strong assumptions on P: Assumption 3.1. The measure P is an element of P(A) for a connected, open and bounded set A ∈ B(R d + ) with a Lipschitz boundary. Furthermore P admits a density ρ : A → (0, ∞) such that there exists α ≥ 1 for which 1/α ≤ ρ(r) ≤ α on A.
Under the above assumption, we have explicit bounds on W ∞ (P,P N ). The case d = 1 follows from Kiefer-Wolfowitz bounds while the case d We let B ∞ ε (P) denote a W ∞ -ball of radius ε around P. The above lemma allows to deduce consistency of the estimator based on such W ∞ balls: πP(g) π P (g), as N → ∞, P ∞ − a.s.
Remark 3.4. The practical use of W ∞ estimators requires a good handle on l N . Its dependence on the set A is mild -from [64] we see that if a bi-Lipschitz homeomorphism φ : A → [0, 1] d exists, then the constant C depends on the domain A only via the Lipschitz constant of φ. However, the knowledge of α requires uniform a priori estimates on the density of P on A. This should be contrasted with W p estimators below, which only require finiteness of certain moments.
Proof of Proposition 3.3. From Lemma 3.2, an application of Borel-Cantelli shows that, P ∞ -a.s., for N large enough P ∈ B ∞ l N (P N ) and, in particular,π ∞ N ≥ π P . Further, as d H (supp(P), supp(P)) ≤ W ∞ (P,P) by (3.1), for allP ∈ B ∞ l N (P N ) we have supp(P) ⊆ supp(P) 2l N , a compact on which g is uniformly continuous. For measures supported on this compact, Proposition 2.8 yields continuity ofP → πP with respect to W ∞ , which in turn implies consistency ofπ ∞ N and concludes the proof.
Thusπ ∞ N is not only consistent but also financially robust. We shall see in Corollary 4.6 below, that it is also statistically robust with respect to W ∞ . However, these results only hold for measures P which satisfy Assumption 3.1. In the next section we introduce a family of estimators which exhibit similar desirable properties for a much larger class of measures P.

3.2.
Wasserstein W p balls and martingale densities. We assume no-arbitrage NA(P) holds and exploit (1.2) to consider estimators of the form for different specifications of the sets of martingale measures Q N based on "balls" aroundP N . In order to guarantee asymptotic consistency we have to ascertain that the true measure P is contained in these balls and that we have some control over the tails of the martingale measures in Q N . Our crucial insight, following recent work of [21], is to work with Wasserstein metrics defined, for p ≥ 1 and P,P ∈ P(R d + ) with a finite p th moment, by where Π(P,P) is the set of probability measures on R d + × R d + with marginals P andP. In case p = 1, [40] showed that Kantorovitch-Rubinstein duality (see [17,Theorem 11.8.2, p.421]) has a particularly nice expression: One's first intuition might be to use Q N = D p ε N ,∞ (P N ) in (3.2). Interestingly, this does not work as the balls are too large. Indeed, Wasserstein distance metrises weak convergence and Lemma A.1 shows that any ball aroundP N includes measures with full support. As it turns out, to obtain a consistent estimator a subtle interplay is required between ε and k in (3.4).

Assumption 3.5.
(1) r 1 , r 2 , . . . are realisations of a time-homogeneous ergodic Markov chain with initial distribution P 1 and unique invariant distribution P such that P 1 P and dP 1 /dP L 2s (P) < ∞ for some s > 3.
We again refer to [55,Corollary A.4] for examples of processes, which satisfy Assumption 3.5. Clearly Assumption 3.5.2 implies 3.5.1. Under this assumption [23], see also [21], used concentration of measure techniques to obtain rates of the decay for P ∞ (W(P,P N ) ≥ ε), see Lemma A.6 and [55,Lemma A.9]. This allows to compute explicitly a function ε N : We say that Assumption 3.5 holds if either Assumption 3.5.1 holds and then ε N is given in (A.8) or Assumption 3.5.2 holds and ε N is then given in (A.1). We state now the main result in this section.
Theorem 3.6. Let g be either Lipschitz continuous and bounded from below or continuous and bounded, p ≥ 1 and P ∈ P(R d + ) satisfying NA(P). Suppose Assumption 3.5 holds and β N ∈ (0, 1) satisfy The above result shows that πQ N is an asymptotically consistent estimator of π P . Note that we assume no arbitrage NA(P) so that, using (1.2), the convergence above is equivalent to We writeQ p N =Q N when we want to stress the dependence on p. As mentioned above, the consistency depends crucially on the choice ofQ N . We discuss this further and motivate the above choice in [55,Section 8]. For p > 1, D p ε,k (P) are weakly compact but D 1 ε,k (P) is not even weakly closed in general, see Lemma A.7. In case of p = 1, taking weak closure ofQ 1 N could destroy the asymptotic consistency of the estimator, e.g., taking g(r) = (r − 1) in the example in the proof of Lemma A.7.
For the particular choice of β N = exp(− √ N ) under Assumption 3.5.2 an explicit computation yields that for N large enough we have However many other choices of β N are feasible. The essential point is that for summable (β N ), a Borel-Cantelli argument implies that for N large enough the true distribution P is within an ε N (β N )-ball aroundP N . This allows us to deduce a sufficient condition for financial robustness of our estimator: Then πQ N (g) ≥ π P (g) for N large enough so that the estimator is asymptotically consistent and converges from above.
The condition (3.7) is motivated by an approximation result, see Lemma A.8. It allows us also to consider the case when we are unsure about the true measure P and instead prefer to superhedge under all measures in its small neighbourhood. Corollary 3.8. In the setup of Theorem 3.6, fix C > 0 and assume there exists Q ∈ M, dQ/dP ∞ < C. Consider C N → C and a fixed ε > 0. Then holds in P ∞ -probability and P ∞ -a.s. whenever ∞ N =1 β N < ∞. We close this section with two examples illustrating that the assumptions on regularity of g in Theorem 3.6 can not be easily relaxed. and Example 3.10 (g bounded, discontinuous). Set g(r) = 1 {r =1} and consider (r N ) N ≥1 i.i.d. from P = δ 1 . Let Let us remark that πQ N (g) acting on infinite dimensional spaces is bounded by a more sophisticated version of the plugin estimator. To see this define the Average Value at risk of g at level 1/k, for k ≥ 1, by In dimension one, d = 1, it can be re-expressed, see [22,Thm. 4.47], as which makes the link with the classical Value-at-Risk apparent. If we now include the ability to trade and optimise the final position, by the translation-invariance of AV @R P 1/k (·) we can write inf which is a superhedging price, where the acceptance cone is given by an AV @R constraint. An analogous representation and bounds for πQ N follow:  (1) (left) and g(r) = (r − 2) + , P = exp(N (0, 1)) (right). Results averaged over 10 3 runs.
Corollary 3.11. In the setup of Theorem 3.6, let g be 1-Lipschitz and P ∈ P(R d on a set of probability greater or equal than 1 − β N .
Note that The proof proceeds by using a min-max argument to interchange the two suprema and the infimum above and uses continuity ofP to W 1 , see [58] and Appendix A.2 for details. We provide a method for the direct calculation of the Wasserstein estimator implemented in TensorFlow 2 , which is based on recent duality results obtained in [18]. As this approximation is computationally quite costly when a large sample size is used, we opt to compute the upper bound of πQ N (g) given in (3.8) instead. This is shown in Figure 1  balls aroundP N . This estimator allowed us to address fundamental shortcomings of the plugin estimator but, as the counterexamples demonstrated, it is only asymptotically consistent under suitable regularity assumptions on g and/or further assumptions on P. To construct an estimator which would be consistent also for discontinuous payoffs while preserving some of the desirable robustness properties of πQ N , it is natural to turn to penalty methods used in risk measures and their representations as non-linear expectations. Namely we use the maximum norm of the Radon-Nikodym derivative, rather than the Wasserstein distance, in the penalisation term.
Theorem 3.12. In the setting of Theorem 2.1, let NA(P) hold and let g : R d + → R + be Borel-measurable and bounded by some constant C > 0. Then for any The direct implementation of (3.9) proves numerically expensive and unstable due to the fraction dQ/dQ ∞ appearing in the penalisation term. Thus, in Figure  2, we show an upper bound on the penalty estimator derived in the proof of Theorem 3.12 in Appendix A.2. We focus on more tractable properties of the plugin and Wasserstein estimator for the rest of the paper.

Statistical robustness of superhedging price estimators
Recall that in Theorem 2.7 we showed that classical robustness in the sense of Hampel can not hold unless the superhedging price is trivial in that (2.6) holds. This is closely related to properties of Lévy-Prokhorov metric, see Lemma A.1, and we are naturally led to consider stronger metrics than d L , which offer a better control on the support. Below, we investigate the use of Wasserstein distances. First we consider W p for p ≥ 1 which is sufficient to establish robustness of the estimator πQ N from Section 3.2. Then we turn to an even stronger metric W ∞ which is needed to study the plugin estimator.

4.1.
Robustness with respect to the Wasserstein-Hausdorff metric. Following [48] we consider Wasserstein-Hausdorff distance, i.e., a Hausdorff distance between subsets of P(R d + ) equipped with W p : Definition 4.1. Let P,P ⊆ P(R d + ). The p-Wasserstein-Hausdorff metric between sets P andP is given by In this generality W p (P,P) can take the value infinity. Properties of this quantity are discussed in [48] assuming compactness and uniform integrability of P and P. We apply this distance to the sets of the formQ is neither compact nor uniformly integrable, see Lemma A.7. We used these sets in Section 3.2 to define consistent estimators πQ N , see (3.2) and Theorem 3.6. The following establishes their robustness: The estimator πQ N studied in Theorem 3.6 is robust with respect to the p-Wasserstein-Hausdorff metric in the sense that Proof. Note that for all g ∈ L 1 and Q i ∈Q i N , i = 1, 2, we have Taking the infimum over all these probability measures γ yields for all p ≥ 1. The claim follows. Remark 4.3. It follows in particular that if P 1 , P 2 admit no arbitrage then lim N →∞ W p (Q 1 N ,Q 2 N ) = 0 implies supp(P 1 ) = supp(P 2 ). Indeed, otherwise there exists a Lipschitz continuous function g such that Robustness with respect to W ∞ and perturbations of the support. We reconsider now robustness of the plugin estimator from Section 2. In analogy to the previous section, it seems natural to simply consider the Hausdorff distance between the supports ofP 1 N andP 2 N . In Proposition 2.8 we established continuity of P(R d + ) P → πP(g) in the pseudometric d H (supp(P), supp(P)) but noted that it was not sufficient for a robustness result. Recalling the Lévy metric in (2.5), if d L (P,P) ≤ ε thenP can be obtained from P by redistributing ε mass to arbitrary points on R d + , while (1 − ε) mass can only be moved in an ε-neighbourhood (in the Euclidean distance) of where P allocated mass. As we have observed before, the former operation causes problems, as it changes the null sets of the measure uncontrollably. This is no longer possible under our pseudometric. However, to obtain robustness, we have to restrict redistribution of mass to an ε-neighbourhood for all sets and not only for the whole support. This is achieved by the W ∞ metric as is clear from the second representation in (3.1). This leads to the following extended notion of robustness: The following asserts robustness of the plugin estimator in the above sense and is the main result in this section.
such that NA(P) holds. Then, for a uniformly continuous g, the plugin estimatorπ The proof of Theorem 4.5 is reported in [55,Section A.3]. There are ways to weaken the continuity assumption on g and obtain robustness on some P ⊆ P(R d + ), see [55,Corollary A.11]. We close this section with a result on robustness ofπ ∞ N (g) from Proposition 3.3.
Corollary 4.6. Let g be a continuous function and P ∈ P(R d + ) satisfying NA(P) and Assumption 3.1. Then, the estimatorπ ∞

Risk measurement estimation
The P-a.s. superhedging price π P (g) is a very conservative assessment of risk of a short position in a liability with payoff g. Instead, we could use a risk measure ρ P for such an assessment, as proposed by [10], leading to Note that we include above the ability to trade in the market in order to (optimally) reduce the risk of g. We consider ρ P with Kusuoka's representation for a set P of probability measures on [0, 1]. This is not very restrictive since this representation, first obtained in [46], holds for any law invariant coherent risk measure, see [37]. Importantly, it enables us to think of ρ P (g) as a function of the underlying measure P. Much like we did for π P (g), we would like to estimate π ρ P (g) directly from the observed stock returns. To this end we introduce the following estimator as the natural equivalent to πQ N (g). In particular, if P = {δ α } where α ∈ [0, 1], we simply have ρ P (g) = AV@R P α (g) and the corresponding estimator π ρ B p ε N (β N ) (P N ) (g) resembles the Wasserstein W p estimator for fixed level 1/k N := α. We have the following consistency result: Proposition 5.1. Assume g satisfies |g(r) − g(r)| ≤ L γ |r −r| γ for some γ ≤ 1 and L γ ∈ R and that sup µ∈P 1 0 µ(dα)/α γ/p < ∞. Then for any P satisfying NA(P) and Assumption 3.5 the limit in P ∞ -probability holds. If Assumption 3.5.2 is satisfied, then the limit also holds P ∞ -a.s.
Finally, we present a simple empirical test of the performance of our estimators. We simulate weekly returns according to a GARCH(1, 1) model: where ω = 0.02, β = 0.1, α = 0.8 and η n is standard student-t distributed with µ = 5 degrees of freedom. We take 1000 samples from the above P ∼ GARCH(1, 1) and calculate the plugin estimator π AV@R0.95 P N ((r − 1) + ) and the Wasserstein es- We compare this to a parametric estimator of π AV@R P 0.95 ((r − 1) + ), where we first estimate the parameters of the GARCH(1,1) model and then compute π AV@RP 0.95 ((r − 1) + ) given the estimated modelP. For each of these estimates, we use a running window of 50 weeks, which is in line with the Basel III regulations for calculating the 10-day AV@R (see [57,MAR33,p.89]), which set the minimum length of the historical observation period to be one year. We also consider the case when the parameters of the model change for observations 330 − 670. The behaviour of the three estimators is shown in Figure 3. Both the Wasserstein and plugin estimator approximate the true value reasonably well -the Wasserstein estimator being the most conservative estimator. The parametric estimator exhibits the most erratic behaviour which is due to the unstable estimation of the GARCH(1,1) parameters with only 50 data points. This shows advantages of our proposed estimators when compared with a parametric approach, even when the model is correctly selected. The last pane in Figure 3 uses S&P500 weekly returns data from 01/01/2006-01/01/2015 with a moving window of 50 weeks. The GARCH(1,1) estimator, for which the model is mis-specified, does not pick up any markets movements, while both the Plugin and Wasserstein estimator detect the financial crisis and its aftermath. For similar plots but with GARCH(1,1) using log-returns, we refer to [55,Section A.4]. While preliminary, we believe that this simple empirical study points to clear advantages of our approach. In particular, it is encouraging to see that in the last pane, the Wasserstein estimators clearly picks up the financial crisis and the Eurozone debt crisis periods. A further in-depth study of comparative performance of different estimators is clearly needed and is left for future research.  The data is from P ∼ GARCH(1, 1) (first pane) and its variant with a change in the parameters for the middle third of the time series. The last pane uses S&P500 returns.

Convergence of Superhedging strategies
Given the consistency results of Sections 2 & 3 establishing convergence of superhedging prices we now turn to the convergence of the corresponding superhedging strategies. We start with the case of the plugin estimator.
Theorem 6.1. Consider P ∈ P(R d + ) satisfying NA(P). Suppose N ≥1 supp(P N ) = supp(P). Then for any sequence of trading strategies (H N ) N ∈N satisfyingπ there exists a subsequence (H N k ) k∈N converging to some H ∈ R d which satisfies Proof. Note that we can assume without loss of generality that H N ∈ lin(supp(P)− 1) for all N ∈ N. NA(P) implies that the sequence (H N ) N ∈N is bounded. Indeed, assume towards a contradiction that H N → ∞. Clearly H N /|H N | →H ∈ R d with |H| = 1. Also sinceπ we concludeH(r − 1) ≥ 0 on supp(P). By NA(P)H(r − 1) = 0 follows P-a.s. and asH ∈ lin(supp(P) − 1)) we haveH = 0, which contradicts |H| = 1. This shows that (H N ) N ∈N is bounded. Thus there exists a subsequence (H N k ) k∈N of (H N ) N ∈N such that H N k → H. Lastly we note that for each r ∈ {r 1 , . . . , } we have Thus the claim follows for continuous g. As in the proof of Theorem 2.1 we can then extend the result to general g using Lusin's theorem. This concludes the proof. In general one cannot replace the claim H N k → H in Theorem 6.1 by H N → H ∈ R d as the following example shows. Example 6.3. Take d = 1 and Then the sequence of superhedging strategies (H N ) N ∈N is uniquely determined as the slope of the line through the points ( which diverges. Next we establish a corresponding result for the Wasserstein estimator πQ N (g). Recall the definition of AV@R given in Section 3.
Theorem 6.4. Let g be Lipschitz continuous and bounded from below. Also let N (β N )) and lim N →∞ πQ N (g) = π P (g). Then for every sequence of trading strategies (H N ) N ∈N satisfying there exists a subsequence (H N k ) k∈N converging to some H ∈ R d which satisfies π P (g) + H(r − 1) ≥ g(r) P-a.s.
Proof. As in Corollary 3.11, we see that Recall lim N →∞ πQ N (g) = π P (g). As g is Lipschitz continuous Lemma A.5 shows that for someC > 0. Note that similarly to the arguments in the proof of Theorem 6.1 we can assume that (H N ) N ∈N is bounded, so a subsequence of (H N ) N ∈N (which we also denote by (H N ) N ∈N for convenience) converges to some H ∈ R d . Fix N ∈ N. Then for all M ≥ N we have Let ε > 0. Then by (6.2) and N ≤ M large enough since we have As ε > 0 was arbitrary, this implies This concludes the proof.
We quote the following result: We can now formulate the following multiperiod analogue of Theorem 2.1: Theorem 7.4. Let P 1 , P be probability measures on R d + and g : (R d + ) T → R be Borel-measurable. Assume that the observations r 1 , r 2 , . . . are i.i.d. samples of P. Then P ∞ -a.s.
..,r N } 0,T = g P 0,T = π P T (g). 3 The more general Markovian case is likely to require a genuinely novel approach, including using a modified definition of the empirical measure, see [1]. We leave it for future research.
Proof. We prove the claim by induction over T ∈ N. The case T = 1 follows from Theorem 2.1. Thus we assume that we have shown that for each By Lusin's theorem there exists a sequence of compact sets K n ⊆ supp(P) such that P(R d + \ K n ) ≤ 1/n and g P 1,2 | Kn is continuous. As in the proof of Theorem 2.1 we have g P 0,2 = (g P 1,2 (·)) P (1) = lim n→∞ (g P 1,2 (·)) Kn (1). As the concave envelope of pointwise increasing functions is increasing, we conclude that g {r1,...,r N } 1,2 (1) is increasing in N ∈ N. Thus P ∞ -a.s.
Interchanging limits (as they are suprema) yields This concludes the proof for T = 2. The general induction step follows analogously.
As g is bounded from below, passing to the limit with N → ∞ gives the result. Now we show the"≤"-inequality. This follows directly from the fact that r T → g(r 1:T −1 , r T ) is 1-Lipschitz, so by the proof of Theorem 3.6 For the induction step we note that for some set C ⊆ P(R d + ), t ∈ 1, . . . T and some 1-Lipschitz function f : In particular the functions sup q T ∈M, q T ∼P R d Remark 7.6. Similar statements are valid for the penalty estimator and bounded functions g as well as for the W ∞ -estimator with continuous g.

Asymptotic Consistency, Arbitrage and Contiguity
We discuss now the consistency of π Q N (g) in (3.2) for general sets of martingale measures Q N . This, in particular, provides a detailed motivation for the construction of our improved estimator in Theorem 3.6. Throughout, we assume NA(P) to be able to use the dual formulation (1.2) and recall that NA(P), by the First Fundamental Theorem of Asset pricing, is equivalent to existence of Q ∈ M, Q ∼ P.
Clearly, a first necessary condition for asymptotic consistency of π Q N (g) is that Q N = ∅ for N large enough. For the plugin estimator, when Q N = {Q ∈ M : Q ∼P N }, this was established in Proposition 2.4 which outlined the relationship between NA(P) and NA(P N ). The proof crucially relied on the fact that supp(P N ) ⊆ supp(P). Indeed, for general measures ν N ⇒ P the affine hull aff(supp(ν N )) could be of higher dimension than aff(supp(P)) and thus there is no relationship between NA(P) and NA(ν N ) as the following example shows: It is thus both natural and necessary to maintain a relationship between Q N and {Q ∈ M : Q ∼P N }. A minimal property seems to be one of an asymptotic inclusion in the sense that for some sequence (k N ) N ∈N with lim N →∞ k N = ∞ we have Then the consistency of the plugin estimator in Theorem 2.1 implies that The main task now is to identify sequences of sets Q N which satisfy (8.1) and for which the reverse inequality holds. For this, we need the Q N to asymptotically decrease to the set {Q ∈ M : Q ∼ P}. More formally, denoting the ε-neighbourhood of a set A by A ε , the following can be seen to be a necessary condition for (8.2) when g is continuous and bounded: Though (8.3) is a necessary condition for asymptotic consistency of π Q N (g) it is not sufficient as the following example shows: Example 8.3. Let P = 1/3(δ 0 + δ 1 + δ 2 ), then sup Q∼P, Q∈M E Q ((1 − r) + ) = 1/2. Note that for ν N = 1/3(1 − 1/N )(δ 0 + δ 1 + δ 2 ) + δ N /N we have ν N ⇒ P and Q N = (1 − 3/N )δ 0 + 1/N (δ 1 + δ 2 + δ N ) ∼ ν N . Furthermore Q N is o-contiguous with respect to ν N . Nevertheless Q N ⇒ δ 0 and thus E Q N ((1 − r) + ) > 3/4 for N > 12, so asymptotic consistency is not satisfied.
The above example shows that even though (Q N ) N ∈N is o-contiguous wrt. (ν N ) N ∈N , Q N converges weakly to a measure Q, which is not a martingale measure. A way to resolve this problem, is to demand uniform integrability of (Q N ) N ∈N as was done in [34] in a continuous-time setting. This allows to establish the following: Lemma 8.4. Assume NA(P) holds. Let P N be such that for every sequence uniformly integrable, then π Q N (g) is asymptotically consistent for all bounded and continuous g.
We omit the proof since relying on uniform integrability is not possible in general: if supp(P) is unbounded (2) Obviously δ 0 is not a martingale measure, so (Q N ) N ∈N cannot be uniformly integrable.
Thus uniform integrability of (Q N ) N ∈N is in general too strict a requirement. In order to resolve this problem we strengthen the assumption, that ν N converge weakly to the true measure P. Instead we look at convergence in Wasserstein distance, which is known to metrize the weak convergence (cf. [67, Theorem 6.9, p. 96]).

Convergence rates for the plugin estimator
We study now in more detail the convergence rates for the plugin estimatorπ N in the one-dimensional case. In particular, we prove Theorem 2.11 and also state some extensions.
Following the approach of Hampel, it would be natural to consider the influence curve of the superhedging functional, which is simply its Gateaux derivative at P in the direction of δ r and represents the marginal influence of an additional observation with value r when the sample size goes to infinity. This however produces trivial results: the function ε → π (1−ε)P+εδr (g) is constant on (0, 1) so the influence curve at (r, P) will be equal to zero or infinity. This happens because (1 − ε)P + εδ r have the same support for all ε ∈ (0, 1). Clearly, to assess sensitivity of π P to changes in P we have to vary the support. To this end we consider as ε → 0 and study its natural normalisation.
The following examples shows that ε is not the correct normalisation.
Example 9.1. Take again g(r) = |r − 1| ∧ 1 and note that instead of considering P = λ [0,2] /2 we can thin out the tails of P by setting dP n dP [0,1] = r n n + 1 and dP n dP [1,2] = (2 − r) n n + 1 Naturally, π P = π Pn but with increasing n the probability mass is less well spread over the support of P n . We calculate for n ≥ 2 that F P n (r) = r n+1 /2 for r ≤ 1 and thus F −1 P n (p) = n+1 √ 2p for p < 1/2. This readily implies The above examples motivates using quantile functions for normalisation as stated in Theorem 2.11 which we now prove.
Proof of Theorem 2.11. As we have noted beforeπ N (g) ≤ π P (g) andπ N (g) is nondecreasing in N . Let us first consider P with bounded support. It suffices to show that sup Q∼P, Q∈M Without loss of generality we assume r 1 ≤ r 2 ≤ · · · ≤ r N for the rest of the proof. By the definition of d N holds for all p ∈ [0, 1]. Thus we note that for all i = 2, . . . , N Next we remark that by definition of π P (g) and Proposition 2.8 there exists C > 0 such that We want to apply a simple Balayage construction to redistribute the mass of Q on the support ofP N . Thus we set A straightforward calculation shows Q N ({r 1 , . . . , r N }) = 1 and E Q N [r] = 1. We have For P ∈ P(R + ) with unbounded support and g bounded by D > 0 we note that .
This concludes the proof.
In order to improve upon, and further specify the results in Theorem 2.11, we recall Lemma A.8 and make the following easy observation: Lemma 9.2. Let C > 0 and P ∈ P(R d + ). Then the set {Q ∈ M | dQ/dP ≤ C} is weakly compact.
the claim follows. Corollary 9.3. Let P ∈ P(R + ) have bounded support and let g : supp(P) → R be continuous such that |g(r) − g(r)| ≤ δ(|r −r|) for some monotone δ : R + → R + with δ(r) → 0 for r → 0. If (3.7) holds then sup Q∼P, Q∈M Proof. Note that by assumption it is sufficient to consider martingale measures Q ∼ P such that dQ/dP ∞ ≤ C for some C > 0. Similarly to the proof of Theorem 2.11 we set N , ∞)). and note that Q N ({r 1 , . . . , r N }) = 1 and We further assume that g is bounded by D. Then (9.1) becomes This concludes the proof.
Remark 9.4. We note that the above asymptotic result can be used to set up a utility based hedging problem to approximate π P . If we let ) for a concave and strictly increasing U and some C ∈ R + then is the value of the utility based hedging problem underP N inf{x ∈ R | ∃H ∈ R s.t. U (x + H(r 1 ) − g(r)) ≥ α N for r = r 1 , . . . , r N }.
For examples of stationary β-mixing Markov chains with exponential decay we refer to Lemma A.3 and Corollary A.4. The above lemma can be used to obtain asymptotic confidence bounds for Corollary 9.3. More precisely, for N → ∞, we obtain for some constant C > 0 and f : R + → R + being the square of the inverse of x → x + δ(x). Similarly, under some regularity assumptions on F −1 P , an analogous but non-asymptotic estimate can be given for Theorem 2.11 using Lemma 2.13.
To close this section we consider the convergence rate in some cases when Theorem 2.11 does not apply. Note that Theorem 2.11 applies for a continuous g whenever P has bounded support or if there exists K > 0 such that Suppose now that no such K exists. Clearly if g(r)/r → ∞ as r → ∞ then π P (g) = ∞ so consider g with linear growth: g(r)/r → c ∈ R as r → ∞. As P necessarily has unbounded support we can take a sequence K n → ∞ and, by the above condition, some (H n ) n∈N such that π P(·|[0,Kn]) (g) + H n (K n − 1) = g(K n ).
Example 9.6. We provide now some examples to illustrate different convergence rates which might be observed in the context of Theorem 2.11 or (9.3) above. We use N = 10 5 realisations of 1000 runs for our numerical illustrations. We start with some examples when Theorem 2.11 applies either directly or because (9.2) holds.
We move now to examples where we can not rely on Theorem 2.11 but use (9.3) instead and show the bound may be sharp. The asymptotic distribution of max i=1,...,N r i can be determined using classical results from extreme value theory. In particular, the scaled maximum of exponential/normal/lognormal random variables converges weakly to a Gumbel distributed random variable Y (see [63,Examples 2,3,p.199]).

Appendix A. Additional results and proofs
A.1. Additional results and proofs for Section 2.
Proof of Theorem 2.1. First note that if A n is a non-decreasing sequence of sets with A = lim n A n = ∪ n A n thenĝ A = lim nĝAn . The "≥" inequality is obvious and the reverse follows sinceĝ An is a non-decreasing sequence of concave functions thus its limit is a concave function dominating g on A.
Proof of Proposition 2.4. Assume first that NA(P) holds. Denote the relative interior of the convex hull of a set A ⊆ R d by ri(A). Furthermore write lin(A) for the linear hull of A and aff(A) for the affine hull. We recall that 1 ∈ ri(supp(P)) if and only if NA(P) holds. Consequently there exist ε > 0, 0 ≤ k ≤ d and r ± i ∈ conv(supp(P)) with r ± i = 1 ± εe i for i = 1, . . . , k, where e 1 , . . . e k are an orthonormal basis of the space lin(supp(P) − 1). Fix some r ± i . Then there exists n ∈ N andr 1 , . . . ,r n ∈ supp(P) such that r ± i can be written as a convex combination ofr 1 , . . . ,r n . Denote by A the finite collection ofr for all r ± i , i = 1, . . . , k. Choosing 0 < δ < ε sufficiently small we have 1 ∈ ri({r ± i | i = 1, . . . , k}) for each choice ofr ± i ∈ B δ (r ± i ) ∩ aff(supp(P)). As supp(P N ) ⊆ supp(P) and P N ⇒ P there exists N 0 ∈ N such that P N (B δ (r) ∩ aff(supp(P))) ≥ P N (B δ (r) ∩ aff(supp(P N ))) = P N (B δ (r)) > 0 for allr ∈ A and N ≥ N 0 . Thus 1 ∈ ri(supp(P N )) and NA(P N ) holds. Conversely if there exists H ∈ R d such that P(H(r − 1) > 0) > 0 and P(H(r − 1) ≥ 0) = 1 then again by continuity H(r − 1) ≥ 0 on supp(P) ⊇ supp(P N ) and the set {r | H(r − 1) > 0} is open. Thus there exists N 0 ∈ N such that Proof of Corollary 2.5. By NA(P) we have As e is continuousP N ⇒ P implies (P) N ⇒P. Next assume supp(P N ) ⊆ supp(P).
Proof of Theorem 2.7. Fix P ∈ P(R d + ) and a sequence P N converging to P. Let {r 1 , r 2 , . . . } be dense in supp(P). Fix n ≥ 1 and note that, for any i ≥ 1, weak convergence implies that P N (B 1/n (r i )) > 0 for all N large enough. In particular there exists r n i ∈ B 1/n (r i ) such thatĝ P N (r n i ) ≥ g(r n i ). Thus, by the same reasoning as in the proof of Theorem 2.1 above, lim inf We conclude using Lemma A.1 since for a sequence with supp(P N ) = R d + , by continuity of g, we have, for all N ≥ 1, For the second part of the theorem, assume P ∈ P(R d + ) is such that Take a sequence (P N ) N ∈N , as above, with supp(P N ) = R d + and P N ⇒ P. Fix ε > 0 such that 2ε < π P N (g) − π P (g).
For every δ > 0 there exists N 0 ∈ N such that for all N ≥ N 0 we have d L (P N , P) ≤ δ. Let T N be an asymptotically consistent estimator of π P (g). Then, for all N large enough Thus T N is not robust at P, which shows the claim.
Proof of Proposition 2.9. We give the following counterexample: Define ω 1 = (1, 1, ...), P n = (1 − 1/n)δ 1 + (1/n)λ [0,2] and g(r) = |1 − r| ∧ 1. Obviously P n ⇒ δ 1 and sup Q∼δ1 E Q [g] = 0 as well as sup Q∼Pn, Q∈M E Q [g] = 1 for all n ∈ N. By consistency we must have T N (ω 1 ) → 0 as N → ∞, in particular we can assume T N0 (ω 1 ) < 1. Thus A.1.1. Proof of Remark 2.2. With regards to Remark 2.2 let us now define the following concepts: Definition A.2. Suppose that {X n | n ∈ N} is a time-homogeneous Markov chain with initial distribution P as its invariant measure and transition kernel K. Then {X n } is called stationary β-mixing with exponential decay if there exist 0 < ρ < 1 and c > 0 such that where · T V defines the total variation norm. The Markov chain {X n | n ∈ N} is called geometrically ergodic if there exists a probability measure P ∈ P(R d + ), a constant 0 < ρ < 1 and a P-integrable nonnegative measurable function g such that The following lemma lists stationarity and ergodicity properties for some common time series models following [51,5,32,8,3,49,25], where for simplicity we only consider the simplest cases and refer to the references above for a more comprehensive picture: Lemma A.3. The following conditions are sufficient for the β-mixing property with exponential decay (if started from the invariant distribution) for the processes {r n | n ∈ N} and {h n | n ∈ N} defined below, as well as for geometric ergodicity of the process {h n | n ∈ N}: • Augmented GARCH(1,1) models, where r n = h n η n , n = 0, 1, . . . n ] = 1 and continuous positive density with respect to the Lebesgue measure, Γ is an increasing function and e n+1 is a measurable function of η n .
As a direct consequence we obtain the following: Corollary A.4. Let the distribution of {r n } n∈N be given by one of the models satisfying the corresponding conditions in Lemma A.3 for stationarity and let P 1 = P 4 . Then the conditions of Theorem 2.1 are satisfied. Furthermore Assumption 3.5.1 is satisfied in the following cases: • LGARCH: There exists an integer q > 3p such that (β + α) q < 1 or E[|η n | 2q ] < ∞ and β + α < 1/(E[|η n | 2q ]) 1/q . 4 We note that GARCH is traditionally used for logarithmic returns log(Sn/S n−1 ) and not raw returns rn = Sn/S n−1 in the econometric literature. As r → log(r) is continuous this is no restriction for Theorem 2.1, however imposing Assumption 3.5 for rn = log(Sn/S n−1 ) would mean that we had to demand exponential moments of the GARCH models. Even in a parametric setting, statistical inference in a heavy-tailed GARCH environment is notoriously difficult, see e.g. [29]. On the other hand, for short time scales we find rn ∼ 1 and the approximation log(rn) ∼ rn − 1 is quite accurate.
Proof. The first claim follows directly from Lemma A.3. To show the second claim it is sufficient to check that there exists q > 6p with E[|r n | q ] < ∞ and (3.5) holds. Using again Lemma A.3, existence of q > 6p with E[|r n | q ] < ∞ is satisfied under the conditions listed in the statement of the corollary. Lastly we note that and thus conclude that for 0 < ρ < 1 in the β-mixing with exponential decay property of (r n ), as given by Lemma A.3, r 1 ∼ P and all f ∞ ≤ 1 we have which is summable.
A test for (strict) stationarity of GARCH models was developed in [24]. Applying their test to daily returns of CAC, DAX, DJA, DJI, DJT, DJU, FTSE, Nasdaq, Nikkei, SMI, and SP500, from January 2, 1990 to January 22, 2009 the authors conclude that non-stationarity of the time series is not plausible.
A.2. Additional results and proofs for Section 3.
Recall also that A is connected, open and bounded and P satisfies Assumption 3.1. It follows that F P (x) + ε ≤ F P (x + αε) for x ∈ R. We thus conclude that with P ∞ -probability greater than (1 − exp(−2 √ N )) we have, for x ∈ R, and in particular W ∞ (P,P N ) ≤ αN −1/4 . Proof of Corollary 3.8. The "≤"-inequality follows as in the proof of Theorem 3.6. For the "≥ "-inequality take Q N ∈ M and ν N ∈ B p ε+ε N (P N ) such that dQ N /dν N ∞ < C N . By assumption there exists a martingale measure Q ∼ P such that dQ/dP ∞ ≤ C − δ for some δ > 0. Define Then λ N ∈ (0, 1), lim N →∞ λ N = 0 and (1 − λ N )ν N + λ N P ∈ B p ε (P) on a set of probability at least 1 − β N . Furthermore (1 − λ N )C N ≤ C − δ N and thus and, as Q N , Q ∈ M and for anyQ ∈ M there exists M > 0 such that where the last inequality holds with P ∞ -probability at least 1 − β N . This concludes the proof.
We recall here the continuity property of AV @R P 1/k N (g): Lemma A.5 ([58], Cor. 11, p.538). For g ∈ L 1 and P,P ∈ P(R d + ) we have This is used for the subsequent proof.
Proof of Corollary 3.11. Recall that We now want to interchange the two suprema and the infimum above, which can be done by the same arguments as in [2, Proofs of Theorem 3.1 and Lemma 3.2]. Indeed, as NA(P) holds and as P ∈ B p ε N (β N ) (P N ) with high probability we conclude by [22,Theorem 1.48,p.29] that there exists N 0 ∈ N such that where ri(A) denotes the relative interior of the convex hull of a set A ∈ B(R d + ). This implies πQ The first inequality in (3.8) is trivial, while the second inequality follows from 2-Lipschitz-continuity of g(r) − H(r − 1) for |H| ≤ 1 and continuity ofP → AV @RP 1/k N (g − H(r − 1)) w.r.t. to W 1 , see [58] or Lemma A.5.
Proof of Theorem 3.12. Clearly It remains to establish the reverse inequality. Let us first assume C N = C. Fix Q ∈ M and considerQ ∼P N ,Q ∈ M such thatQ Q. Define ε ≥ 0 through Note that we can write Q = f (r)Q + ν for some Radon-Nykodym derivative f and some measure ν singular toQ, where ν(R d In particular a contradiction. Furthermore, by the definition of ε, so there is no gain from shifting mass and in particular sup which is non-negative for ε ≤ 1 − C N /C. Finally, the case whenQ is not absolutely continuous with respect to Q is trivial since then dQ/dQ ∞ = ∞.
A.2.1. Proof of Theorem 3.6. For simplicity of exposition we first prove Theorem 3.6 under Assumption 3.5.2. We adopt the notation of Section 3.
Lemma A.6 ([23] Theorem 2). Under Assumption 3.5.2 we have for N ≥ 1, d = 2p and ε > 0, where c 1 , c 2 are positive constants that only depend on p, d, a, c and E P [exp(c|r| a )]. Thus for some confidence level β ∈ (0, 1) we can choose Proof. Take a sequence Q n ∈ D p ε N (β N ),k N (P) such that Q n ⇒ Q. For every n ∈ N there exists ν n ∈ B p ε N (β N ) (P) such that dQ n /dν n ∞ ≤ k N . First observe that for any n we have It follows, by weak convergence and monotone convergence theorem, that it is enough to consider the case d = 1. Then Consider now the terms on the RHS. The third term can be made arbitrarily small by taking large K since Q admits first moment. The fourth term can be bounded, in analogy to (A.2), as follows: where we took K large enough. Finally, for a fixed K, the difference between the first two terms can be made small by taking n large due to weak convergence of measures. The bound If p > 1, the ball B p ε N (P) is weakly compact (see [67, Def. 6.8., p.96]), in particular uniformly integrable, and it follows that also the fourth term on the RHS of (A.3) converges to zero, uniformly in n, as K → ∞ Further, possibly on a subsequence, (ν n ) n∈N converges weakly to a limit ν ∈ B p ε N (P). By regularity of probability measures it is sufficient to test dQ/dν against bounded continuous functions and we easily conclude that dQ/dν ∞ ≤ k N . This shows compactness of D p ε N (β N ),k N (P).
Consider now p = 1. We give the following counterexample: take P = δ 1 and set for r n ≥ 2 Furthermore let Then W 1 (ν n , δ 1 ) ≤ ε N (β N ) and Assume lim n→∞ r n = ∞. Then Taking the closure of D 1 ε N (β N ),k N (P) would ensure compactness, but consistency of the estimator sup Q∈Q N E Q [g] in Theorem 3.6 would be lost in general since the closure might include non-martingale measures. To see this, take for instance g(r) = (r − 1) in the example in the proof of Lemma A.7.
Proof of Theorem 3.6 under Assumption 3.5.2. Let us assume that Assumption 3.5.2 is satisfied. Note that the "≥"-inequality follows from Lemma A.8. Indeed, Lemma A.8 implies that for all N ∈ N there exists a martingale measure Q N ∼ P with dQ N /dP ∞ ≤ k N such that sup Q∼P, Q∈M with a N → 0 as N → ∞. With P ∞ -probability (1 − β N ) we have P ∈ B p ε N (β N ) (P N ) and hence Q N ∈Q N . This gives We recall that NA(P) gives π P (g) = sup Q∼P, Q∈M E Q [g] and hence, for any ε > 0 and N large enough, P ∞ (πQ N (g) − π P (g) ≤ −ε) ≤ β N → 0 as N → ∞. For the "≤"-inequality, we assume g is Lipschitz continuous bounded from below and we take H ∈ R d such that Then and the second term on the RHS vanishes since Q N is a martingale measure. To treat the last term on the RHS consider the functioñ g := (g − H(r − 1) − π P (g)) ∨ 0 which is non-negative, C-Lipschitz for some C > 0 and {g > 0} = A c . Since P(A c ) = 0 we have in particular P). We conclude that, for any ε > 0, This establishes the convergence of πQ N (g) to π P (g) in P ∞ -probability. Further, whenever ∞ N =1 β N < ∞, a simple application of Borel-Cantelli lemma, similarly as in [21,Lemma 3.7], shows that the convergence holds P ∞ -a.s. This concludes the proof in the case of Lipschitz continuous g bounded from below and under Assumption 3.5.2.
It remains to argue the "≤"-inequality when g is bounded and continuous. We fix a small δ > 0 and definẽ Similarly to Corollary 3.11 we see that for K > 1 large enough Thus there exists a sequence H N ∈ R d such that for all N ∈ N. Take K large enough so that d/K < δ. For notational simplicity we assume that P has full support. Recall that g is bounded, |g| ≤ C. We now show that H i N is bounded. Let us first look at the lower bound: for this, we take i ∈ {1, . . . , d} and suppose that H i N ≤ 0 (otherwise we trivially bound H i N from below by 0). We work on the set {P ∈ B p ε N (β N ) (P N )}, which has P ∞ -probability at least 1 − β N . Now we take N large enough such that k N · P (s j r j ≥ 1 for all j = i, K/2 < r i < K) ≥ 1 for all (s 1 , . . . , s i−1 , s i+1 , . . . , s d ) ∈ {−1, 1} d−1 . Then defining s j = −sign(H j N ) for all j = i and dQ s1,...,si−1,si+1,s d dP := 1 {sj rj ≥1 for all j =i, K/2<ri<K} P(s j r j ≥ 1 for all j = i, K/2 < r i < K) we have by (A.5) and thus H i N ≥ −(2(C + 1))/(K/2 − 1). On the other hand, assuming H i N ≥ 0 and setting dQ s1,...,si−1,si+1,s d dP := 1 {sj rj ≥1 for all j =i, 0<ri<1/2} P(s j r j ≥ 1 for all j = i, 0 < r i < 1/2) we got by the same arguments as above ] + πQ N (g) ≤ 2C + 1 for a possibly larger N . In conclusion thus −(2(C + 1))/( Consider now Q N ∈Q N . Note that by Markov's inequality we have Then where we used (A.6) and (A.7) as well as |H j N |Q N (([0, K] d ) c ) ≤ d(4C + 2)/K for j = i in the second inequality. In the last inequality we used (A.5). By definition of πQ N (g), and since E Q N [H N (r − 1)] = 0, this shows We obtain the same result in the case of P without full support with a possible addition of the term constant times ε N k N on the RHS. Note that we either have r i = 1 P-a.s., or else by NA(P ), P puts mass on r i on either side of 1. If r i is bounded under P, we can takeH i N = H i N and the arguments remain the same but we have the additional error term as Q N can put small mass on unbounded r 1 .
Exchanging the roles of r and r finally yields |ĝ(r) −ĝ(r )| ≤ (L ∨ 2C)|r − r | for r, r ∈ R d + . Thus we can argue as in the proof of Theorem 3.6 for Lipschitz continuous claims (and using the same notation) to obtain for someĈ > 0. In particular As ε was arbitrary and ε N k N → 0 the claim follows.
We have used the following lemma: Before we finish the proof of Theorem 3.6 we first give the convergence rates forP N under Assumption 3.5.1. 5 This is a slight modification of the result in [23] but the proof is essentially the same as that of Theorem 15 therein and is hence omitted. Proof of Theorem 4.5. Let ε > 0 and fix P 0 ∈ P(R d + ) such that NA(P 0 ) holds. As g is uniformly continuous, its P-concave envelope operator is continuous at P 0 , see Proposition 2.8, and we can take δ small enough so that for allP,P ∈ P(R d + ) with d H (supp(P), supp(P)) ≤ 2δ we have |πP(g) − πP(g)| ≤ ε/9.
This concludes the proof.
Corollary A.11. Let P ∈ P(R d + ) such that NA(P) holds. (i) Let g be continuous and P ⊆ P(R d + ). If P ∈ P and for all δ > 0, there exists a compact set K ⊆ R d + such that sup P∈P πP(g) − πP (·|K) (g) ≤ δ, thenπ N (g) is robust at P wrt. W ∞ on P. (ii) Let g be a continuous function of linear growth and P ⊆ P(R + ). If P ∈ P, P is uniformly integrable and for all δ > 0 there exists C > 0 such that sup P∈P πP(g) − sup dQ/dP ∞ ≤C, Q∈M E Q [g] ≤ δ, (A.13) thenπ N is robust at P wrt. W ∞ on P.
The rest of the proof follows as in the proof of Theorem 4.5 above using uniform continuity of g on K. Otherwise there exist sequences (L n ) n∈N , (H n ) n∈N and (P n ) n∈N such that H n → ∞. Note that by uniform integrability ofP n we haveP n ⇒P and alsoP n (·|[0, L n ] d ) ⇒ P, where NA(P) holds. TakeH n := H n /|H n |, then after possibly taking a subse-quenceH n →H with |H| = 1 andH(r − 1) ≥ 0P-a.s., which leads toH = 0 by NA(P), a contradiction. Take  Thus again we can restrict to K = [0, L] d as before and proceed as in the proof of Theorem 4.5.
Proof of Corollary 4.6. Recall that P is compactly supported, say supp(P) ⊆ B R (0) for some R > 0, so that we can assume that g is uniformly continuous. Note that we also have that supp(P N ) ⊆ B R (0), P ∞ -a.s. Theorem 4.5 then shows robustness ofπ N . Thus, for any ε > 0, there exists δ > 0 and N 0 ∈ N such that for all N ≥ N 0 and allP ∈ P(R d + ) we have: W ∞ (P, P) ≤ δ ⇒ d L (LP(π N ), L P (π N )) ≤ ε/3.
A.4. Additional results for Section 5. In this Section we present some simulations complementing Figure 3 in the main article. More specifically we compare estimates for the quantity π AV@R P 0.95 (g), where we set g(r) = (r − 1) + (Figure 10) for historical gold price (WGC/GOLD DAILY USD) returns and Apple (AAPL) returns, and g(r) = |r − 1| (Figure 11) for historical S&P500 and DAX30 returns. Finally, Figure 12 reproduces the case of S&P500 and DAX30 returns and g(r) = (r − 1) + from Figure 3 but this time the GARCH(1,1)-estimator uses logreturns instead of simple returns. While before it ignored market crisis now it overreacts to it. We use 50 or 100 data points to build the estimates and plot the average of the last 5 or 10 running estimates. Figure 10. Comparison of estimates for π AV@R P 0.95 ((r−1) + ). The first two panes show estimates with a rolling window of 100 data points and we plot the average of the last 10 estimates. The first pane uses Apple returns while the second one uses gold returns. The last pane uses gold returns for estimates with a rolling window of 50 data points where we plot the average of the last 5 estimates. Figure 11. Comparison of estimates for π AV@R P 0.95 (|r − 1|). Estimates use a rolling window of 50 data points and we plot the average of the last 5 estimates. The first pane uses DAX30 returns while the second one uses S&P500 returns. Figure 12. Comparison of estimates for π AV@R P 0.95 ((r − 1) + ) estimated on log-returns for the GARCH(1,1) model. Estimates use a rolling window of 50 data points and we plot the average of the last 5 estimates. The first pane uses DAX30 returns while the second one uses S&P500 returns.


Implementation of Tourism Facilities Expansion to Hedonism and Community Imaging

This study examines "Implementation of Expansion of Tourism Facility to Hedonism and Community Imagery in Kutuh Village Ubud Gianyar. The rapid development of tourism in this region led to land conversion is not unstoppable. The landowners in the village are encouraged to rent to investors. Also, the internal dimension consists of economic capital, is the most important element because economic capital can drive all business activities, whether small, medium, and large capital business. Another aspect is imagery, illustrated through perceptions, cognitions, motivations, and attitudes that take place in individual consumers. Like building luxury homes and showing off luxury cars as one of the imaging of the new rich. On the other hand, the external dimension of rapid tourism progress has implications for hedonists and community imaging Kutuh Ubud Gianyar. The consumer society revolves around symbols and signs. Mass media become miracles in the liturgy of objects, enjoyment is understood as the embodiment of freedom, and the human body is the main consumer object. Another aspect of globalization, the discourse of globalization as a process characterized by the rapid development of science, information technology, and telecommunications will be able to change the world fundamentally. Besides Maslow's hierarchy of needs is an interpretation of the level of human needs represented in the form of a pyramid with a more basic need at the bottom. These needs have a level or hierarchy, ranging from the lowest to the highest (self-actualization). Seeing the problem of life of the community I need Ubud Gianyar as something that needs to be faced and solved, not avoided. The last aspect is the government policy, in the regulation is defined as one of the areas of tourism area in the village of Kutuh Ubud Gianyar. So it is not denied the opportunity of lively people to rent the land for the expansion of tourism and other economic driving business.

Introduction
Globalization cannot be separated with the progress of science, information technology, and telecommunications holistically. The influence of global culture flows at Kutuh Ubud Gianyar has implications for cultural practices of capitalism, such as the emergence of a cultural industry that refers to the touristification of cultural forms, as the entertainment industry, mass culture, popular culture, and consumer culture. Cultural tours are in line with the services offered by the tourism industry by selling standard, standard, deluxe and suite room services, as well as performing arts, attractions, architecture, cuisine, as tour packages. The common feature of the accommodation facility is to utilize terraces, cliffs and scenic river as the location of the hotel. This accommodation facility complements the image of Kutuh Village Ubud Gianyar as a typical tourist village regarding structure, nature, and modern tastes. As a blend of high touch and high tech simultaneously offers a distinctive charm and image as an international village.
The working community has a specific purpose that is to meet the needs. Needs cannot be separated from everyday life, during human life requires a variety of needs, such as food, clothing, housing, education, and health. Needs are influenced by culture, environment, time, and religion. The higher the cultural level of a society, the higher the need to be fulfilled. Physiological needs are the most basic and greatest potential for all the fulfillment of the above needs. Hungry people will always be motivated to eat, not to make friends or be rewarded. The Kutuh Ubud Gianyar community will ignore or suppress all other needs until their physiological needs are satisfied. If the Kutuh Ubud Gianyar community is already established, the need to satisfy hunger as a lifestyle. They usually have enough food, but when they say hunger then what they think is the taste of food to be chosen, not the hunger he feels. Someone who is hungry will not be too concerned with the taste, smell, temperature or texture of food.
Hedonism according to Aristotle in Russell (2004: 243) pleasure is different from happiness because there can be no happiness without pleasure. Said three views of pleasure: (1) all pleasures are not good; (2) some good pleasures, but mostly bad; (3) good enjoyment, but not the best. Aristotle rejected the first opinion because suffering is certainly bad, so pleasure is certainly good. Rightly he says that it is absurd to say that humans can be happy in suffering: good luck whose nature is outward, to some extent, necessary for the realization of happiness. He also denies the view that all pleasures are physical; everything contains a spiritual element, and pleasure contains so many possibilities for the enjoyment of ever-present pleasurable pleasures. Furthermore, it is said to be a bad pleasure, but it is not the pleasure felt by good people, it may be different enjoyment. Lifestyle affairs are not just a monopoly of moneyed people either; not even the poor can still use certain hedonistic models even if they are just playing, imitating, or pretending, just as the rich can also pretend to be poor. Not because adhere saving or ascetic ideology is not because of the influence of simple lifestyle campaigns, but because of lifestyle choices. For the community of Kutuh Ubud Gianyar modern matter of imaging and lifestyle is not clear any more limits. In the discourse of capitalism, everything produced by capitalism will eventually be deconstructed by subsequent new production under the laws of progress and novelty.

Literature Review
Until now, a special study on the implementation of tourism Expansion of Hedonism and Imaging of Kutuh Ubud and Bali communities is still very minimal. However, some of the results of studies that have been conducted that are relevant and can be used as references in this study can be put forward as follows. Based on the study of studies on Ubud more focused on the problem of painting and tourism. Several studies have been conducted in the village of Ubud in the last fifteen years, among others (Picard, 2006) with a book entitled Tourism Culture and Tourism Culture (Ubud Case Study). The results of this study indicate that the area of Ubud attracts tourists who want to see Bali in its artistic and ceremonial image. They want to see the original Balinese village. Every day Ubud is flooded with hundreds of special tourists coming from Nusa Dua, Sanur, and Kuta to see the village painter. Also, many visitors to learn the arts that have become a source of fame Ubud Village, the art of dance, sculpture, and painting. So the focus of this study describes the development of tourism and the increasing number of local artists interested in pursuing her ancestral art heritage to remain sustainable.
Tjokorda Arta Ardana Sukawati (2010) conducted research entitled "Social Change of Ubud Traditional Village, Gianyar, Bali in the Era of Globalization of Cultural Studies." In that research, the socio-cultural and socio-economic situation of the readiness of the people of Ubud in anticipation of the impact caused by the lack of distance between countries in the world due to the increasingly widespread globalization swept the world. So, the focus of this study is ready to be unprepared by the Ubud community already in the era of globalization that is rapidly sweeping across all aspects of community life including the social changes caused by the current globalization.

Research Method
The data used in this research are: a. Primary data, i.e., data, and information obtained directly from sources/respondents. This data was obtained by interviewing the Kutuh Ubud Gianyar community home which rented land for tourism expansion, as a research sampling unit. The population of this study is all of Kutuh Ubud Gianyar community who rent the land by using proportional sampling technique with a sample of 27 people. b. Secondary data, namely data and supporting information obtained from the office of Kutuh Village Ubud Gianyar and other relevant data.

Research Sites
Research location in Kutuh Village Ubud Gianyar by observing and mapping the community who rent the land for the expansion of tourism. The observation time as an object of analysis is the condition of 1 (one) year running in 2017. 2017 usage considerations are taken to facilitate the recording of information to the respondents, i.e., the Kutuh Ubud community who rent the land for the expansion of tourism.

Research instrument
The instruments used in this study are qualitative descriptive research through literature surveys and field surveys and in-depth interviews. The interviews were open-ended questions and closed-ended questions.

Data analysis method
The method of analysis used in this study includes several methods by the objectives of the study are as follows.
1) The theory of Hedonism used to analyze society into hedonists, divide consumption behavior into three periods with the social conditions of each surrounding. These three periods, namely the classical period, the emergence of consumption sociology, and the postmodernist period, connect consumption with social symbols in society. 2) Hierarchy Theory The need to describe the real condition that needs at a low level must be fulfilled or at least sufficiently met first before higher-level needs become motivating. Maslow also said that the need for love includes the love that gives and the love that it receives. Society must understand love, must be able to teach it, create it and predict it. Otherwise, the world will drift into a wave of animosity and hatred. a. Observation Techniques and Interviews Observation is conducted by observing directly in the field so that it can see and observe more carefully about the condition of the location that is rented as the expansion of tourism. In-depth interviews (In-depth Interview) with some informants who have set 27 people.

b. Documentation Study
Documentation study is data collection techniques by taking some documents or mass notes ago, both on the people who rent the land in Kutuh Ubud Gianyar and relevant agencies.

Results and Analysis 3.1 Epicurus and Russell Hedonism
Theory Hedonism is a worldview that assumes that people will be happy by seeking as much happiness as possible and avoiding painful feelings as much as possible. Hedonism is the doctrine that pleasure or pleasure is the goal of human life and action. Hedonism has become the grip of the community Kutuh Ubud Gianyar many noble values of humanity fade, even disappear. Their social sensitivity threatened to be displaced when they always consider the profit and loss in socializing. The Kutuh Ubud Gianyar community looks like a living mummy that is of no use to them. They seem to be a royal guard of pleasure that no one can sniff let alone taste it. Other people can only gawk at their establishment. They do not care. As a result, when there are people who need a helping hand, they hide and are reluctant to make sacrifices. The cause of hedonism due to external factors namely the rapid flow of industrialization and globalization that attacks the community Kutuh Ubud Gianyar is an inevitable factor. Values that were once considered taboo are now considered normal. Media communication, especially the advertising media is very tangent to ethical and moral issues. Through imaginative symbols, mass communication media obviously take into account and take advantage of passions, feelings, and desires. Currently, hedonists promote a variety of offers of human needs until the life of the world sparkling night porn action through the action of television media, advertising and other print media (Budi Santoso, 2016: 69).
Hedonistic behavior cannot be separated from the association of fellow in big cities who prefer pleasure and enjoyment. In associating there is always pressure from within the community of Kutuh Ubud Gianyar to do the same with the group. If someone lives in an environment that pleasure-loving life, the pursuit of pleasure, then naturally the person will follow a lifestyle that has been invested in the social environment. Theo Huijbers sometimes say because of being pressured problem demanding economic needs; the metropolitan community can drift increasingly hedonistic consumerism. While the internal factors, the weakness of one's religious beliefs also affect the behavior of some communities Kutuh Ubud Gianyar which glorify the fun and mere hura-rah. Binzar Situmorang (2016:72) stated that "Spirituality someone become a benchmark in everyday life, especially for those who love the pursuit of pleasure.
Hedonism is a lifestyle that emphasizes and elevates material-physical pleasure during life on earth. Ubud Kiani Gianyar people who have a hedonistic lifestyle always assume that heaps of property, position, rank, position, and the like that is outwardly is happiness and joy. The era of globalization of many people Kutuh Ubud Gianyar follows the hedonic lifestyle because at this time many media that offer a fun life with modern style, consumptive, and genset. Lifestyle ideas that are demanded and pursued as perpetrators of modern life are free life without boundaries, both ethics of ethics, moral, and morals. It is also motivated by the tightness of the world of competition, especially in the field of economy and the principles of the fulfillment of needs and unlimited desires. The hedonic lifestyle of this society is not sensitive to the circumstances, but always do things according to desire, not care about the surrounding environment.
In an uncertain condition, they are living a life spree with the luxury of a sophisticated facility. Ironically, all the goods are not necessarily bought from their own money. Some holders of power use vehicles and facilities from corrupt money, while most of those who rent the land in Kutuh Ubud have a luxurious lifestyle of debt or from something unclear. This attitude is what experts call the hedonic lifestyle. Ubud Kiani Gianyar people spend more money for activities that are not clear or all amusement. Important considerations at least activity is no longer a priority, more importantly, it can bring happiness for a moment. Look how the community is competing with luxury, homes, vehicles, mobile phones, and other luxury facilities. All that often is not based solely on necessity, but the emphasis is more on the value of prestige or pride alone. As a result, for those who cannot afford to be tempted to do everything to follow the trend of today. Hedonis like this is an illusion. The illusion gives false happiness. Happiness in a dream that will eventually disappoint itself. It is said so because without realizing it, the cost of living is too great to be wasted (Emanuel Subangun, 2004: 96). Affirmed that this action as something deadly. How not, luxury goods enjoyed by the community Kutuh Ubud Gianyar is imported goods. That is why all spending spent is not for the progress of the nation itself, but for other nations that are economical, socially, politically colonized. That is the most profound ignorance in the individual.
Hedonist and materialist attitudes have also created unhealthy competition amongst people. The existence of this attitude creates an unhealthy social gap. On the other hand, there is a class full of luxury, while others live within limitations. With concern and simplicity, there needs to be a re-awareness of these commendable characteristics. It is important because over time this society will be destroyed due to the attitude of its people who do not support development in this country. The current orientation of wealth will be greatly influenced by the lifestyles of future generations, even at the leader level in the country. Operationally hedonism is the lifestyle of the community Kutuh Ubud Gianyar which emphasizes and elevates the material-physical pleasure during life on earth. Hedonistic groups always assume that heaps of property, position, rank, position, and the like that is outward is happiness and joy. With a hedonist there is selfishness, regard everything as a matter of personal affairs as a human right. So it will not care what others think about what he does. Always indifferent about his or her lifestyle judgment that is good or not. Other factors that cause the tendency of hedonism is the culture of Kutuh Ubud Gianyar community and the surrounding environment. In a community of Kutuh Ubud Gianyar that has a luxurious living culture, the tendency to luxury will dominate all members of society. In this case, the luxuries of officials and public figures will have a tremendous impact on this lifestyle. The mass media factor in this contemporary and digital era of advertising contained in various media means helped to create a culture of hedonism. These media in many cases advertise products that are not needed. These ads also leave behind the psychological impacts on the hedonists. Many of the harm caused by hedonism. First, the disappearance of wealth, increasing the gap between poor and rich poverty, bankruptcy and debt in the midst of small communities. Ibn Khaldun historian and Muslim sociologist in this regard said: To what extent society is immersed in hedonism, so far they will approach the limit of destruction (Ibn Khaldun, 2016: 37). The process of destruction will occur because hedonism will slowly cause poverty in society and the state.
The extent to which hedonism epidemic, so far poverty will spread in the community. On the other hand, wasting a fortune to buy expensive things that are only meant to be proud of, will gradually drag a country to a foreign party. It is what is happening right now in the world. Many countries in the world depend on the West who at all times market new products for consumption. Though the work, effort and effort to seek treasure, can lead a person and his community to progress but do not let that fall into the valley of hedonism.

Abraham Maslow's Hierarchy of Needs Theory
The concept of the basic needs hierarchy of the community Kutuh Ubud Gianyar, takes precedence over the needs of others. As the need for clothing, food and home, people feel thirsty, then people will tend to try to quench thirst. People can live without food for weeks. But without water, the individual can live for only a few days because the need for water is stronger than the need for food. These needs are often called Abraham Maslow as the basic needs described as a hierarchy or ladder depicting the level of need. There are five levels of needs, namely: (1) Physiological Needs (Physiological Needs) or the most basic needs of every person Kutuh Ubud Gianyar namely the need to maintain his life physically. These needs are like the need for food, drink, shelter, sleep and oxygen (clothing, food, boards). Physiological needs are the most basic and greatest potential for all the fulfillment of the above needs. Hungry man will always be motivated to eat, not to make friends or appreciate. Man will ignore or suppress all other needs until his physiological needs are satisfied. In the already established community of Kutuh Ubud Gianyar, the need to satisfy hunger is a lifestyle. They usually have enough food, but when they say hunger then what they really think is the taste of food to be chosen, not the hunger he feels. Someone who is really hungry will not be too concerned with the taste, smell, temperature or texture of food. Physiological needs differ from other IJSSH Vol. 2 No. 1, April 2018, pages: 33~42 needs in two ways. First, the physiological need is the only requirement that can be fully satisfied or minimally insurmountable. Humans can feel enough in the activity of eating so that at this point, the driving force to eat will be lost. For someone who has just completed a great meal, and then imagining a food again is enough to make her nauseous. Second, that is typical in physiological needs is the essence of repetition. After humans eat, they will eventually become hungry again and will continue to search for food and water again. While the needs of the higher levels do not continue to emerge. For example, a person who is at least partially fulfilled their need to be loved and appreciated will still be confident that they can maintain the fulfillment of those needs without having to look for them again. (2) The Need for Safety (Security / Security Needs). Once the physiological needs are adequately satisfied, what Maslow calls the needs of security. These security needs include physical security, stability, dependence, protection and freedom from threatening powers such as crime, war, terrorism, disease, fear, anxiety, anger, unrest and natural disasters. As well as psychic needs that threaten psychiatric conditions such as not mocked, not humiliated, not stressed, and so forth. The need for security is different from physiological needs because this need cannot be completely fulfilled. Humans can never be protected completely from the threats of meteors, fire, flood or other dangerous behavior. According to Maslow, unsafe people will behave just like unsafe children. They will behave as if always in great danger. An insecure person has a need for excessive order and stability and will try hard to avoid things that are foreign and unexpected. (3) The Need for Sense of Ownership and Social Needs. If the physiological needs and the need for security have been met, then the need for love, affection and belonging has arisen. These needs include encouragement to be needed by others to be regarded as citizens of their social community. This form of fulfillment needs such as friendship, the desire to have a partner and heredity, the need to be close to family and interpersonal needs such as the need to give and receive love. Someone whose love needs have been relatively fulfilled since childhood will not feel panic when rejecting love. He will have great confidence that he will be accepted by people who are important to him. When someone else rejects him, he will not feel crushed. For Maslow, love involves a healthy and loving relationship between two people, including mutual trust. Often love becomes corrupted if one party fears its weaknesses and mistakes. Maslow also said that the need for love includes the love that gives and the love that it receives. The community of Kutuh Ubud Gianyar must understand love, must be able to teach it, create it and predict it. Otherwise, the world will drift into a wave of animosity and hatred. (4) Esteem Needs. Once the needs are loved and possessed sufficient, then man will be free to pursue his ego needs on the desire to excel and have prestige. Maslow found that everyone in Kutuh Ubud Gianyar had two categories of reward needs, namely lower and higher needs. Low needs are the need to respect others, the need for status, fame, glory, recognition, attention, reputation, appreciation, dignity, even domination. High demand is the need for self-esteem including feelings, beliefs, competencies, achievements, mastery, independence and freedom. Once humans can meet the need to be respected, they are ready to enter the gates of self-actualization. (5) The Need for Self-actualization Needs (Maslow in Sobur, 2003: 72). The last level of Abraham Maslow's basic needs is self-actualization, namely the need to prove and show himself to others. At this stage, the community of Ubud Gianyar Kutuh develop as much as possible all the potential he has. The need for self-actualization is a need that does not involve balance, but involves a constant desire to fulfill potential. Maslow describes this need as a desire to become fully self-worth, to be anything according to his ability. Maslow hypothesized that after the community of Kutuh Ubud Gianyar satisfies the needs at the lowest level, society will satisfy the needs at the next level. If at the highest level but basic needs are not satisfied, then the community can return to the level of previous needs. According to Maslow, the satisfaction of various needs is driven by two forces namely deficiency motivation and growth motivation (Maslow in Sobur, 2003: 89). The motivation of deficiency aims to overcome the problem of human tension due to various shortcomings. While the growth motivation is based on the capacity of every human being to grow and develop. Capacity is the nature of every human being.

Conclusion
To understand the various variations of the implementation of the expansion of tourism facilities to hedonism and community imaging in Kutuh Ubud Gianyar after the lease of land, the dimensions of the times, the globalization of tourism, and the powers of change, both internal and external. Another thing is the power, money, and actors who play so that efforts to analyze the hedonic dimension and image of society that leases the land and the power of capitalism internally and externally needs to be done. Accordingly, the conclusions drawn from the results of the study are as follows. Firstly, Hedonism among the Kutuh Ubud Gianyar community has grown tremendously with the development of the times, a mindset that only concerns the pleasure of making this society lulled into a sometimes unrealistic life. Another factor that causes the tendency of hedonism is the culture of Kutuh Ubud Gianyar community and the surrounding environment. In a community of Kutuh Ubud Gianyar that has a luxurious living culture, the tendency to luxury will dominate all members of society. In this case, the luxuries of officials and public figures will have a tremendous impact on this lifestyle. The mass media factor in this contemporary and digital era of advertising contained in social media helped create a culture of hedonism. These media in many cases advertise products that are not needed.
Second, Imaging practices subsequently become a community lifestyle. Imagery becomes a new worldview of society. Along with the continuous operation of cross-country industry and the growth of supermarkets, hypermarkets, and malls. In fact, with its beautiful strategy, consumer goods are adapted to the experiences and philosophical views of the local community (Fordism). Human needs because of its nature is endless then prioritize the basic needs first. After that just meet the needs of a secondary or tertiary. For some unmet needs required hard work in sufficient The fulfillment of needs according to this Maslow is expected to make the life of the community Kutuh Ubud Gianyar more peaceful. The construction of a society of capitalism which views production as the driving force of the economy must now be rebuilt that capitalist society sees to consume as an economic driver. Globalization as a process is characterized by the absence of borders between countries in the world.

Recommendation
First, to fortify themselves from hedonism and imagery that offers only momentary pleasure, it must start from the self as well as the support of others. For the people who rent the land in Kutuh Ubud Gianyar should increase self-control so that the lease of land can be used for more productive purposes. Embed the attitude of life sparingly among families, many useful things that can be done without having a spree.
Secondly, this study is a preliminary study on hedonism and community imaging in Kutuh Ubud Gianyar. This study traces the behavior of hedonism and community imaging that receives land rent for expansion of tourism facilities. In this research has been disclosed behavior of hedonism and community imaging in brief, but has not been studied in depth so that needed further research comprehensively.