# <i>L<sub>p</sub></i>-Nested Symmetric Distributions

Fabian Sinz, Matthias Bethge.

Year: 2010, Volume: 11, Issue: 111, Pages: 3409−3451

#### Abstract

In this paper, we introduce a new family of probability densities called *L _{p}*-nested symmetric distributions. The common property, shared by all members of the new class, is the same functional form

*ρ(*, where

**x**) = ~ρ(f(**x**))*f*is a nested cascade of

*L*-norms

_{p}*||*.

**x**||_{p}= (∑ |x_{i}|^{p})^{1/p}*L*-nested symmetric distributions thereby are a special case of

_{p}*ν*-spherical distributions for which

*f*is only required to be positively homogeneous of degree one. While both,

*ν*-spherical and

*L*-nested symmetric distributions, contain many widely used families of probability models such as the Gaussian, spherically and elliptically symmetric distributions,

_{p}*L*-spherically symmetric distributions, and certain types of independent component analysis (ICA) and independent subspace analysis (ISA) models,

_{p}*ν*-spherical distributions are usually computationally intractable. Here we demonstrate that

*L*-nested symmetric distributions are still computationally feasible by deriving an analytic expression for its normalization constant, gradients for maximum likelihood estimation, analytic expressions for certain types of marginals, as well as an exact and efficient sampling algorithm. We discuss the tight links of

_{p}*L*-nested symmetric distributions to well known machine learning methods such as ICA, ISA and mixed norm regularizers, and introduce the nested radial factorization algorithm (NRF), which is a form of non-linear ICA that transforms any linearly mixed, non-factorial

_{p}*L*-nested symmetric source into statistically independent signals. As a corollary, we also introduce the uniform distribution on the

_{p}*L*-nested unit sphere.

_{p}