Growth function

From HandWiki

The growth function, also called the shatter coefficient or the shattering number, measures the richness of a set family or class of function. It is especially used in the context of statistical learning theory, where it is used to study properties of statistical learning methods. The term 'growth function' was coined by Vapnik and Chervonenkis in their 1968 paper, where they also proved many of its properties.[1] It is a basic concept in machine learning.[2] [3]

Definitions

Set-family definition

Let [math]\displaystyle{ H }[/math] be a set family (a set of sets) and [math]\displaystyle{ C }[/math] a set. Their intersection is defined as the following set-family:

[math]\displaystyle{ H\cap C := \{h\cap C\mid h\in H\} }[/math]

The intersection-size (also called the index) of [math]\displaystyle{ H }[/math] with respect to [math]\displaystyle{ C }[/math] is [math]\displaystyle{ |H\cap C| }[/math]. If a set [math]\displaystyle{ C_m }[/math] has [math]\displaystyle{ m }[/math] elements then the index is at most [math]\displaystyle{ 2^m }[/math]. If the index is exactly 2m then the set [math]\displaystyle{ C }[/math] is said to be shattered by [math]\displaystyle{ H }[/math], because [math]\displaystyle{ H\cap C }[/math] contains all the subsets of [math]\displaystyle{ C }[/math], i.e.:

[math]\displaystyle{ |H\cap C|=2^{|C|}, }[/math]

The growth function measures the size of [math]\displaystyle{ H\cap C }[/math] as a function of [math]\displaystyle{ |C| }[/math]. Formally:

[math]\displaystyle{ \operatorname{Growth}(H,m) := \max_{C: |C|=m} |H\cap C| }[/math]

Hypothesis-class definition

Equivalently, let [math]\displaystyle{ H }[/math] be a hypothesis-class (a set of binary functions) and [math]\displaystyle{ C }[/math] a set with [math]\displaystyle{ m }[/math] elements. The restriction of [math]\displaystyle{ H }[/math] to [math]\displaystyle{ C }[/math] is the set of binary functions on [math]\displaystyle{ C }[/math] that can be derived from [math]\displaystyle{ H }[/math]:[3]:45

[math]\displaystyle{ H_{C} := \{(h(x_1),\ldots,h(x_m))\mid h\in H, x_i\in C\} }[/math]

The growth function measures the size of [math]\displaystyle{ H_C }[/math] as a function of [math]\displaystyle{ |C| }[/math]:[3]:49

[math]\displaystyle{ \operatorname{Growth}(H,m) := \max_{C: |C|=m} |H_C| }[/math]

Examples

1. The domain is the real line [math]\displaystyle{ \mathbb{R} }[/math]. The set-family [math]\displaystyle{ H }[/math] contains all the half-lines (rays) from a given number to positive infinity, i.e., all sets of the form [math]\displaystyle{ \{x \gt x_0 \mid x\in \mathbb{R}\} }[/math] for some [math]\displaystyle{ x_0\in \mathbb{R} }[/math]. For any set [math]\displaystyle{ C }[/math] of [math]\displaystyle{ m }[/math] real numbers, the intersection [math]\displaystyle{ H\cap C }[/math] contains [math]\displaystyle{ m+1 }[/math] sets: the empty set, the set containing the largest element of [math]\displaystyle{ C }[/math], the set containing the two largest elements of [math]\displaystyle{ C }[/math], and so on. Therefore: [math]\displaystyle{ \operatorname{Growth}(H,m)=m+1 }[/math].[1]:Ex.1 The same is true whether [math]\displaystyle{ H }[/math] contains open half-lines, closed half-lines, or both.

2. The domain is the segment [math]\displaystyle{ [0,1] }[/math]. The set-family [math]\displaystyle{ H }[/math] contains all the open sets. For any finite set [math]\displaystyle{ C }[/math] of [math]\displaystyle{ m }[/math] real numbers, the intersection [math]\displaystyle{ H\cap C }[/math] contains all possible subsets of [math]\displaystyle{ C }[/math]. There are [math]\displaystyle{ 2^m }[/math] such subsets, so [math]\displaystyle{ \operatorname{Growth}(H,m)=2^m }[/math]. [1]:Ex.2

3. The domain is the Euclidean space [math]\displaystyle{ \mathbb{R}^n }[/math]. The set-family [math]\displaystyle{ H }[/math] contains all the half-spaces of the form: [math]\displaystyle{ x\cdot \phi \geq 1 }[/math], where [math]\displaystyle{ \phi }[/math] is a fixed vector. Then [math]\displaystyle{ \operatorname{Growth}(H,m)=\operatorname{Comp}(n,m) }[/math], where Comp is the number of components in a partitioning of an n-dimensional space by m hyperplanes.[1]:Ex.3

4. The domain is the real line [math]\displaystyle{ \mathbb{R} }[/math]. The set-family [math]\displaystyle{ H }[/math] contains all the real intervals, i.e., all sets of the form [math]\displaystyle{ \{x \in [x_0,x_1] | x\in \mathbb{R}\} }[/math] for some [math]\displaystyle{ x_0,x_1\in \mathbb{R} }[/math]. For any set [math]\displaystyle{ C }[/math] of [math]\displaystyle{ m }[/math] real numbers, the intersection [math]\displaystyle{ H\cap C }[/math] contains all runs of between 0 and [math]\displaystyle{ m }[/math] consecutive elements of [math]\displaystyle{ C }[/math]. The number of such runs is [math]\displaystyle{ {m+1 \choose 2}+1 }[/math], so [math]\displaystyle{ \operatorname{Growth}(H,m)={m+1 \choose 2}+1 }[/math].

Polynomial or exponential

The main property that makes the growth function interesting is that it can be either polynomial or exponential - nothing in-between.

The following is a property of the intersection-size:[1]:Lem.1

  • If, for some set [math]\displaystyle{ C_m }[/math] of size [math]\displaystyle{ m }[/math], and for some number [math]\displaystyle{ n\leq m }[/math], [math]\displaystyle{ |H\cap C_m| \geq \operatorname{Comp}(n,m) }[/math] -
  • then, there exists a subset [math]\displaystyle{ C_n\subseteq C_m }[/math] of size [math]\displaystyle{ n }[/math] such that [math]\displaystyle{ |H\cap C_n| = 2^n }[/math].

This implies the following property of the Growth function.[1]:Th.1 For every family [math]\displaystyle{ H }[/math] there are two cases:

  • The exponential case: [math]\displaystyle{ \operatorname{Growth}(H,m) = 2^m }[/math] identically.
  • The polynomial case: [math]\displaystyle{ \operatorname{Growth}(H,m) }[/math] is majorized by [math]\displaystyle{ \operatorname{Comp}(n,m) \leq m^n+1 }[/math], where [math]\displaystyle{ n }[/math] is the smallest integer for which [math]\displaystyle{ \operatorname{Growth}(H,n) \lt 2^n }[/math].

Other properties

Trivial upper bound

For any finite [math]\displaystyle{ H }[/math]:

[math]\displaystyle{ \operatorname{Growth}(H,m) \leq |H| }[/math]

since for every [math]\displaystyle{ C }[/math], the number of elements in [math]\displaystyle{ H\cap C }[/math] is at most [math]\displaystyle{ |H| }[/math]. Therefore, the growth function is mainly interesting when [math]\displaystyle{ H }[/math] is infinite.

Exponential upper bound

For any nonempty [math]\displaystyle{ H }[/math]:

[math]\displaystyle{ \operatorname{Growth}(H,m) \leq 2^m }[/math]

I.e, the growth function has an exponential upper-bound.

We say that a set-family [math]\displaystyle{ H }[/math] shatters a set [math]\displaystyle{ C }[/math] if their intersection contains all possible subsets of [math]\displaystyle{ C }[/math], i.e. [math]\displaystyle{ H\cap C = 2^C }[/math]. If [math]\displaystyle{ H }[/math] shatters [math]\displaystyle{ C }[/math] of size [math]\displaystyle{ m }[/math], then [math]\displaystyle{ \operatorname{Growth}(H,C)=2^{m} }[/math], which is the upper bound.

Cartesian intersection

Define the Cartesian intersection of two set-families as:

[math]\displaystyle{ H_1\bigotimes H_2 := \{h_1\cap h_2 \mid h_1\in H_1, h_2\in H_2\} }[/math].

Then:[2]:57

[math]\displaystyle{ \operatorname{Growth}(H_1\bigotimes H_2, m) \leq \operatorname{Growth}(H_1, m)\cdot \operatorname{Growth}(H_2, m) }[/math]

Union

For every two set-families:[2]:58

[math]\displaystyle{ \operatorname{Growth}(H_1\cup H_2, m) \leq \operatorname{Growth}(H_1, m) + \operatorname{Growth}(H_2, m) }[/math]

VC dimension

The VC dimension of [math]\displaystyle{ H }[/math] is defined according to these two cases:

  • In the polynomial case, [math]\displaystyle{ \operatorname{VCDim}(H) = n-1 }[/math] = the largest integer [math]\displaystyle{ d }[/math] for which [math]\displaystyle{ \operatorname{Growth}(H,d) = 2^d }[/math].
  • In the exponential case [math]\displaystyle{ \operatorname{VCDim}(H) = \infty }[/math].

So [math]\displaystyle{ \operatorname{VCDim}(H)\geq d }[/math] if-and-only-if [math]\displaystyle{ \operatorname{Growth}(H,d)=2^d }[/math].

The growth function can be regarded as a refinement of the concept of VC dimension. The VC dimension only tells us whether [math]\displaystyle{ \operatorname{Growth}(H,d) }[/math] is equal to or smaller than [math]\displaystyle{ 2^d }[/math], while the growth function tells us exactly how [math]\displaystyle{ \operatorname{Growth}(H,m) }[/math] changes as a function of [math]\displaystyle{ m }[/math].

Another connection between the growth function and the VC dimension is given by the Sauer–Shelah lemma:[3]:49

If [math]\displaystyle{ \operatorname{VCDim}(H)=d }[/math], then:
for all [math]\displaystyle{ m }[/math]: [math]\displaystyle{ \operatorname{Growth}(H,m)\leq \sum_{i=0}^{d} {m\choose i} }[/math]

In particular,

for all [math]\displaystyle{ m\gt d+1 }[/math]: [math]\displaystyle{ \operatorname{Growth}(H,m)\leq (em/d)^d = O(m^d) }[/math]
so when the VC dimension is finite, the growth function grows polynomially with [math]\displaystyle{ m }[/math].

This upper bound is tight, i.e., for all [math]\displaystyle{ m\gt d }[/math] there exists [math]\displaystyle{ H }[/math] with VC dimension [math]\displaystyle{ d }[/math] such that:[2]:56

[math]\displaystyle{ \operatorname{Growth}(H,m)= \sum_{i=0}^{d} {m\choose i} }[/math]

Entropy

While the growth-function is related to the maximum intersection-size, the entropy is related to the average intersection size:[1]:272–273

[math]\displaystyle{ \operatorname{Entropy}(H, m) = E_{|C_m|=m}\big[\log_2(|H\cap C_m|)\big] }[/math]

The intersection-size has the following property. For every set-family [math]\displaystyle{ H }[/math]:

[math]\displaystyle{ |H\cap (C_1 \cup C_2)| \leq |H\cap C_1|\cdot |H \cap C_2| }[/math]

Hence:

[math]\displaystyle{ \operatorname{Entropy}(H, m_1+m_2) \leq \operatorname{Entropy}(H, m_1) + \operatorname{Entropy}(H, m_2) }[/math]

Moreover, the sequence [math]\displaystyle{ \operatorname{Entropy}(H, m)/m }[/math] converges to a constant [math]\displaystyle{ c\in [0,1] }[/math] when [math]\displaystyle{ m\to\infty }[/math].

Moreover, the random-variable [math]\displaystyle{ \log_2{|H\cap C_m|/m} }[/math] is concentrated near [math]\displaystyle{ c }[/math].

Applications in probability theory

Let [math]\displaystyle{ \Omega }[/math] be a set on which a probability measure [math]\displaystyle{ \Pr }[/math] is defined. Let [math]\displaystyle{ H }[/math] be family of subsets of [math]\displaystyle{ \Omega }[/math] (= a family of events).

Suppose we choose a set [math]\displaystyle{ C_m }[/math] that contains [math]\displaystyle{ m }[/math] elements of [math]\displaystyle{ \Omega }[/math], where each element is chosen at random according to the probability measure [math]\displaystyle{ P }[/math], independently of the others (i.e., with replacements). For each event [math]\displaystyle{ h\in H }[/math], we compare the following two quantities:

  • Its relative frequency in [math]\displaystyle{ C_m }[/math], i.e., [math]\displaystyle{ |h\cap C_m|/m }[/math];
  • Its probability [math]\displaystyle{ \Pr[h] }[/math].

We are interested in the difference, [math]\displaystyle{ D(h,C_m) := \big||h\cap C_m|/m - \Pr[h]\big| }[/math]. This difference satisfies the following upper bound:

[math]\displaystyle{ \Pr\left[\forall h\in H: D(h,C_m) \leq \sqrt{8(\ln\operatorname{Growth}(H, 2m) + \ln(4/\delta)) \over m} \right] ~~~~ \gt ~~~~ 1 - \delta }[/math]

which is equivalent to:[1]:Th.2

[math]\displaystyle{ \Pr\big[\forall h\in H: D(h,C_m) \leq \varepsilon\big] ~~~~ \gt ~~~~ 1 - 4\cdot \operatorname{Growth}(H, 2m)\cdot \exp(-\varepsilon^2\cdot m/8) }[/math]

In words: the probability that for all events in [math]\displaystyle{ H }[/math], the relative-frequency is near the probability, is lower-bounded by an expression that depends on the growth-function of [math]\displaystyle{ H }[/math].

A corollary of this is that, if the growth function is polynomial in [math]\displaystyle{ m }[/math] (i.e., there exists some [math]\displaystyle{ n }[/math] such that [math]\displaystyle{ \operatorname{Growth}(H,m)\leq m^n+1 }[/math]), then the above probability approaches 1 as [math]\displaystyle{ m\to\infty }[/math]. I.e, the family [math]\displaystyle{ H }[/math] enjoys uniform convergence in probability.

References

  1. 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 Vapnik, V. N.; Chervonenkis, A. Ya. (1971). "On the Uniform Convergence of Relative Frequencies of Events to Their Probabilities". Theory of Probability & Its Applications 16 (2): 264. doi:10.1137/1116025.  This is an English translation, by B. Seckler, of the Russian paper: "On the Uniform Convergence of Relative Frequencies of Events to Their Probabilities". Dokl. Akad. Nauk 181 (4): 781. 1968.  The translation was reproduced as: Vapnik, V. N.; Chervonenkis, A. Ya. (2015). "On the Uniform Convergence of Relative Frequencies of Events to Their Probabilities". Measures of Complexity. pp. 11. doi:10.1007/978-3-319-21852-6_3. ISBN 978-3-319-21851-9. 
  2. 2.0 2.1 2.2 2.3 Mohri, Mehryar; Rostamizadeh, Afshin; Talwalkar, Ameet (2012). Foundations of Machine Learning. USA, Massachusetts: MIT Press. ISBN 9780262018258. , especially Section 3.2
  3. 3.0 3.1 3.2 3.3 Shalev-Shwartz, Shai; Ben-David, Shai (2014). Understanding Machine Learning – from Theory to Algorithms. Cambridge University Press. ISBN 9781107057135.