L evy Computability of Probability Distribution Functions Takakazu - - PowerPoint PPT Presentation

l evy computability of probability distribution functions
SMART_READER_LITE
LIVE PREVIEW

L evy Computability of Probability Distribution Functions Takakazu - - PowerPoint PPT Presentation

< 13Nancy > 1 L evy Computability of Probability Distribution Functions Takakazu Mori Yoshiki Tsujii Mariko Yasugi Kyoto-Sangyo University ~ } < 13Nancy > 2 0. Introduction


slide-1
SLIDE 1

<13Nancy> 1

L´ evy Computability

  • f Probability Distribution Functions

Takakazu Mori Yoshiki Tsujii Mariko Yasugi

森 隆一 辻井 芳樹 八杉 満利子

Kyoto-Sangyo University

~ }

slide-2
SLIDE 2

<13Nancy> 2

0. Introduction

★ ✧ ✥ ✦ P ★ ✧ ✥ ✦ F ★ ✧ ✥ ✦ CL ★ ✧ ✥ ✦ CF Fine case ✻ µ((a, b]) = F (b) − F (a) ❄ F (x) = µ((−∞, x]) ✘✘✘✘ ✘ ✿ Bochner ✘ ✘ ✘ ✘ ✘ ✾

  • R eitxµ(dx)

❅ ❅ ❅ ❅ ❅ ❅ ■ ❅ ❅ ❅ ❅ ❅ ❅ ❘ Inversion Formula ✘✘ ✘ ✿ ✘ ✘ ✘ ✾ L´ evy ❅ ❅ ❅ ❅ ❅ ❅ ❘ ❅ ❅ ❅ ❅ ❅ ❅ ❅ ❅ ■ vague conv.

set of all probability measures on R

  • unif. conv.

positive definite uniformly continuous ’(0) = 1

  • unif. conv.

increasing Lipschitz f(`1) = 0; f(+1) = 1

pointwase conv. at points of continuity

set of all probability distribution functions on R Discontinuous

Figure 1: The spaces P, F, CL

Notations: µ ↔ F ↔ LF = f (unless oterwise stated)

~ }

slide-3
SLIDE 3

<13Nancy> 3

Introduction 1 Classical Theory 2 L´ evy height and L´ evy metric 3 Computability on P 4 L´ evy computabilities 5 L´ evy Computability and Fine Computability for Probability Distribution Functions 6 Effective L´ evy Convergence

~ }

slide-4
SLIDE 4

<13Nancy> 4

Computabilities on F We had treated Fine computability and effective Fine convergence We seek: computability and effective convergence on F relations to Fine computabilities on F relations to computability and effective convergence on P Translated Dirac measure δ 1

3 is computable, but the corresponding

probability distribution functions is not Fine computable. We have an example of probability distribution function F such that F (0) is not computable, but the corresponding µ is computable (Example 4.2) We propose L´ evy computability and effective L´ evy convergence on F.

~

slide-5
SLIDE 5

<13Nancy> 5

Convergence on P motivation convergence ↔ determining class ↔ computability a class A of functions or of sets is a determining class if µ(η) = ν(η) for ∀η ∈ A implies µ = ν

~ }

slide-6
SLIDE 6

<13Nancy> 6

1. Classical Theory about probability measures on R (Summary) Probability Distribution Functions are characterized by (F-i) monotonically increasing (F-ii) right-continuous (F-iii) F (∞) = 1, F (−∞) = 0 Notations: F: the set of all probability distribution functions Cκ: the set of all continuous functions with compact support Cb: the set of all bounded continuous functions µ(f) =

  • R f(x)µ(dx)

{gℓ}: a dense sequence of Cκ {gℓ}-metric : d{gℓ}(µ, ν) =

  • ℓ=1

2−ℓ(|µ(gℓ) − ν(gℓ)| ∧ 1)

~ }

slide-7
SLIDE 7

<13Nancy> 7

The following converences are equivalent ([1], [6]). (i) µm(f) → µ(f) ∀f ∈ Cκ (vague) (ii) µm(f) → µ(f) ∀f ∈ Cb (weak) (iii) lim infm µm(G) ≥ µ(G) ∀G: open domain theory S2007[23] (iv) d{gℓ}(µm, µ) → 0 {gℓ}: dense in Cκ W1999[27][0,1] [0, 1] (iii′) lim infm µm(I) ≥ µ(I) ∀I: open interval W1999[27] Intervals SS2005[22] R (v) ϕm(t) → ϕ(t) for any t (vi) Fm(x) → F (x) ∀x a point of continuity of F ⇔ µ({x}) = 0, {(−∞, x]} (vii) dL(Fm, F ) → 0 L´ evy (other metrices: cf. [2])

~ }

slide-8
SLIDE 8

<13Nancy> 8

2. L´ evy height and L´ evy metric L´ evy [11], Ito [6] Notations: GF = {(x, y) | F−(x) ≤ y ≤ F (x)), x ∈ R} (F−(x) = lim

y↑x F (y) )

∀t ∈ R, (x(t), y(t)) = the unique crossing point of X + Y = t with GF Definition 2.1 L´ evy height: LF (t) = y(t) Definition 2.2 L´ evy metric (distance): dL(F, G) = sup

t∈R

  • LF (t) − LG(t)
  • Definition 2.3

L´ evy convergence: dL(Fn, F ) → 0 Remark 2.4 dL(F, G) is equal to the L´ evy(-Prokhorov) metric, that is, dL(F, G) = inf{ǫ > 0 | F (x − ǫ) − ǫ ≤ G(x) ≤ F (x + ǫ) + ǫ, for all x}

~ }

slide-9
SLIDE 9

<13Nancy> 9 F

a b x = t ` f(t) f(t) = LF (t) ˜ F`(x) t ˜ F (x)

LF

˜ F (a) ˜ F`(b) ˜ F (b) ˜ F`(x) t Dx ˜ F (x)

Figure 2: F and LF

Put LF (t) = f(t). Dx = {t | x = t − f(t)} sup{f(t) | t ∈ Dx} − inf{f(t) | t ∈ Dx} = the jump of F at x

~

slide-10
SLIDE 10

<13Nancy> 10

Properties of L´ evy height Let CL be the set of functions which satisfy the following: (Li) 0 ≤ f(t) − f(s) ≤ t − s if s ≤ t. (Lii) limt→−∞ f(t) = 0. (Liii) limt→∞ f(t) = 1. Proposition 2.5 L: F → CL

  • ne-to-one and onto

CL : closed convex subset of Cb w.r.t. the sup-norm (distance) d∞ Hence, (CL, d∞) is complete

~ }

slide-11
SLIDE 11

<13Nancy> 11

Properties of ˜ F , ˜ F−, LF Dx = [ ˜ F−(x), ˜ F (x)]

  • ˜

F (x) = x + F (x) is strictly increasing and continuous LF (t) = F ( ˜ F −1(t))

  • [ ˜

F−(x), ˜ F (x)) ∩ Range(F ) = φ

  • ˜

F −1(t) is computable if F is computable

  • F (x) = f(sup Dx)
  • ˆ

f(t) = t − f(t) is a nondecreasing continuous function from R onto R

slide-12
SLIDE 12

<13Nancy> 12

Example 2.6 Dirac measures δa, Da: prof. dist. function

  • Da(x) = 0 if x < a and = 1 if x ≥ a
  • LDa(t) =

     if t ≤ a t − a if a ≤ t ≤ 1 + a 1 if t ≥ 1 + a

  • dL(Da, Db) = d∞(LDa, LDb) = |a − b| ∧ 1
slide-13
SLIDE 13

<13Nancy> 13

3. Computability on P We employ computable notions on R by Pour-El and Richards. Definition 3.1 {µm} is computable

⇌ {µm(fn)} is computable for any computable sequence fn

with compact support (L(n) s.t. fn(x) = 0 if |x| ≥ L(n) for some recursive L(n)) Definition 3.2 {µm} converges effectively to µ

⇌ {µm(fn)} converges effectively to {µ(fn)}

for any computable {fn} with recursive compact support Notation: {µm}

e

− → µ

~ }

slide-14
SLIDE 14

<13Nancy> 14

weakly: computable sequence fn with compact support → effectively bounded computable {µ(fn)} ∃M(n): recursive such that |fn(x)| ≤ M(n) Notation: {µm}

ew

− → µ Proposition 3.3 (1) {µm} is computable ⇔ {µm} is weakly computable (2) Assume {µm} and µ are computable. Then, {µm}

e

− → µ ⇔ {µm}

ew

− → µ

~

slide-15
SLIDE 15

<13Nancy> 15

4. L´ evy computabilities Definition 4.1 {Fn} is said to be L´ evy computable

⇌ {LFn} is computable.

Intuitively, L´ evy computability means that we can draw the graph GF effectively. That is, ( ˆ f(t), f(t)) is a one parametric representation of GF, in the sense of Skotokhod [24]. ( ˆ f(t) = t − f(t))

(Section 6, P. 21)

Next Example shows that there exists a probability measure ν: ν is computable G(0) is not a computable real

~ }

slide-16
SLIDE 16

<13Nancy> 16

Example 4.2 (Example 3.4 in [16]) α: one-to-one recursive with non-recursive range. d = ∞

i=1 2−α(i) < 1, dn = n i=1 2−α(i) (d0 = 0)

ν = (1 − d)δ0 +

  • i=1

2−α(i)δ2−(i−1), νn = (1 − dn)δ0 +

n

  • i=1

2−α(i)δ2−(i−1) G and {Gn}: the corresponding probability distribution functions

1 1 ` d 2`¸(1) 2`¸(2) 2`¸(3) 2 1 ` d 1 ` d1 1 ` d2 1 ` d3 2`¸(1) 2`¸(2) 2`¸(3)

Figure 3: Graph of G and LG

~ }

slide-17
SLIDE 17

<13Nancy> 17

The followings holds:

  • {Gn}: monotonically decreasing w.r.t. n, and converges to G
  • {LGn}: monotonically decreasing w.r.t. n, and converges to LG
  • LG(t) = LGn(t) if t ≤ 1 − d or t ≥ 1 − dn+1 + 2−(n+1)

= on (1 − d, 1 − dn+1 + 2−(n+1))

  • {LGn} is computable

We can prove that {LGn} converges to LG effectively uniformly. This implies that LG is computable and G is L´ evy computable. G is not continuous, Fine continuous, not Fine computable, L´ evy computable.

~

slide-18
SLIDE 18

<13Nancy> 18

Theorem 4.3 Assume that {Fm} is continuous. Then {Fm} is L´ evy computable ⇔ {Fm} is computable. Proposition 4.4 µ(g) =

  • R g(x)dF (x) =
  • R g( ˆ

f(t))d f(t) for all g ∈ Cκ Theorem 4.5 If {Fm} is L´ evy computable, then {µm} is computable. Lemma 4.6 Da: the probability distribution function of δa LF is computable, p is positive omputable, a is computable, then (1) L(pDa) is computable. (2) L(pDa + F ) is computable. (3) LF (t) ≤ L(pDa + F )(t) ≤ p + LF (t).

~

slide-19
SLIDE 19

<13Nancy> 19

P → CL µ(g) =

  • R g(x)dF (x) =
  • R g( ˆ

f(t))d f(t) P → F

❇ ❇ ❇ ❇ ❇ ❇ ❇ ❇ ❇ ❇ ❇ ❇

1 c c ` 2`n c + 2`n ˜ w+

c;2−n

˜ w`

c;2−n

µ( ˜ w−

c,n) ≤ F (c) = µ(χ(−∞.c]) ≤ µ( ˜

w+

c,n)

Figure 4: ˜ w+

c,n(x) and ˜

w−

c,n(x)

If µ is computable and c is computable, then µ( ˜ w−

c,n) and µ( ˜

w+

c,n) are computable.

F (c) is right (lower) computable but we cannot derive the computability of F (c) (+ continuity)

slide-20
SLIDE 20

<13Nancy> 20

5. L´ evy Computability and Fine Computability for Probability Distribution Functions Fine topology is generated by {I(k, i) = [ i

2k, i+1 2k ) | k ∈ N, i ∈ Z}

J(x, k) is the unique I(k, i) which contains x. Fine computability is defined with respect to this Fine topology. {ei} is an effective enumaration of dyadic rationals. Definition 5.1 A sequence of functions {fn} is said to be Fine computable if it satisfies (i) (Sequential Fine computability) {fn(xm)} is computable for any Fine computable {xm} (ii) (Effective Fine Continuity) There exists a recursive function α(n, k, i) such that (ii-a) x ∈ J(ei, α(n, k, i)) ⇒ |fn(x) − fn(ei)| < 2−k (ii-b) ∞

i=1 J(ei, α(n, k, i)) = [0, 1) for each n, k

~ }

slide-21
SLIDE 21

<13Nancy> 21

uniformly Fine computable if α(n, k, i) does not depend on i Theorem 5.2 {Fm} is Fine computable ⇒ {Fm} is L´ evy computable First, we prove the following special case. Proposition 5.3 {Fm} is uniformly Fine computable ⇒ {Fm} is L´ evy computable. Outline of the proof of Proposition 5.3 Lemma 5.4 If {µm} is computable, then there exists a recursive function L(m, k) such that µm(wn) > 1 − 2−k, equivalently µ(wc

n) < 2−k, for all n ≥ L(m, k).

(For a single F ) α(k): mopdulus of effective uniform Fine continuity note: F−(x) is computable if x is Fine computable

~

slide-22
SLIDE 22

<13Nancy> 22

Let {tn} be computable ∀k, the set of fite open intervals (−∞, −L(k) + 2−k), (L(k) − 2−k, ∞), ( ˜ F (−L(k) + i2−α(k) − 2−k, ˜ F−(−L(k) + (i + 1)2−α(k) + 2 · 2−k) (0 ≤ i ≤ 2 · L(k)2α(k)) ( ˜ F−(−L(k) + i2−α(k)) − 2−k, ˜ F (−L(k) + i2−α(k)) + 2−k) ( ˜ F (−L(k) + i2−α(k)) − ˜ F−(−L(k) + i2−α(k)) > 2−k) is an open covering of R We can define effectively a sequence {rn,k} such that |f(tn) − rn,k| < 2−k

slide-23
SLIDE 23

<13Nancy> 23

Outline of the proof of Theorem 5.2 Proposition 5.5 ([15]) Let {fm} be a Fine computable sequences of

  • functions. Define

ϕm,n(x) = 2n−1

j=0 fm(j2−n)χI(n,j)(x).

Then, {ϕm,n} Fine converges effectively to {fm}. We can prove that {Lϕn} converges effectively uniformly to LF .

slide-24
SLIDE 24

<13Nancy> 24

6. Effective L´ evy Convergence Definition 6.1 {Fm} is said to L´ evy converge effectively to F

⇌ dL(Fm, F ) converges effectively to zero

{Fm}

eL

− → F Special case of Skorokhod M1-convergence for GADLAC A pair of functions (λ(t), ξ(t)) is said to be a parametric representation

  • f the graph of GF = {(x, z) | F−(x) ≤ z ≤ F (x), x ∈ R}

if GF = {(λ(t), ξ(t) | t ∈ R}, ξ(t) is continuous and λ(t) is continuous and monotonically increasing. (( ˆ f(t), f(t))) {Fm} M1 converges to F : if there exist a parametric representation (λ(t), ξ(t)) of GF and a sequence of parametric representations (λm(t), ξm(t)) of {GFm} respectively, such that lim

n→∞

  • supt |ξm(t) − ξ(t)| + supt |λm(t) − λ(t)|
  • = 0

~ }

slide-25
SLIDE 25

<13Nancy> 25

Theorem 6.2 Let {µm}, µ ∈ P. Then {Fm}

eL

− → F ⇒ {µm}

e

− → µ use Proposition 4.4 µ(g) =

  • R g(x)dF (x) =
  • R g( ˆ

f(t))d f(t) Theorem 6.3 Let {µm}, µ ∈ P and {Fm} be L` evy computable. Then {Fm}

eL

− → F ⇒ µ is computable

~ }

slide-26
SLIDE 26

<13Nancy> 26

Definition 6.4 (Effective d-irrationally pointwise Fine convergence, [15]) Let {Fm}, F be sequentially Fine computable. {Fm} converges effectively d-irrationally pointwise Fine to F

⇌ {Fm(xn)} converges effectively to {F (xn)}

for any Fine computable d-irrational sequence {xn}. Theorem 6.5 ([16]) Let {Fm}, F be sequentially Fine computable. Assume further that F is effectively Fine continuous. Then, effective convergence of {µm} to µ is equivalent to effective d-irrationally pointwise Fine convergence of {Fm} to F . Theorem 6.6 Let {Fm}, F be Fine computable. Then, {Fm}

eL

− → F ⇔ {Fm} effectively d-irrationally pointwise converges to F

~ }

slide-27
SLIDE 27

<13Nancy> 27

7. {gℓ}-metric on P Proposition 7.1 Let {gℓ} be a computable sequence in Cκ and {µm}, {νm} be computable sequences of probability measures. Then, {d{gℓ}(µm, νn)} is computable (double) sequence of reals. Proposition 7.2 Let {gℓ} be an effective separating set. Then, effective convergence is equivalent to effective {gn}-convergence. Theorem 7.3 Let S be the set of all computable sequence of probability

  • measures. Then, (P, d{gℓ}, S) is a metric space with a computability

structure in the sense of (cf. Definition [13]).

~

slide-28
SLIDE 28

<13Nancy> 28

Example 7.4 {gℓ}: an example of an effective separating set in Cκ.

i2`m

(i ` 1)2`m (i + 1)2`m

(tm,i = 2−mum,i Weihrauch) Figure 5: Graph of um,i

an effective enumeration of the set of all finite linear combinations with rational coefficients of {um,i}m∈N,i∈Z. Example 7.5 Example of an effective separating set in (P, d{gℓ}) an effective enumeration of the set of all finite convex combinations with rational coefficients of {δi2−m} Proposition 7.6 Define µm = 2m 22m

i=0

µ(um,−m+i2−m)δ−m+i2−m. Then, {µm} converges to µ. Moreover, {µm} is computable and converges effectively to µ, if µ is computable.

slide-29
SLIDE 29

<13Nancy> 29

Definition 7.7 (Effective compactness, [13]) (X, d, S) is said to be effectively totally bounded if there exist an effective separating set {en} and a recursive function α such that X =

α(p)

  • n=1

BX(en, 2−p) for all p. If (X, d, S) is effectively totally bounded and effectively complete, then we say that (X, d, S) is effectively compact. Proposition 7.8 (P([0, 1]), d{gℓ}) is effectively compact. Weihrauch [27] had used tn,m = 2−nun,m and defined the representation δ′′

  • m. The metric space (P([0, 1]), ρ) is equivalent to (P([0, 1]), dv,{gp}).

He also defined representations δm and δ′

m and proved these

representations are equivalent. (Theorem 4.2 in [27]). By Kamo [8], d{gp}-computability is equivalent to δ′′

m-computability.

~ }

slide-30
SLIDE 30

<13Nancy> 30

8. Comments on Proofs Proposition If {µm} is computable and converges effectively to µ, then µ is computable. The following Lemmas and Proposition are used many times. Lemma 8.1 (Monotone Lemma, [19]) Let {xn,k} be a computable sequence of reals which converges monotonically to {xn} as k tends to infinity for each n. Then, {xn} is computable if and only if the convergence is effective.

✡ ✡ ✡ ✡ ✡ ❏ ❏ ❏ ❏ ❏

1 `n ` 1 `n n n + 1 wn

Figure 6: Graph of wn

slide-31
SLIDE 31

<13Nancy> 31

The following Proposition is fundamental. Proposition 8.2 (Effective tightness of an effectively convergent sequence, Effectivization of Lemma 15.4 in [25]) If a computable {µm} effectively converges to µ, then there exists a recursive function α(k) such that µm(wc

α(k)) < 2−k

for all m. It also hold that µm([−α(k) − 1, α(k) + 1]C) < 2−k for all m.

slide-32
SLIDE 32

<13Nancy> 32

References

[1] Billingsley, P., Convergence of Probability Measures. John Wiley & Sons, Inc. 1968. [2] Deza, M. M. and E. Deza, Encyclopedia of Distances (Chap. 14). Springer, 2009. [3] Edgar, G. A., Integral, Probability and Fractal Measures. Springer, 1998. [4] G´ acs, Peter. Lecture notes on descriptional complexity and randomness. http://www.cs.bu.edu/faculty/gacs/papers/ait-notes.pdf. [5] Gold, E.M., Limiting recursion, JSL, 30-1 (1965), 28-48 [6] Ito, K. An Introduction to Probability Theory. Cambridge University Press, 1978. (Probability Theory. Iwanami Shotenn, 1991. (in Japanese)) [7] Kamo, H. Effective Dini’s Theorem on Effectively Compact Metric Spaces. ENTCS 120, 73-82, 2005. [8] Kamo, H. Private communication. [9] Kolmogorov, A. N. On Skorokhod Convergence. Teor. Prob. Appl., 1 (1956), 215- 222. [10] Lawson, J. D., Domains, integration and ‘positive analysis’, Math. Struct. in com. Science 14(2004), 815-832. [11] L´ evy, P., Th´ eorie de L’Addition Variables Al´

  • eatoires. Gautheir-Villars, 1954.
slide-33
SLIDE 33

<13Nancy> 33

[12] Lu, H. and K. Weihrauch, Computable Riesz Representation for the Dual of C[0; 1], ENTCS 167(2007), 157-177. [13] Mori, T., Y. Tsujii and M. Yasugi. Computability Structures on Metric Spaces. Combinatorics, Complexity and Logic (Proceedings of DMTCS’96), ed. by Bridges et al., 351-362. Springer, 1996. [14] Mori, T., Y. Tsujii and M. Yasugi. Computability of probability distributions and probability distribution functions. Proceedings of the Sixth International Confer- ence on Computability and Complexity in Analysis (DROPS 20), 185-196, 2009. [15] Mori, T., M. Yasugi and Y. Tsujii. Fine convergence of functions and its effec-

  • tivization. Automata, Formal Languages and Algebraic Systems, World Scientific,

2010. [16] Mori, T., Y. Tsujii and M. Yasugi. Fine Computability of Probability Distribu- tion Functions and Computability of Probability Distributions on the Real Line. (manuscript) [17] Mori, T., Y. Tsujii and M. Yasugi. Computability of Probability Distributions and Characteristic Functions. (manuscript) [18] Parthasarathy, K. R. Probability Measures on Metric Spaces. Academic, 1967. [19] Pour-El, M.B. and J. I. Richards. Computability in Analysis and Physics. Springer, 1988. [20] Prokhorov, Yu. V. Convergence of random processes and limit theorems in prob- ability theory. Teor. Prob. Appl., 1 (1956), 157-214.

slide-34
SLIDE 34

<13Nancy> 34

[21] Rachev, S. T., S. V. Stoyanov and F. J. Fabozzi. A Probability Metrics Approach to Financial Risk Measures. Wiley-Blackwell, 2011. [22] Schr¨

  • der, M. and A. Simpson. Representing Probability Measures using Proba-

bilistic Processes. CCA 2005 (Kyoto), 211-226. [23] Schr¨

  • der,

M., Admissible Representations of Probability Measures. ENTCS 167(2007), 61-78. [24] Skorokhod, A. V. Limit theorems for stochastic processes. Teor. Prob. Appl., 1 (1956), 261-290. [25] Tsurumi, S., Probability Theory. Shibunndou, 1964. (in Japanese) [26] Tsujii, Y., M. Yasugi and T. Mori. Some Properties of the Effective Uniform Topological Space. Computability and Complexity in Analysis, (Lecture Notes in Computer Science 2064), ed. by Blanck, J. et al., 336-356. Springer, 2001. [27] Weihrauch, K. Computability on the probability measures on the Borel sets of the unit interval. TCS 219(1999), 421-437. [28] Weihrauch, K., Computable Analysis. Springer, 2000. [29] Yasugi, M., T. Mori and Y. Tsujii. Effective properties of sets and functions in metric spaces with computability structure. Theoretical Computer Science, 219:467-486, 1999. [30] Yasugi, M., T. Mori and Y. Tsujii. Sequential computability of a function -Effective Fine space and limiting recursion-. JUCS 11-12, 2179-2191, 2005.

slide-35
SLIDE 35

<13Nancy> 35

9. Characteristic Functions (CCA2012)

Theorem 9.1 If {µm} is computable, then {ϕm} is uniformly computable. Theorem 9.2 (Effective Glivenko, cf, Theorem 2.6.4 in Ito [6]) Let {ϕm} and ϕ be computable. Then, {µm} converges effectively to µ if {ϕm} converges effectively to ϕ . Theorem 9.3 (Effectivization of Theorem 2.6.3 in Ito [6]) Let {µm} and µ be

  • computable. Then, {ϕm} converges effectively (compact-uniformly) to ϕ if {µm}

converges effectively to µ. Theorem 9.4 {µm} is computable if {ϕm} is computable. Theorem 9.4 is the converse to Theorem 9.1. So, we obtain: Theorem {µm} is computable if and only if {ϕm} is computable. Theorem 9.5 (Effective Bochner’s theorem) In order for ϕ(t) to be a characteristic function of a computable probability measure, it is necessary and sufficient that the following three conditions holds. (i) ϕ is positive definite. (ii) ϕ is computable. (iii) ϕ(0) = 1.

~ }

slide-36
SLIDE 36

<13Nancy> 36

10. De Moivre-Laplace Central Limit Theorem (CCA2012)

The Central Limit Theorem is one of important theorems in probability theory and in statitics. Let (Ω, B, P, {Xm}) be a realization of Coin Tossing (Bernoulli Trials) with success probability p. Theorem 10.1 (Effective de Moivre-Laplace) If p is a computable real number, then the sequence of probability measures of random variables Ym = X1 + · · · Xm − mp

  • mp(1 − p)

=

m

  • ‘=1

X‘ − p √mpq converges effectively to the standard Gaussian probability measure. ψm(t) = E(eitYm) = m

‘=1 E(eit X‘`p

pmpq ) = (pe itpq pmp + qe` itpp pmq )m

By Theorems 10.1 and Theorem 9.2, the following hold: E(f(Ym)) → 1 √ 2π

  • R

f(t)x2` t2

2 dt

effectively if f is bounded computable.

slide-37
SLIDE 37

<13Nancy> 37

11. Graphs

  • ✟✟

  • ✟✟

✟ ✟ ✟

✟ Figure 7:

1 2δ0 + 1 2δ 1

2

slide-38
SLIDE 38

<13Nancy> 38

  • ✟✟

  • ✟✟

✟ ✟✟✟✟ ✟ Figure 8:

1 2δ0 + 1 2δ1

slide-39
SLIDE 39

<13Nancy> 39

  • ✟✟

  • ✟✟

✟ ✟✟ ✟ ✟✟ ✟ Figure 9:

1 2δ0 + 1 2δ2

slide-40
SLIDE 40

<13Nancy> 40 Fm

x ˜ Fm(x) (x; Fm(x)) (zm; ym) (x; F (x))

Fm

x ˜ Fm(x) (x; Fm(x)) (zm; ym) (x; F (x))

Figure 10: Distance between F (x) and Fm(x)

slide-41
SLIDE 41

<13Nancy> 41 xk,ℓ xk,ℓ+1 rk,ℓ rk,ℓ+1 rk,ℓ+2 Ik,ℓ Ik,ℓ+1

F (xk;‘) ` 2`k F (xk;‘) + 2`k F (xk;‘+1) ` 2`k F (xk;‘+1) + 2`k

˜ F (xk,ℓ) ˜ F−(rk,ℓ+1) ˜ F (xk,ℓ+1) ˜ F (rk,ℓ+1) Figure 11: L´ evy convergence

slide-42
SLIDE 42

<13Nancy> 42 h1 h2 h2 Figure 12: Graph of singular G

slide-43
SLIDE 43

<13Nancy> 43

1 d 2`¸(1) 2`¸(2) 2`¸(2) 2 d 2`¸(1) 2`¸(2) 2`¸(3) d1 d2 d3 LF1 LF2

Figure 13: Graph of F and LF

slide-44
SLIDE 44

<13Nancy> 44 U ✲ J1 ✚ ✚ ❃ ❩ ❩ ⑦ M1 J2 ❩ ❩ ⑦ ✚ ✚ ❃ M2 ✁ ✁ ✁ ✁ ✕

uniform convergence at continuity point of x(t)

pointwise convergence at continuity point of x(t)

✻ ❄

convergence of measure

+ suitable condition on

modulus of (dis)continuity

✻ ❄

Ji , Mi - convergence

Figure 14: Relations between Convergences by Skorokhod (case R)

M2 : H-metric between GF and GG M1 + J2 = J1

slide-45
SLIDE 45

<13Nancy> 45

s3 s4 s5 d c b d c b a s1 s2 s3 s4 s5 s6 s7s8

x(t) y(t) ν[a,c]

[0,1][y(t)] = 3,

ν[b,c]

[0,1][y(t)] = 3

ν[a,d]

[0,1] [y(t)] = 0

ν[a,c]

[0,1][x(t)] = ν[b,c] [0,1][x(t)]

= ν[a,d]

[0,1] [x(t)] = 1

Figure 15: Examples of functions in F or D

slide-46
SLIDE 46

<13Nancy> 46

12. Computability on F and F ↔ P

(classical) convergence of {Fm} to F Fm(x) → F (x) for ∀x: point of continuity of F Postulate I (Computability of pdf F ) there is a computable sequence {xn} such that it is dense in R F is effectively continuous at xn: ∃α(n, k) s.t. |y − xn| < 2`¸(n;k) ⇒ |F (y) − F (xn)| < 2`k {F (xn)} is computable Postulate II (Effective convergene of {Fm}) there is a computable sequence {xn} such that it is dense in R each Fm is effectively continuous at xn {Fm(xn)} converges effectively uniformly at {xn} Proposition 12.1 Suppose that there is a computable sequence {xn} such that it is dense in R and {Fm(xn)} is computable. Then {µm} is computable. Outline of the Proof. for a single F , µ and a single f ∈ C». Take α(k) an effective modulus of uniform continuity of f and an integer L such that f(x) = 0 for |x| ≥ L.

slide-47
SLIDE 47

<13Nancy> 47

Let yk;i = −L + i2`¸(k+1) (0 ≤ i ≤ 2L2¸(k+1) + 1). We can find effectively n(k, i) such that xn(k;i) ∈ (yk;i`1, yk;i). Then Sk = 2L2¸(k+1)+1

i=1

f(xn(k;i))(F (xn(k;i)) − F (xn(k;i`1)) converges effectively to

  • R f(x)µ(dx).

|µ(f) − Sk| =

  • 2L2¸(k+1)+1

i=1

  • (yk;i`1;yk;i](f(x) − f(xn(k;i)))µ(dx)
  • < 2`k.

Proposition 12.2 {µm} is computable and {xn} is computable. If Fm is continuous at {xn}, then Fm(xn) is computable.

  • Proof. F (x − 2`n) ≤ µ( ˜

w`

c;n) ≤ F (x) ≤ µ( ˜

w+

c;n) ≤ F (x − 2`n)

{µ( ˜ w`

c;n)} and {µ( ˜

w+

c;n)} are computable

µ( ˜ w`

c;n) − µ( ˜

w+

c;n) ↓ 0 ⇒ convergence is effective

Proposition 12.3 Let LF be computable, x be computable and F be continuous at x. Then F (x) is computable. Outline of the Proof: On a neighborhood of x, ˜ F is strictly increasing and ˆ f(t) = t − f(t) is its inverse. (f = LF ) w`

x;n, w+ x;n, and Proposition 4.4 may be useful.

µ(g) =

  • R g(x)dF (x) =
  • R g( ˆ

f(t))d f(t) for all g ∈ C». Lipschitz continuity implies ≤

  • R g( ˆ

f)(t)dt if g is nonnegative.

slide-48
SLIDE 48

<13Nancy> 48

Take g = w`

x;n, w+ x;n

  • R w`

x;n( ˆ

f(t))d f(t) = µ(w`

x;n) ≤ F (x) ≤ µ(w+ x;n) =

  • R w+

x;n( ˆ

f(t))d f(t) 0 ≤

  • R w+

x;n( ˆ

f(t))d f(t) −

  • R w`

x;n( ˆ

f(t))d f(t) =

  • R{w+

x;n( ˆ

f(t)) − w`

x;n( ˆ

f(t))}d f(t) ≤

  • R{w+

x;n( ˆ

f(t)) − w`

x;n( ˆ

f(t))}dt = ˜

F (x+2`n) ˜ F (x`2`n) {w+ x;n( ˆ

f(t)) − w`

x;n( ˆ

f(t))}dt → 0 effectively (?) Conjecture 12.4 Suppose that there is a computable sequence {xn} such that it is dense in R and {Fm(xn)} is computable. Then {LFm} is computable. If this conjecture is false, then we can say that L´ evy computability is stronger than computability of the correspondint probability measures.

slide-49
SLIDE 49

<13Nancy> 49

Other Facts Fact 12.5 Let {µm} and µ be computable, x be computable, F be continuous at x, and {µm} converge effectively to µ. Then {Fm(x)} converges effectively uniformly to F at x. Or converges effectively continuously at x. Fact 12.6 Let {Fm} and F be L´ evy computable and x be computable. If {LFm} converges effectively to LF then {Fm(x)} converges effectively uniformly to F at x. Need to show

  • R g( ˆ

f(t))d f(t) is computable if g is computable. Put g = un;i. effective convergence ⇔ effective {gn}-convergence ⇒ ?? Fact 12.7 (classical) {Fm} : a sequence of probability distribution functions, {Fm(x)} : converges (pointwise) to F (x) at every continuity point of F Then, {Fm(x)} converges uniformly at every continuity point of F (∵) ∀ǫ > 0 ∃y1 < x < y2 points of continuite of F , |F (y2) − F (y1)| < ǫ Fm(y1) → F (y1), Fm(y2) → F (y2), ∃N s.t. m ≥ N, |Fm(y1) − F (y1)| < ǫ, |Fm(y2) − F (y2)| < ǫ m ≥ N and |y − x| < min{y2 − x, x − y1}, |Fm(y) − F (x)| < 3ǫ

slide-50
SLIDE 50

<13Nancy> 50

Other Conjectures Conjecture 12.8 Let {fm} is a computable sequence in CL and converges monotonically to f ∈ CL. Then the convergence is effective and f is computable. Effective Dini Theorem 12.9 below by Kamo [7] Theorem 12.9 Kamo [7]) Let (X, d, S) be an effectively compact metric space. Let {fn} be a computable sequence of real-valued functions on X and f a computable real-valued function on X. If {fn} converges pointwise monotonically to f as n → ∞, then {fn} converges effectively uniformly to f. Conjecture 12.10 If X is a effectively compact metric space, then (P(X), dv;fgng) is effectively compact for an effective separating set {gn} of P(X). Theorem 12.11 (Effective decomposition of unity, Theorem 3 in [29]) Let {On} be an effective local finite r.e. covering of X. Then there exists a computable sequence of functions {fn} which satisfies (i) fn(x) ≥ 0, (ii) fn(x) = 0 if x / ∈ On, (iii) 1

n=1 fn(x) = 1.

slide-51
SLIDE 51

<13Nancy> 51

13. GADLAC

Definition 13.1 A pair of functions (λ(t), ξ(t)) is said to be a parametric representation of the graph of GF = {(x, z) | F`(x) ≤ z ≤ F (x), x ∈ R} if GF = {(λ(t), ξ(t) | t ∈ R}, ξ(t) is continuous and λ(t) is continuous and monotonically increasing. Definition 13.2 (Skorokhod) {Fm} SM1 converges to F : if there exist a parametric representation (λ(t), ξ(t)) of GF and a sequence of parametric representations (λm(t), ξm(t)) of {GFm} respectively, such that lim

n!1

  • supt |ξm(t) − ξ(t)| + supt |λm(t) − λ(t)|
  • = 0

Postulate (1) A gadlac F is said to be SM1- computable if there exist computable functions (λ(t), ξ(t)) which consist a parametric representation of GF . (computable parametric representation) (2) {Fm} SM1-converges effectively to F : if there exist a computable parametric representation (λ(t), ξ(t)) of GF and a computable sequence of parametric representations (λm(t), ξm(t)) of {GFm} respectively, such that λm(t) and ξm(t)) converges effectively to λ(t) and ξ(t)) respectively.

slide-52
SLIDE 52

<13Nancy> 52

{Fn} J1-converges to F : There exists a sequence of continuous one-to-one and onto mappings {λ(x)} such that lim

n!1 sup x

|Fn(x) − F (λn(x))| = 0, lim

n!1 sup x

|λn(x) − x| = 0. (1) This topology is the well known Skorokhod convergence. A equivalent metric which make complete the space of all gadlaces is discussed in Prokhorov [20] Appendix 1 and in Kolmogorov [9] (see also [1], [18]). Since F (λ(x)) is not continuous, it seems difficult to define the corresponding notion

  • f computability.

{Fm} SJ1 converges effectively to F : if there exist computable one-to-one and onto mapping {λm(t)} such that supt |λm(t) − t| and supt |Fm − F (λm(t))| converge effectively to zero. {Fn} is said to converge uniformly to F at x: For all ǫ > 0, there exists a δ > 0 such that lim sup

n!1

sup

jy`xj<‹

|Fn(x) − F (x)| < ǫ.

slide-53
SLIDE 53

<13Nancy> 53

Well known properties (1) If {Fn} converges uniformly to F at any x in some closed finite interval [a, b], then the convergence is uniform on [a, b]. (2) (λ1(t), ξ1(t)) and (λ2(t), ξ2(t)) are parametric representations of some GF , then there exists a monotonically increasing function u(t) such that λ1(t) = λ2(u(t)) and ξ1(t) = ξ(u(t)).

? not continuous, generalized inverse

(3) If {Fn} converges to F with respect to one of the four convergences, then {Fn} converges uniformly o F at any point of continuity of F . (4) Let Λ be the set of all continuous one-to-one and onto mappings from [0, 1] to [0, 1]. Then, dS(F, G) = inf

  • ǫ
  • ∃λ ∈ Λ, sup

t

|t − λ(t)| < ǫ, sup

t

|F (t) − G(λ(t))| < ǫ

  • (=)

inf

–2Λ

  • sup

t

|F (t) − G(λ(t))| + sup

t

|t − λ(t)|

  • is called the Skorokhod metric. dS-convergence and SJ1-convergence are
  • equivalent. (D([0, 1]), dS) is not complete.

(Skorokhod [24])

slide-54
SLIDE 54

<13Nancy> 54

(5) Let ||λ|| = sups6=t

  • log –(t)`–(s)

t`s

  • and

˜ dS(F, G) = inf

  • ǫ
  • ∃λ, ||λ|| < ǫ, sup

t

|F (t) − G(λ(t))| < ǫ

  • .

Then, dS and ˜ dS are equivalent. (D([0, 1]), ˜ dS) is complete. (Prokhorov [20], Billingsley [1]) (6) Let Fn(t) =      if t < 1

2 − 1 n n 2 (t − 1 2) + 1 2

if

1 2 − 1 n ≤ t < 1 2 + 1 n

1 if t ≥ 1

2 + 1 n

. Then, {Fn} SM1-converges to D 1

2, but the convergence is not SJ2. So, {Fn} does

not SJ1-converges to D 1

2.

Conjecture 13.3 Effective SM1 convergence implies effective uniform convergence at any computable continuity point of F .

slide-55
SLIDE 55

<13Nancy> 55

14. Two dimensional probability measure

Necessary to handle with the Wasserstein Metric Definition 14.1 The L´ evy (L´ evy-Prokhorov) metric ˜ dL(µ, ν) = ˜ dL(F, G) is defined by ˜ dL(F, G) = inf{ǫ > 0 | F (x − ǫ, y − ǫ) − ǫ ≤ G(x, y) ≤ F (x + ǫ, y + ǫ) + ǫ, for all x} GF = {(x, y, z) | F (x`, y`) ≤ z ≤ F (x, y)} ℓs;t = {z + x = s} ∩ {z + y = t}, (s, t, v) = GF ∩ ℓs;t LF (s.t) := v