Publications · Conference

Approximation by Finite Mixtures of Continuous Density Functions That Vanish at Infinity

TrungTin Nguyen, Hien D. Nguyen, Faicel Chamroukhi, Geoffrey J. McLachlan

Cogent Math. & Stats. · Journal Cogent Mathematics & Statistics. Open-access journal article (Vol. 7, 2020).

Abstract

Given sufficiently many components, it is often cited that finite mixture models can approximate any other probability density function (pdf) to an arbitrary degree of accuracy. Unfortunately, the nature of this approximation result is often left unclear. We prove that finite mixture models constructed from pdfs in 𝒞0 can be used to conduct approximation of various classes of approximands in a number of different modes. That is, we prove approximands in 𝒞0 can be uniformly approximated, approximands in 𝒞b can be uniformly approximated on compact sets, and approximands in p can be approximated with respect to the p, for p[1,). Furthermore, we also prove that measurable functions can be approximated, almost everywhere.

Keywords— Approximation theory, probability density functions, finite mixture models, Riemann summation, uniform approximation.

1 Introduction

Let x be an element in the Euclidean space, defined by n and the norm 2, for some n. Let f:n be a function, such that f0, everywhere, and fdλ=1, where λ is the Lesbegue measure. We say that f is a probability density function (pdf) on the domain n (an expression that we will drop, from hereon in). Let g:n be another pdf, and for each m, define the functional class:

mg = {h:h(x)=i=1mci1σing(xμiσi),μinσi+c𝕊m1i[m]},

where c=(c1,,cm), +=(0,),

𝕊m1={cm:i=1mci=1 and ci0,i[m]},

[m]={1,,m}, and () is the matrix transposition operator. We say that any hmg is a m-component location-scale finite mixture of the pdf g.

The study of pdfs in the class mg is an evergreen area of applied and technical research, in statistics. We point the interested reader to the many comprehensive books on the topic, such as [10],[35], [22], [20], [23], [14], [33], [24], and [15].

Much of the popularity of finite mixture models stem from the folk theorem, which states that for any density f, there exists an hmg, for some sufficiently large number of components m, such that h approximates f arbitrarily closely, in some sense. Examples of this folk theorem come in statements such as: “provided the number of component densities is not bounded above, certain forms of mixture can be used to provide arbitrarily close approximation to a given probability distribution” [35, p. 50], “the [mixture] model forms can fit any distribution and significantly increase model fit” [37, p. 173], and “a mixture model can approximate almost any distribution” [39, p. 500]. Other statements conveying the same sentiment are reported in [28]. There is a sense of vagary in the reported statements, and little is ever made clear regarding the technical nature of the folk theorem.

In order to proceed, we require the following definitions. We say that f is compactly supported on 𝕂n, if 𝕂 is compact and if 𝟏𝕂f=0, where 𝟏𝕏 is the indicator function that takes value 1 when x𝕏 and 0, elsewhere, and () is the set complement operator (i.e., 𝕏=n\𝕏). Here, 𝕏 is a generic subset of n. Furthermore, we say that fp(𝕏) for any 1p<, if

fp(𝕏)=(|𝟏𝕏f|pdλ)1/p<,

and for p=, if

f(𝕏)=inf{a0:λ({x𝕏:|f(x)|>a})=0}<,

where we call p(𝕏) the p-norm on 𝕏. When 𝕏=n, we shall write p(n)=p. In addition, we define the so-called Kullback-Leibler divergence, see [17], between any two pdfs f and g on 𝕏 as

KL𝕏(f,g)=𝟏𝕏flog(fg)dλ.

In [28], the approximation of pdfs f by the class mg was explored in a restrictive setting. Let {hmg} be a sequence of functions that draw elements from the nested sequence of sets {mg} (i.e., h1g1g,h2g2g,). The following result of [40] was presented in [28], along with a collection of its implications, such as the results of from [19] and [30].

Theorem 1 (Zeevi and Meir, 1997).

If

f{f:𝟏𝕂fββ>0}2(𝕂)

and g are pdfs and 𝕂 is compact, then there exists a sequence {hmg} such that

limmfhmg2(𝕂)=0 and limmKL𝕂(f,hmg)=0.

Although powerful, this result is restrictive in the sense that it only permits approximation in the 2 norm on compact sets 𝕂, and that the result only allows for approximation of functions f that are strictly positive on 𝕂. In general, other modes of approximation are desirable, in particular approximation in p-norm for p=1 or p= are of interest, where the latter case is generally referred to as uniform approximation. Furthermore, the strict-positivity assumption, and the restriction on compact sets limits the scope of applicability of Theorem 1. An example of an interesting application of extensions beyond Theorem 1 is within the 1-norm approximation framework of [8].

Let g:n again be a pdf. Then, for each m, we define

𝒩mg ={h:h(x)=i=1mci1σing(xμiσi)μinσi+cii[m]},

which we call the set of m-component location-scale linear combinations of the pdf g. In the past, results regarding approximations of pdfs f via functions η𝒩mg have been more forthcoming. For example, in the case of g=ϕ, where

ϕ(x)=(2π)n/2exp(x22/2), (1)

is the standard normal pdf. Denoting the class of continuous functions with support on n by 𝒞. We have the result that for every pdf f, compact set 𝕂n, and ϵ>0, there exists an m and h𝒩mϕ, such that fh(𝕂)<ϵ [32, Lem. 1]. Furthermore, upon defining the set of continuous functions that vanish at infinity by

𝒞0 ={f𝒞:ϵ>0, a compact 𝕂n, such thatf(𝕂)<ϵ},

we also have the result: for every pdf f𝒞0 and ϵ>0, there exists an m and h𝒩mϕ, such that fh<ϵ [32, Thm. 2]. Both of the results from [32] are simple implications of the famous Stone-Weierstrass theorem (cf. [34] and [7]).

To the best of our knowledge, the strongest available claim that is made regarding the folk theorem, within a probabilistic or statistical context, is that of [6, Thm. 33.2]. Let {ηmg} be a sequence of functions that draw elements from the nested sequence of sets {𝒩mg}, in the same manner as {hmg}. We paraphrase the claim without loss of fidelity, as follows.

Claim 1.

If f,g𝒞 are pdfs and 𝕂n is compact, then there exists a sequence {ηmg}, such that

limmfηmg(𝕂)=0.

Unfortunately, the proof of Claim 1 is not provided within [6]. The only reference of the result is to an undisclosed location in [4], which, upon investigation, can be inferred to be Theorem 5 of [4, Ch. 20]. It is further notable that there is no proof provided for the theorem. Instead, it is stated that the proof is similar to that of Theorem 1 in [4, Ch. 24], which is a reproduction of the proof for [38, Lem. 3.1].

There is a major problem in applying the proof technique of [38, Lem. 3.1] in order to prove Claim 1. The proof of [38, Lem. 3.1] critically depends upon the statement that “there is no loss of generality in assuming that f(x)=0 for xn\2𝕂”. Here, for a+, a𝕂={xn:x=ayy𝕂}. The assumption is necessary in order to write any convolution with f and an arbitrary continuous function as an integral over a compact domain, and then to use a Riemann sum to approximate such an integral. Subsequently, such a proof technique does not work outside the class of continuous functions that are compactly supported on a𝕂. Thus, one cannot verify Claim 1 from the materials of [38], [4], and [6], alone.

Some recent results in the spirit of Claim 1 have been obtained by [27] and [26], using methods from the study of universal series (see for example in [25]).

Let

𝒲={f𝒞0:ynsupx[0,1]n |f(x+y)|<}

denote the so-called Wiener’s algebra (see, e.g., [11]) and let

𝒱 ={f𝒞0:xn|f(x)|β(1+x2)nθβ,θ+}

be a class of functions with tails decaying at a faster rate than o(x2n). In [26], it is noted that 𝒱𝒲. Further, let

𝒞c={f𝒞: a compact set 𝕂, such that 𝟏𝕂f=0},

denote the set of compactly supported continuous functions. The following theorem was proved in [27].

Theorem 2 (Nestoridis and Stefanopoulos, 2007, Thm. 3.2).

If g𝒱, then the following statements hold.

  • (a)

    For any f𝒞c, there exists a sequence {ηmg} (ηmg𝒩mg), such that

    limmfηmg1+fηmg=0.
  • (b)

    For any f𝒞0, there exists a sequence {ηmg} (ηmg𝒩mg), such that

    limmfηmg=0.
  • (c)

    For any 1p< and fp, there exists a sequence {ηmg} (ηmg𝒩mg), such that

    limmfηmgp=0.
  • (d)

    For any measurable f, there exists a sequence {ηmg} (ηmg𝒩mg), such that

    limmηmg=f, almost everywhere.
  • (e)

    If ν is a σ-finite Borel measure on n, then for any ν-measurable f, there exists a sequence {ηmg} (ηmg𝒩mg), such that

    limmηmg=f,

    almost everywhere, with respect to ν.

The result was then improved upon, in [26], whereupon the more general space 𝒲 was taken as a replacement for 𝒱, in Theorem 2. Denote the class of bounded continuous functions by 𝒞b=𝒞. The following theorem was proved in [26].

Theorem 3 (Nestoridis et al., 2011, Thm. 3.2).

If g𝒲, then the following statements are true.

  • (a)

    The conclusion of Theorem 2(a) holds, with 𝒞c replaced by 𝒞01.

  • (b)

    The conclusions of Theorem 2(b)–(e) hold.

  • (c)

    For any f𝒞b and compact 𝕂n, there exists a sequence {ηmg}, such that

    limmfηmg(𝕂)=0.

Utilizing the techniques from [27], [1] proved a similar set of results to Theorem 2, under the restriction that f is a non-negative function with support , using g=ϕ (i.e. g has form (1), where n=1) and taking {hmϕ} as the approximating sequence, instead of {ηmg}. That is, the following result is obtained.

Theorem 4 (Bacharoglou, 2010, Cor. 2.5).

If f:+{0}, then the following statements are true.

  • (a)

    For any pdf f𝒞c, there exists a sequence {hmϕ} (hmϕmϕ), such that

    limmfhmϕ1+fhmϕ=0.
  • (b)

    For any f𝒞0, such that f11, there exists a sequence {hmϕ} (hmϕmϕ), such that

    limmfhmϕ=0.
  • (c)

    For any 1<p< and f𝒞p, such that f11, there exists a sequence {hmϕ} (hmϕmϕ), such that

    limmfhmϕp=0.
  • (d)

    For any measurable f, there exists a sequence {hmϕ} (hmϕmϕ), such that

    limmhmϕ=f, almost everywhere.
  • (e)

    For any pdf f𝒞, there exists a sequence {hmϕ} (hmϕmϕ), such that

    limmfhmϕ1=0.

To the best of our knowledge, Theorem 4 is the most complete characterization of the approximating capabilities of the mixture of normal distributions. However, it is restrictive in two ways. First, it does not permit characterization of approximation via the class mg for any g except the normal pdf ϕ. Although ϕ is traditionally the most common choice for g in practice, the modern mixture model literature has seen the use of many more exotic component pdfs, such as the student-t pdf and its skew and modified variants (see, e.g., [29], [13], and [18]). Thus, its use is somewhat limited in the modern context. Furthermore, modern applications tend to call for n>1, further restricting the impact of the result as a theoretical bulwark for finite mixture modeling in practice. A remark in [1] states that the result can generalized to the case where g𝒱 instead of g=ϕ. However, no suggestions were proposed, regarding the generalization of Theorem 4 to the case of n>1.

In this article, we prove a novel set of results that largely generalize Theorem 4. Using techniques inspired by [9] and [4], we are able to obtain a set of results regarding the approximation capability of the class of m-component mixture models mg, when g𝒞0 or g𝒱, and for any n. By definition of 𝒱, the majority of our results extend beyond the proposed possible generalizations of Theorem 4.

The article proceeds as follows. Our main theorem is stated and its seperate parts are proved in the Section 2. Comments and discussion are provided in Section 3. Necessary technical lemmas and results are also included, for reference, in the Appendix.

2 Main result

The remainder of the article is devoted to proving the following theorem.

Theorem 5 (Main result).

If we assume that f and g are pdfs and that g𝒞0, then the following statements are true.

  • (a)

    For any f𝒞0, there exists a sequence {hmg} (hmgmg), such that

    limmfhmg=0.
  • (b)

    For any f𝒞b and compact 𝕂n, there exists a sequence {hmg} (hmgmg), such that

    limmfhmg(𝕂)=0.
  • (c)

    For any 1<p< and fp, there exists a sequence {hmg} (hmgmg), such that

    limmfhmgp=0.
  • (d)

    For any measurable f, there exists a sequence {hmg} (hmgmg), such that

    limmhmg=f, almost everywhere.
  • (e)

    If ν is a σ-finite Borel measure on n, then for any ν-measurable f, there exists a sequence {hmg} (hmgmg), such that

    limmhmg=f,

    almost everywhere, with respect to ν.

If we assume instead that g𝒱, then the following statement is also true.

  • (f)

    For any f𝒞, there exists a sequence {hmg} (hmgmg), such that

    limmfhmg1=0.

2.1 Technical preliminaries

Before we begin to prove the main theorem, we establish some technical results regarding our class of component densities 𝒞0. Let f,g1 and denote the convolution of f and g by fg=gf. Further, we denote the sequence of dilates of g by {gk:gk(x)=kng(kx),k}. The following result is an alternative to Lemma 5 and Corollary 1. Here, we replace a boundedness assumption on the approximand, in the aforementioned theorem by a vanishing at infinity assumption, instead.

Lemma 1.

Let g be a pdf and f𝒞0, such that f>0. Then,

limkgkff=0.
Proof.

It suffices to show that for any ϵ>0, there exists a k(ϵ), such that gkff<ϵ, for all kk(ϵ). By Lemma 6, f𝒞b, and thus f<. By making the substitution z=kx, we obtain for each k

gk(x)dλ=kng(kx)dλ=g(z)dλ=1.

By Corollary 1, we obtain limk𝟏{x:x2>δ}gkdλ=0 and thus we can choose a k(ϵ), such that

𝟏{x:x2>δ}gkdλ<ϵ4f.

Since g is a pdf, we have

|(gkf)(x)f(x)| =|gk(y)[f(xy)f(x)]dλ(y)|
gk(y)|f(xy)f(x)|dλ(y).

By uniform continuity, for any ϵ>0, there exists a δ(ϵ)>0 such that |f(xy)f(x)|<ϵ/2, for any x,yn, such that y2<δ(ϵ) (Lemma 6). Thus, on the one hand, for any δ(ϵ), we can pick a k(ϵ) such that

𝟏{y:y2>δ(ϵ)}gk(y)|f(xy)f(x)|dλ(y)
2f𝟏{y:y2>δ(ϵ)}gkdλ
2f×ϵ4f=ϵ2, (2)

and on the other hand

𝟏{y:y2δ(ϵ)}gk(y)|f(xy)f(x)|dλ(y)
ϵ2𝟏{y:y2δ(ϵ)}gkdλ
ϵ2×1=ϵ2. (3)

The proof is completed by summing (2) and (3). ∎

Lemma 2.

If f𝒞0 is such that f0, and ϵ>0, then there exists a h𝒞c, such that 0hf, and

fh<ϵ
Proof.

Since f𝒞0, there exists a compact 𝕂n such that f(𝕂)<ϵ/2. By Lemma 7, there exists some g𝒞c, such that 0g1 and 𝟏𝕂g=1. Let h=gf, which implies that h0 and 0hf. Furthermore, notice that 𝟏𝕂(fh)=0 and hf, by construction. The proof is completed by observing that

fh =fh(𝕂)
f(𝕂)+h(𝕂)
2f(𝕂)<ϵ.

For any δ>0, uniformly continuous function f, let

w(f,δ)=sup{x,yn:xy2δ}|f(x)f(y)|

denote the modulus of continuity of f. Furthermore, define the diameter of a set 𝕏n by diam(𝕏)=supx,y𝕏xy2 and denote an open ball, centered at xn with radius r>0 by 𝔹(x,r)={yn:xy2<r}.

Notice that the class mg can be parameterized as

mg ={h:h(x)=i=1mciking(kixzi),
 zinki+c𝕊m1i[m]},

where ki=1/σi and zi=μi/σi. The following result is the primary mechanism that permits us to construct finite mixture approximations for convolutions of form gkf. The argument motivated by the approaches taken in Theorem 1 in [4, Ch. 24], [27, Lem. 3.1], and [26, Thm. 3.1].

Lemma 3.

Let f𝒞 and g𝒞0 be pdfs. Furthermore, let 𝕂n be compact and h𝒞c, where 𝟏𝕂h=0 and 0hf. Then for any k, there exists a sequence {hmg}, such that

limmgkhhmg=0.
Proof.

It suffices to show that for any k and ϵ>0, there exists a sufficiently large enough m(ϵ) so that for all mm(ϵ),hmgmg such that

gkhhmg<ϵ. (4)

For any k, we can write

(gkh)(x) =gk(xy)h(y)dλ(y)
=𝟏{y:y𝕂}gk(xy)h(y)dλ(y)
=𝟏{y:y𝕂}kng(kxky)h(y)dλ(y)
=𝟏{z:zk𝕂}g(kxz)h(zk)dλ(z).

Here, k𝕂 is continuous image of a compact set, and hence is compact (cf. [31, Thm. 4.14]). By Lemma 8, for any δ>0, there exists κin (i[m1], m), such that k𝕂i=1m1𝔹(κi,δ/2). Further, if 𝔹iδ=k𝕂𝔹(κi,δ/2), then we have k𝕂=i=1m1𝔹iδ. We can obtain a disjoint covering of k𝕂 by taking 𝔸1δ=𝔹1 and 𝔸iδ=𝔹iδ\j=1i1𝔹jδ (i[m1]) and noting that k𝕂=i=1m1𝔸iδ, by construction (cf. [4, Ch. 24]). Furthermore, each 𝔸iδ is a Borel set and diam(𝔸iδ)δ.

For convenience, let Πmδ={𝔸iδ:i[m1]} denote the disjoint covering, or partition, of k𝕂. We seek to show that there exists an m and Πmδ, such that

gkhi=1mciking(kixzi)<ϵ,

where ki=k,

ci=kn𝟏{z:z𝔸iδ}h(z/k)dλ(z),

and zi𝔸iδ, for i[m1].

Further, zm𝔸m1δ and cm=1i=1m1ci, with km chosen as follows. By Lemma 6, gC< for some positive C. Then, cmkmng(kmxzm)cmkmnC. We may choose km so that kmn=ϵ/(2cmC), so that

cmkmng(kmxzm)ϵ2.

Since 0hf, the sum of ci (i[m1]) satisfies the inequality

i=1m1ci =kni=1m1𝟏{z:z𝔸iδ}h(zk)dλ
=kn𝟏{z:zk𝕂}h(zk)dλ
=𝟏{x:x𝕂}hdλ𝟏{x:x𝕂}fdλfdλ=1.

Thus, 0cm1, and our construction implies that hmgmg, where

hmg(x)=i=1mciking(kixzi)xn.

We can bound the left-hand side of (4) as follows:

gkhhgm
(gkh)(x)i=1m1ciking(kixzi)
 +cmkmng(kmxzm)
(gkh)(x)i=1m1ciking(kixzi)+ϵ2
=𝟏{z:zk𝕂}g(kxz)h(zk)dλ(z)
 i=1m1𝟏{z:z𝔸iδ}g(kxzi)h(zk)dλ(z)+ϵ2
i=1m1𝟏{z:z𝔸iδ}g(kxz)g(kxzi)h(zk)dλ(z)+ϵ2. (5)

Since

kxz(kxzi)2=zzi2diam(𝔸iδ)δ,

we have |g(kxz)g(kxzi)|w(g,δ), for each i[m1]. Since limδ0w(g,δ)=0 (cf. [21, Thm. 4.7.3]), we may choose a δ(ϵ)>0 so that w(g,δ(ϵ))<ϵ/(2kn). We may proceed from (2.1) as follows:

gkhhgm w(g,δ(ϵ))𝟏{z:zk𝕂}h(zk)dλ+ϵ2
=w(g,δ(ϵ))knhdλ+ϵ2
w(g,δ(ϵ))kn+ϵ2
<ϵ2+ϵ2=ϵ. (6)

To conclude the proof, it suffices to choose an appropriate sequence of partitions Πmδ(ϵ),mm(ϵ), for some large but finite m(ϵ), so that (2.1) and (6) hold, which is possible by Lemma 8. ∎

For any r, let 𝔹¯r={xn:x2r} be a closed ball of radius r, centered at the origin.

Lemma 4.

If f1, such that f0, then

limrf𝟏𝔹¯rf1=0.
Proof.

By construction, each element of the sequence {𝟏𝔹¯rf} (r) is measurable, 0𝟏𝔹¯rff, and

limr𝟏𝔹¯rf=f,

point-wise. We obtain our conclusion via the Lesbegue dominated convergence theorem. ∎

2.2 Proof of Theorem 5(a)

We now proceed to prove each of the parts of Theorem 5. To prove Theorem 5(a) it suffices to show that for every ϵ>0, there exists a hmgmg, such that fhmg<ϵ.

Start by applying Lemma 2 to obtain h𝒞c, such that 0hf and fh<ϵ/2. Then, we have

fhmg fh+hhmg
<ϵ2+hhmg. (7)

The goal is to find a hmg, such that hhmg<ϵ/2. Since h𝒞c, we may find a compact 𝕂n such that h(𝕂)=0. Apply Lemma 1 to show the existence of a k(ϵ), such that

hgkh<ϵ4,

for all kk(ϵ). With a fixed k=k(ϵ), apply Lemma 3 to show that there exists a hmgmg, such that

gk(ϵ)hhmg<ϵ4.

By the triangle inequality, we have

hhmg hgk(ϵ)h+gk(ϵ)hhmg
<ϵ4+ϵ4=ϵ2. (8)

The proof is complete by substitution of (8) into (7).

2.3 Proof of Theorem 5(b)

For any ϵ>0 and compact 𝕂n, it suffices to show that there exists a sufficiently large enough m(ϵ) so that for all mm(ϵ),hmgmg, such that fhmg(𝕂)<ϵ.

By Lemma 5, we can find a k(ϵ,𝕂), such that

fgkf(𝕂)<ϵ3, (9)

for every kk(ϵ,𝕂). Since g𝒞0, gC< for some positive C, by Lemma 6. For any k,r, via Young’s convolution inequality:

gkfgk(𝟏𝔹¯rf)knC(𝟏𝔹¯rf)dλ=knCf𝟏𝔹¯rf1. (10)

For fixed k, we may choose r(ϵ,𝕂), using Lemma 4, so that f𝟏𝔹¯rf1ϵ/(3knC) and thus the final term of (10) is bounded from above by ϵ/3 for all rr(ϵ,𝕂). Thus, for k=k(ϵ,𝕂) and, rr(ϵ,𝕂)

gk(ϵ,𝕂)fgk(ϵ,𝕂)(𝟏𝔹¯r(ϵ,𝕂)f)ϵ3. (11)

Using Lemma 3, with approximand 𝟏𝔹¯r(ϵ,𝕂)f, component density g, compact set 𝔹¯r(ϵ,𝕂), h=𝟏𝔹¯r(ϵ,𝕂)f, and with k=k(ϵ,𝕂) fixed, we have the existence of a density hmgmg,mm(ϵ), such that

gk(ϵ,𝕂)(𝟏𝔹¯r(ϵ,𝕂)f)hmgϵ3. (12)

We obtain the desired result by combining (9), (11), and (12), via the triangle inequality.

2.4 Proof of Theorem 5(c)

The technique used to prove Theorem 5(c) is different to those used in the previous sections. Here, we use a result of [9] that generalizes the classic Barron-Jones Hilbert space approximation result (cf. [16] and [2]) to Banach spaces.

To prove Theorem 5(c), it suffices to show that for every ϵ>0, there exists a sufficiently large enough m(ϵ) so that for all mm(ϵ),hmgmg such that fhmgp<ϵ. Begin by applying Corollary 1 to obtain a k(ϵ), such that

fgkfp<ϵ2 (13)

for all kk(ϵ).

For some pdf g and fixed k, let us define the class

𝒢gk={h:h(x)=kng(kxkμ)μn},

write the m-point convex hull of 𝒢gk as

Convm(𝒢gk)={h:h=i=1mcigigi𝒢gkc𝕊m1i[m]},

and call Conv(𝒢gk)=Conv(𝒢gk) the convex hull of 𝒢gk. We further say that Conv¯(𝒢gk) is the closure of Conv(𝒢gk).

Because g is a pdf, g𝒞0𝒞b, and 𝒞b, we observe that g1. Thus, gp, for any 1<p<, by Lemma 9. Since g is a pdf and fp, we have the existence of gkf and the fact that gkfp is finite.

Furthermore, for any ψ𝒢gk, since gp and by definition of 𝒢gk, we have ψpkn/pgp. Thus, we have

ψgkfpψp+gkfpK, (14)

by choosing K=kn/pgp+gkfp>0.

Following [36], we can write the closure of 𝒢gk as

Conv¯(Ggk)={h:h(x)=kng(kxkμ)f(μ)dλ(μ),f is a pdf},

and thus we immediately have gkfConv¯(Ggk). Combined with (14), we can apply Lemma 11 to obtain the conclusion that there exists a function hmgConvm(𝒢gk(ϵ))mg, such that

hmggk(ϵ)fpKCpm11/α,

where α=min{p,2} and Cp is a finite constant. Since p>1, m11/α is strictly increasing, and hence we can choose an m(ϵ), such that for all mm(ϵ),

hmggk(ϵ)fpϵ2. (15)

The proof is then completed by combining (13) and (15) via the triangle inequality.

2.5 Proof of Theorem 5(d) and Theorem 5(e)

By Theorem 5(a), there exists a sequence {hmg} that uniformly converges to f, as m. Thus, by Lemma 12, {hmg} almost uniformly converges to f and also converges almost everywhere, to f, with respect to any measure ν. We prove Theorem 5(d) by setting ν=λ, and we prove Theorem 5(e) by not specifying ν.

2.6 Proof of Theorem 5(f)

It suffices to show that for any ϵ>0, there exists a sufficiently large enough m(ϵ) so that for all mm(ϵ),hmgmg, where g𝒱, such that fhmg1<ϵ. Begin by applying Lemma 4 in order to find a r(ϵ), for any ϵ>0, such that for all rr(ϵ),

f𝟏𝔹¯rf1ϵ24<ϵ2, (16)

where 0𝟏𝔹¯rff, and 𝟏𝔹¯rf𝒞c with compact support 𝔹¯r.

Let 𝕂=𝔹¯r and apply the triangle inequality to obtain

fhmg1 f𝟏𝕂f1+𝟏𝕂fhmg1
ϵ2+𝟏𝕂fhmg1.

Hence we need to show that there exists a function hmgmg, such that

𝟏𝕂fhmg1ϵ2.

Since g𝒱 and gk(x)=kng(kx), by substitution, we have

gk(x)βkθ(k1+x2)n+θ, (17)

where β,θ>0 are independent of k. By Lemma 5 and Corollary 1, we can obtain a k1(ϵ), such that for all kk1(ϵ),

𝟏𝕂fgk(𝟏𝕂f)1ϵ4. (18)

Suppose that γ>1 and let

𝕂k={xn:dist(x,𝕂)kγ},

where

dist(x,𝕏)=inf{xy2:y𝕏}.

By construction, λ(𝕂k)=λ(𝕂)+O(kγ) and thus there exists a k2 such that λ(𝕂k)λ(𝕂)+1, for any kk2.

For any k>k2, we can show that

gk(𝟏𝕂f)hm1g1(𝕂k)<ϵ8. (19)

To do so, firstly, for any xn,

gk(𝟏𝕂f) =𝟏𝕂gk(xy)f(y)dλ(y)
=𝟏k𝕂g(kxz)f(zk)dλ(z).

To obtain a Riemann sum approximation of gk(𝟏𝕂f), we use an argument analogous to that of Lemma 3. That is, we partition k𝕂 into m1 disjoint Borel sets Πm={𝔸1,,𝔸m1}, and we approximate gk(𝟏𝕂f) by a hm1gm1g, where for each i[m1], ki=k, zi𝔸i, and

ci=kn𝟏𝔸if(zk)dλ(z).

Define km+, zmn, and cm=1i=1m1ci, where

cm=fdλ𝟏𝕂fdλ=f𝟏𝕂f1ϵ24 (20)

by (16). Then, by a similar argument to Lemma 3, ci0 for all i[m] and i=1mci=1. Thus, we may define an element hmgmg via the parameters above.

For sufficiently large kk2, we use Lemma 3 to show that

gk(𝟏𝕂f)hm1g(𝕂k)<ϵ8(λ(𝕂)+1),

which implies

gk(𝟏𝕂f)hm1g1(𝕂k) <𝟏𝕂kϵ8(λ(𝕂)+1)dλ
<ϵλ(𝕂k)8(λ(𝕂)+1)<ϵ8, (21)

and thus (19) is proved. Using (19), we write

gk(𝟏𝕂f)hmg1
= gk(𝟏𝕂f)hm1gcmkmng(kmxzm)1
gk(𝟏𝕂f)hm1g1(𝕂k)
+gk(𝟏𝕂f)hm1g1(𝕂k)
+cmkmng(kmxzm)1
ϵ8+cm+gk(𝟏𝕂f)1(𝕂k)+hm1g1(𝕂k),

where cmkmng(kmxzm)1cm since kmng(kmxzm) is a pdf. The aim is now to prove that

gk(𝟏𝕂f)1(𝕂k)<ϵ24 and hm1g1(𝕂k)<ϵ24.

Using polar coordinates and (17), we have

𝟏{x:xy2>kγ}gk(xy)dλ(x)
𝟏{x:xy2>kγ}βkθ(k1+xy2)n+θdλ(x)
=βAnkθ𝟏(kγ,)rn1(k1+r)n+θdλ(r)
βAnkθ𝟏(kγ,)rθ1dλ(r)
=βAnkθ(γ1)/θ,

where An is the surface area of a unit sphere embedded in n. We then have

gk(𝟏𝕂f)1(𝕂k)
=𝟏{y𝕂}𝟏{x𝕂k}f(y)gk(xy)dλ(x)dλ(y)
𝟏𝕂f𝟏{y𝕂}𝟏{x:xy2>kγ}βkθ(k1+xy2)n+θdλ(x)dλ(y)
𝟏𝕂fλ(𝕂)βAnkθ(γ1)/θ,

which implies that we can choose a k3, such that for all kk3,

gk(𝟏𝕂f)1(𝕂k)<ϵ24. (22)

Lastly, we write

hm1g1(𝕂k)
=𝟏𝕂ki=1m1cikng(kxzi)dλ
=i=1m1𝟏𝕂k[kn𝟏𝔸if(zk)dλ(z)]kng(kxzi)dλ(x)
𝟏𝕂fi=1m1knλ(𝔸i)𝟏𝕂kgk(xzik)dλ
𝟏𝕂fi=1m1knλ(𝔸i)βAnkθ(γ1)θ
𝟏𝕂fλ(𝕂)βAnkθ(γ1)/θ,

which implies that we can choose the same k3 as above to obtain the bound

hm1g1(𝕂k)<ϵ24, (23)

for any kk3.

Thus, we obtain the bound 𝟏𝕂fhmg1<ϵ/2, for all kmax{k1,k2,k3}, by combining (18), (19), (20), (21), (22), and (23), via the triangle inequality. The result is proved by combing the bound above, with (16), for an appropriately large r(ϵ).

3 Comments and discussion

3.1 Relationship to Theorem 1

In the proof of Theorem 1, the famous Hilbert space approximation result of [16] and [2] was used to bound the 2 norm between any approximand f2 and a convex combination of bounded functions in 2. This approximation theorem is exactly the p=2 case of the more general theorem of [9], as presented in Lemma 11. Thus, one can view Theorem 5(c) as the p(1,) generalization of Theorem 1.

3.2 The class 𝒲 is a proper subset of the class 𝒞0

Here, we comment on the nature of class 𝒲, which was investigated by [1] and [26]. We recall that [1] conjectured that Theorem 4 generalizes from g=ϕ to g𝒱. In Theorem 5(a)–(e), we assume that g𝒞0. We can demonstrate that g𝒞0 is a strictly weaker condition than g𝒱 or g𝒲.

For example, consider the function in g: such that g(x)=0 if x<0 and

g(x) =i=122ii[(xi+1)2i𝟏{i1x<i1/2}+(xi)2i𝟏{i1/2x<i}] if x0,

and note that

𝟏(1/2,1/2)(2x)2iidλ=12i2+i<1i2.

Since i=1(1/i2)=π2/6, g1. Furthermore, g is continuous since all stationary points of g are continuous. In , g𝒞0 if

limx±g(x)=0.

For x0, we observe that g=0 and thus the left limit is satisfied. On the right, for any 1/ϵ>0, we have x(ϵ)ϵ1/2, so that g(x)<1/ϵ, for all x>x(ϵ), where is the ceiling operator. Therefore, g𝒞0.

Within each interval i1x<i, we observe that g is locally maximized at x=i1/2. The local maximum corresponding to each of these points is 1/i. Thus g𝒲, since

i=11i<ysupx[0,1] |g(x+y)|,

where i=1(1/i)=. Furthermore, g𝒱 since 𝒱𝒲.

3.3 Convergence in measure

Along with the conclusions of Theorem 5(d) and (e), Lemma 12 also implies convergence in measure. That is, if ν is a σ-finite Borel measure on n, then for any ν-measurable f, there exists a sequence {hmg}, such that for any ϵ>0,

limmυ({xn:|f(x)hmg(x)|ϵ})=0.

Appendix A Technical results

Throughout the main text, we utilize a number of established technical results. For the convenience of the reader, we append these results within this Appendix. Sources from which we draw the unproved results are provided at the end of the section.

Lemma 5.

Let {gk} be a sequence of pdfs in 1 and for every δ>0

limk𝟏{x:x2>δ}gkdλ=0.

Then, for all fp and 1p<,

limkgkffp=0.

Furthermore, for all f𝒞b and any compact 𝕂n,

limkgkff(𝕂)=0.

The sequences {gk} from Lemma 5 are often called approximate identities or approximations of the identity. A simple construction of approximate identities is by taking dilations gk(x)=kng(kx), which yields the following corollary.

Corollary 1.

Let g be a pdf. Then the sequence of dilations {gk:gk(x)=kng(kx)}, satisfies the hypothesis of Lemma 5 and hence permits its conclusion.

Lemma 6.

The class 𝒞0 is a subset of 𝒞b. Furthermore, if f𝒞0, then f is uniformly continuous.

Lemma 7 (Urysohn’s Lemma).

If 𝕂n is compact, then there exists some g𝒞c, such that 0g1 and 𝟏𝕂g=1.

Lemma 8.

If 𝕏n is bounded, then for any r>0, 𝕏 can be covered by i=1m𝔹(xi,r) for some finite m, where xin and i[m].

Lemma 9.

If 0<p<q<r, then prq.

Let Γ: be the usual gamma function, defined as Γ(z)=𝟏(0,)xz1exp(x)dλ.

Lemma 10.

If fp and g1, for 1p, then fg exists and we have fgpg1fp.

Lemma 11.

Let 𝒢p, for some 1p<, and let fConv¯(𝒢). For any K>0, such that fαp<K, for all α𝒢, there exists a hmConvm(𝒢), such that

fhmpCpKm11/α,

where α=min{p,2}, and

Cp={1if 1p2,2[πΓ(p+12)]1/pif p>2.
Lemma 12.

In any measure ν, uniform convergence implies almost uniform convergence, and almost uniform convergence implies almost everywhere convergence and convergence in measure, with respect to ν.

Appendix B Sources of results

Lemma 5 is reported as Theorem 9.3.3 in [21] (see also Theorem 2 of [4, Ch. 20]). The proof of Corollary 1 can be taken from that of Theorem 4 of [4, Ch. 20]. Lemma 6 appears in [5], as Proposition 1.4.5. Lemma 7 is taken from Corollary 1.2.9 of [5]. Lemma 8 appears as Theorem 1.2.2 in [5]. Lemma 9 can be found in [12, Prop. 6.10]. Lemma 10 can be found in [21, Thm. 9.3.1]. Lemma 11 appears as Corollary 2.6 in [9]. Lemma 12 can be obtained from the definition of almost uniform convergence, Lemma 7.10, and Theorem 7.11 of [3].

Acknowledgment

HDN is personally funded by Australian Research Council (ARC) grant DE170101134. HDN and GJM are supported by ARC grant DP180101192. FC is supported by Agence Nationale de la Recherche (ANR) grant SMILES ANR-18-CE40-0014 and by Région Normandie grant RIN AStERiCs.

References

  • [1] A G Bacharoglou. Approximation of probability distributions by convex mixtures of Gaussian measures. Proceedings of the American Mathematical Society, 138:2619–2628, 2010. a b c d
  • [2] A R Barron. Universal approximation bound for superpositions of a sigmoidal function. IEEE Transactions on Information Theory, 39:930–945, 1993. a b
  • [3] R G Bartle. The Elements of Integration and Lesbegue Measure. Wiley, New York, 1995. a
  • [4] W Cheney and W Light. A Course in Approximation Theory. Brooks/Cole, Pacific Grove, 2000. a b c d e f g h i
  • [5] J B Conway. A Course in Abstract Analysis. American Mathematical Society, Providence, 2012. a b c
  • [6] A DasGupta. Asymptotic Theory Of Statistics And Probability. Springer, New York, 2008. a b c
  • [7] L De Branges. The Stone-Weierstrass theorem. Proceedings of the American Mathematical Society, 10:822–824, 1959. a
  • [8] L Devroye and G Lugosi. Combinatorial Methods in Density Estimation. Springer, New York, 2000. a
  • [9] M J Donahue, L Gurvits, C Darken, and E Sontag. Rates of convex approximation in non-Hilbert spaces. Constructive Approximation, 13:187–220, 1997. a b c d
  • [10] B S Everitt and D J Hand. Finite Mixture Distributions. Chapman and Hall, London, 1981. a
  • [11] H G Feichtinger. A characterization of Wiener’s algebra on locally compact groups. Archiv der Mathematik, 29:136–140, 1977. a
  • [12] G B Folland. Real Analysis: Modern Techniques and Their Applications. Wiley, New York, 1999. a
  • [13] F Forbes and D Wraith. A new family of multivariate heavy-tailed distributions with variable marginal amounts of tailweights. Statistics and Computing, 24:971–984, 2013. a
  • [14] S Fruwirth-Schnatter. Finite Mixture and Markov Switching Models. Springer, New York, 2006. a
  • [15] S Fruwirth-Schnatter, G Celeux, and C P Robert, editors. Handbook of Mixture Analysis. CRC Press, Boca Raton, 2019. a
  • [16] L K Jones. A simple lemma on greedy approximation in Hilbert space and convergene rates for projection pursuit regression and neural network training. Annals of Statistics, 20:608–613, 1992. a b
  • [17] S Kullback and R A Leibler. On information and sufficiency. Annals of Mathematical Statistics, 22:79–86, 1951. a
  • [18] S X Lee and G J McLachlan. Finite mixtures of canonical fundamental skew t-distributions: the unification of the trstricted and unrestricted skew t-mixture models. Statistics and Computing, 26:573–589, 2016. a
  • [19] J Q Li and A R Barron. Mixture density estimation. In S A Solla, T K Leen, and K R Mueller, editors, Advances in Neural Information Processing Systems, volume 12, Cambridge, 1999. MIT Press. a
  • [20] B G Lindsay. Mixture models: theory, geometry and applications. In NSF-CBMS Regional Conference Series in Probability and Statistics, 1995. a
  • [21] B Makarov and A Podkorytov. Real Analysis: Measures, Integrals and Applications. Springer, London, 2013. a b c
  • [22] G J McLachlan and K E Basford. Mixture Models: Inference and Applications to Clustering. Marcel Dekker, New York, 1988. a
  • [23] G J McLachlan and D Peel. Finite Mixture Models. Wiley, New York, 2000. a
  • [24] Kerrie L Mengersen, Christian Robert, and Mike Titterington. Mixtures: Estimation and Applications. Wiley, New York, 2011. a
  • [25] V Nestoridis and C Papadimitropoulos. Abstract theory of universal series and an application to Dirichlet series. Comptes Rendus Academy of Science Paris Series I, 341:530–543, 2005. a
  • [26] V Nestoridis, S Schmutzhard, and V Stefanopoulos. Universal series induced by approximate identities and some relevant applications. Journal of Approximation Theory, 163:1783–1797, 2011. a b c d e f
  • [27] V Nestoridis and V Stefanopoulos. Universial series and approximate identities. Technical report, Department of Mathematics and Statistics, University of Cyprus, 2007. a b c d
  • [28] H D Nguyen and G J McLachlan. On approximations via convolution-defined mixture models. Communications in Statistics - Theory and Methods, 2019. In press. a b c
  • [29] D Peel and G J McLachlan. Robust mixture modelling using the t distribution. Statistics and Computing, 10:339–348, 2000. a
  • [30] A Rakhlin, D Panchenko, and S Mukherjee. Risk bounds for mixture density estimation. ESAIM: Probability and Statistics, 9:220–229, 2005. a
  • [31] W Rudin. Principles of Mathematical Analysis. McGraw-Hill, Singapore, 1976. a
  • [32] I W Sandberg. Gaussian radial basis functions and inner product sspace. Circuits, Systems and Signal Processing, 20:635–642, 2001. a b c
  • [33] P Schlattmann. Medical Applications of Finite Mixture Models. Springer, Berlin, 2009. a
  • [34] M H Stone. The generalized Weierstrass approximation theorem. Mathematical Magazine, 21:237–254, 1948. a
  • [35] D M Titterington, A F M Smith, and U E Makov. Statistical Analysis of Finite Mixture Distributions. Wiley, New York, 1985. a b
  • [36] S van de Geer. Asymptotic theory for maximum likelihood in nonparametric mixture models. Computational Statistics and Data Analysis, 41:453–464, 2003. a
  • [37] J L Walker and M Ben-Akiva. Advances in discrete choice: mixture models. In Andre De Palma, Robin Lindsey, Emile Quinet, and Roger Vickerman, editors, A Handbook of Transport Economics, pages 160–187. Edward Edgar, Cheltenham, 2011. a
  • [38] Y Xu, W A Light, and E W Cheney. Constructive methods of approximation by ridge functions and radial functions. Numerical Algorithms, 4:205–223, 1993. a b c d
  • [39] G Yona. Introduction to Computational Proteomics. CRC Press, Boca Raton, 2011. a
  • [40] A J Zeevi and R Meir. Density estimation through convex combinations of densities: approximation and estimation bounds. Neural Computation, 10:99–109, 1997. a

Cite this paper

Please cite the published version. Venue: Cogent Mathematics & Statistics, Open-access journal article (Vol. 7, 2020). DOI: 10.1080/25742558.2020.1750861. Official record: Cogent Math. Stat..

BibTeX
@article{nguyen2020approximation,
  title     = {Approximation by finite mixtures of continuous density functions that vanish at infinity},
  author    = {{TrungTin Nguyen} and {Hien D. Nguyen} and {Faicel Chamroukhi} and {Geoffrey J. McLachlan}},
  journal   = {Cogent Mathematics \& Statistics},
  volume    = {7}, number = {1}, pages = {1750861},
  year      = {2020}, publisher = {Taylor \& Francis},
  doi       = {10.1080/25742558.2020.1750861},
}