Neighbourhoods of Randomness and Independence

Similar documents
Information geometry for bivariate distribution control

AN AFFINE EMBEDDING OF THE GAMMA MANIFOLD

Information Geometry of Bivariate Gamma Exponential Distributions

Topics in Information Geometry. Dodson, CTJ. MIMS EPrint: Manchester Institute for Mathematical Sciences School of Mathematics

Information geometry of the power inverse Gaussian distribution

Universal connection and curvature for statistical manifold geometry

On the entropy flows to disorder

On the entropy flows to disorder

arxiv: v2 [math-ph] 22 May 2009

Universal connection and curvature for statistical manifold geometry

On the entropy flows to disorder

Some illustrations of information geometry in biology and physics

On random sequence spacing distributions

Random Variables and Their Distributions

Including colour, pdf graphics and hyperlinks

Lecture Notes in Mathematics Editors: J.-M. Morel, Cachan F. Takens, Groningen B. Teissier, Paris

Chapter 5 continued. Chapter 5 sections

LECTURE: KOBORDISMENTHEORIE, WINTER TERM 2011/12; SUMMARY AND LITERATURE

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =

Statistics Ph.D. Qualifying Exam: Part I October 18, 2003

Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011

Chapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves

3.2 Frobenius Theorem

4.7 The Levi-Civita connection and parallel transport

Example 2 (new version): the generators are. Example 4: the generators are

Orientation transport

INTERSECTION THEORY. Note that. If we identify the bidual V with V in a canonical fashion, then we deduce (U # ) # = U.

Open Questions in Quantum Information Geometry

Bordism and the Pontryagin-Thom Theorem

Information Geometric view of Belief Propagation

Introduction to Information Geometry

Changing coordinates to adapt to a map of constant rank

Chapter 5. Chapter 5 sections

Comparing Distance Measures Between Bivariate Gamma Processes

0.1 Complex Analogues 1

A CONSTRUCTION OF TRANSVERSE SUBMANIFOLDS

LECTURE 2. (TEXED): IN CLASS: PROBABLY LECTURE 3. MANIFOLDS 1. FALL TANGENT VECTORS.

History (Minimal Dynamics Theorem, Meeks-Perez-Ros 2005, CMC Dynamics Theorem, Meeks-Tinaglia 2008) Giuseppe Tinaglia

1 The classification of Monotone sets in the Heisenberg group

We have the following immediate corollary. 1

ON r-neighbourly SUBMANIFOLDS IN R N. Victor A. Vassiliev. 1. Introduction. Let M be a k-dimensional manifold, and r a natural number.

STA 4322 Exam I Name: Introduction to Statistics Theory

Lecture 3. Probability - Part 2. Luigi Freda. ALCOR Lab DIAG University of Rome La Sapienza. October 19, 2016

HILBERT S THEOREM ON THE HYPERBOLIC PLANE

Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds

The existence of light-like homogeneous geodesics in homogeneous Lorentzian manifolds. Sao Paulo, 2013

Information theoretic analysis of protein sequences shows that amino acids self cluster

Geometric inequalities for black holes

Geometry 9: Serre-Swan theorem

ON THE FOLIATION OF SPACE-TIME BY CONSTANT MEAN CURVATURE HYPERSURFACES

Cayley Graphs of Finitely Generated Groups

Unsupervised Learning Techniques Class 07, 1 March 2006 Andrea Caponnetto

Qualifying Exams I, Jan where µ is the Lebesgue measure on [0,1]. In this problems, all functions are assumed to be in L 1 [0,1].

arxiv: v4 [math.dg] 10 Feb 2018

Available at: off IC/2001/96 United Nations Educational Scientific and Cultural Organization and International Atomic

Poisson geometry of b-manifolds. Eva Miranda

Foliations of Three Dimensional Manifolds

Finite Presentations of Hyperbolic Groups

zi z i, zi+1 z i,, zn z i. z j, zj+1 z j,, zj 1 z j,, zn

Hyperbolic Geometry on Geometric Surfaces

EXOTIC SMOOTH STRUCTURES ON TOPOLOGICAL FIBRE BUNDLES B

ELEMENTS OF PROBABILITY THEORY

Visual Tracking via Geometric Particle Filtering on the Affine Group with Optimal Importance Functions

Contents. O-minimal geometry. Tobias Kaiser. Universität Passau. 19. Juli O-minimal geometry

Fisher Information in Gaussian Graphical Models

Math 6455 Oct 10, Differential Geometry I Fall 2006, Georgia Tech. Integration on Manifolds, Volume, and Partitions of Unity

MATH DIFFERENTIAL GEOMETRY. Contents

LECTURE 28: VECTOR BUNDLES AND FIBER BUNDLES

Micro-support of sheaves

arxiv: v1 [math.gt] 23 Apr 2014

Cup product and intersection

LECTURE 9: THE WHITNEY EMBEDDING THEOREM

The Minimal Element Theorem

GEOMETRY OF GEODESIC SPHERES IN A COMPLEX PROJECTIVE SPACE IN TERMS OF THEIR GEODESICS

Representation of Markov chains. Stochastic stability AIMS 2014

Notes on quotients and group actions

t x 1 e t dt, and simplify the answer when possible (for example, when r is a positive even number). In particular, confirm that EX 4 = 3.

CONSTRAINED PERCOLATION ON Z 2

0.1 Spec of a monoid

This exam is closed book and closed notes. (You will have access to a copy of the Table of Common Distributions given in the back of the text.

Lipschitz matchbox manifolds

1 )(y 0) {1}. Thus, the total count of points in (F 1 (y)) is equal to deg y0

Lecture 1: August 28

CS Lecture 19. Exponential Families & Expectation Propagation

MFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015

Differential Geometry, Lie Groups, and Symmetric Spaces

Lecture 8. Connections

1 Contact structures and Reeb vector fields

Multivariate Distributions

Dimensionality reduction for classification of stochastic fibre radiographs

Biconservative surfaces in Riemannian manifolds

Advanced Machine Learning & Perception

PICARD S THEOREM STEFAN FRIEDL

Diffusion/Inference geometries of data features, situational awareness and visualization. Ronald R Coifman Mathematics Yale University

On Consistent Hypotheses Testing

Dimensionality Reduction for Information Geometric Characterization of Surface Topographies. Dodson, CTJ and Mettanen, M and Sampson, WW

GSI Geometric Science of Information, Paris, August Dimensionality reduction for classification of stochastic fibre radiographs

Transversality. Abhishek Khetan. December 13, Basics 1. 2 The Transversality Theorem 1. 3 Transversality and Homotopy 2

Weighted Exponential Distribution and Process

Algebraic Topology European Mathematical Society Zürich 2008 Tammo tom Dieck Georg-August-Universität

Transcription:

Neighbourhoods of Randomness and Independence C.T.J. Dodson School of Mathematics, Manchester University Augment information geometric measures in spaces of distributions, via explicit geometric representations of neighbourhoods for these important states for stochastic processes: randomness, independence, uniformity. Significant theoretically because very general, and practically because topological so stable under perturbations.

Proposition Every neighbourhood of a random process contains a neighbourhood of stochastic processes subordinate to gamma distributions. Proof Dodson and Matsuzoe provided affine immersion in R 3 for G, manifold of gamma pdfs. Natural coordinates (µ = α/β, α) Then G is the graph of the affine immersion {h, ξ} where ξ is a transversal vector field along h h : G R 3 : ( ) µ α µ α log Γ(α) α log µ Submanifold of exponential pdfs is curve, ξ =. (, ) R 3 : µ {µ,, log µ }

Tubular neighbourhood in R 3.6 cos θ.6 µ cos θ {µ,.6 sin θ, log µ} + µ + µ θ [, π) contains all immersions for small enough perturbations of exponential distributions. This tubular neighbourhood intersects with the gamma manifold immersion to yield the required neighbourhood of gamma distributions.

µ 3 - α 3 Gamma manifold affine immersion in natural coordinates µ = α/β, α as a surface in R 3. Tubular neighbourhood surrounds all exponential distributions these lie on the curve α = in the surface.

α 3 3 Continuous image of gamma manifold affine immersion as surface in R 3 using standard coordinates. β - - Tubular neighbourhood surrounds all exponential distributions these lie on the curve α = in the surface.

Log-gamma manifold L and neighbourhoods of uniformity The pdfs for random variable N [, ] given by g(n, µ, β) = for µ > and β >. β µ N ( β µ )β (log N )β Γ(β) Uniform distribution lim g(n, µ, ) = g(n,, ) =. µ

Log-gamma pdfs g(n; µ, β), N [, ], with central mean < N >=.5, and β =.5,,, 5. Cases β < correspond in gamma distributions to clustering in an underlying spatial process; conversely, β > corresponds to dispersion and greater evenness than random. g(n; µ, β).5 β = 5.5.5 β = β = β =.5...6.8 N

Family of log-gamma pdfs p ν (N), N [, ] with mean 3 8.. 6 ν N.6.8 Uniform distribution for ν = : p (N) =, N ie random (Poisson) case; clustering ν < ; smoothed ν >.

Log-gamma manifold L L isometric with gamma manifold, G. Hence, immersion of G in R 3 above, represents also the log-gamma manifold L. Then, since the isometry sends the exponential distribution to the uniform distribution on [, ], we obtain a general deduction Proposition Every neighbourhood of the uniform distribution contains a neighbourhood of log-gamma distributions. Equivalently, Proposition Every neighbourhood of a uniform stochastic process contains a neighbourhood of stochastic processes subordinate to log-gamma distributions.

Freund manifold F and neighbourhoods of independence Let F be the manifold of Freund bivariate mixture exponential distributions, so with positive parameters α i, β i, F {f f(x, y; α, β, α, β ) = { α β e β y (α +α β )x for x < y α β e β x (α +α β )y for y x }. Submanifold F F : α = α, β = β The pdfs are of form : f(x, y; α, β ) = { α β e β y ( α β )x for < x < y α β e β x ( α β )y for < y < x Parameters α, β >. The covariance, correlation coefficient and marginal pdfs, for X and Y are given by : Cov(X, Y ) = ( α ), β ρ(x, Y ) = α 3 α + β,

f X (x) = ( ) α α β f Y (y) = ( ) α α β with x, y β e β x + β e β y + ( α β α β ( α β α β ) ) ( α ) e α x ( α ) e α y Independence case ρ(x, Y ) = = σ when α = β F an exponential family with parameters (α, β ) and potential function ψ = log(α β )

Proposition In the affine embedding of the Freund submanifold F in R 3, a tubular neighbourhood of the curve α = β will contain all immersions of bivariate exponential processes sufficiently close to the case of independence.

α - - β Affine immersion in natural coordinates (α, β ) as a surface in R 3 for the Freund submanifold F. Tubular neighbourhood surrounds the curve (α = β in the surface) consisting of all bivariate distributions having identical exponential marginals and zero covariance.

Neighbourhoods of independence for Gaussian processes Proposition Let N be the bivariate Gaussian manifold with the Fisher metric g and the exponential connection (). Denote by (θ i ) a natural coordinate system. Then N can be realized in R 6 by the affine immersion {f, ξ} f : Θ R 6 : [ [ ] ] θi θ i, ξ = ϕ(θ) with potential function ϕ(θ) = log( π ) ( θ θ 3 θ θ θ + θ ) θ 5.,

Proposition The bivariate Gaussian 5-manifold admits a -dimensional submanifold through which can be provided a neighbourhood of bivariate distributions having identical Gaussian marginals and zero covariance containing the independent case. Outline proof Bivariate Gaussians with zero means (µ = µ = ) and identical standard deviation σ = σ = σ is represented by the surface N R 3 : (θ, θ ) (θ, θ, ϕ(θ)), ξ = (,, ). where (θ, θ ) = ( σ, σ ) = σ σ ϕ(θ) = log( π σ).

Affine immersion for Gaussian processes Submanifold of the independent case with zero means and identical standard deviation is curve (, ) R 3 : ( ) (,, log( π σ)) σ σ Tubular neighbourhood represents departures from independence.

Proposition In the affine immersion as a surface in R 3 for the bivariate Gaussian distributions with zero means and identical standard deviation σ, a tubular neighbourhood of the case of zero covariance will contain all immersions of bivariate Gaussian processes sufficiently close to the independence case. Figure shows affine embedding of -submanifold as a surface in R 3, and an R 3 -tubular neighbourhood, of the curve σ = in the surface.

Affine immersion in natural coordinates (θ, θ ) = ( σ, σ ) as surface in R3 for bivariate Gaussian distributions with zero means and identical standard deviation σ. Tubular neighbourhood surrounds curve (σ = in the surface) representing all bivariate distributions having identical Gaussian marginals and zero covariance. Tubular neighourhood contains all small departures from independence. θ - - -3 -.5 -.5 - -.5 θ.5

Continuous image of the affine immersion as a surface in R 3 using standard coordinates for the bivariate Gaussian distributions with zero means and identical standard deviation σ. Tubular neighbourhood surrounds curve of independent states, (σ =,) and contains all small departures from independence. σ 3 3 - -.5 σ.5

Results Explicit representations in R 3 of tubular neighbourhoods that provide the following: A. All processes sufficiently close to a Poisson process B. All processes sufficiently close to a uniform process C. All bivariate processes sufficiently close to the independent bivariate Poisson 3 process D. All bivariate processes sufficiently close to the independent bivariate Gaussian process In an information geometric immersion Via unique associated exponential process 3 Marginals having same mean Marginals both N(, σ)