Abstract

The circular kernel density estimator, with the wrapped Cauchy kernel, is derived from the empirical version of Carathéodory function that is used in the literature on orthogonal polynomials on the unit circle. An equivalence between the resulting circular kernel density estimator, to Fourier series density estimator, has also been established. This adds further weight to the considerable role of the wrapped Cauchy distribution in circular statistics.

1. Introduction

Consider an absolutely continuous (with respect to the Lebesgue measure) circular density , ; that is, is -periodic,In the literature on modeling circular data, starting from the classical text of Mardia [1], the standard texts such as Fisher [2], Jammalamadaka and SenGupta [3], and Mardia and Jupp [4] cover parametric models along with many inference problems. More recently various alternatives to these classical parametric models, exhibiting asymmetry and multimodality, have been investigated with respect to their mathematical properties and goodness of fit to some real data; see Abe and Pewsey [5], Jones and Pewsey [6], Kato and Jones [7], Kato and Jones [8], Minh and Farnum [9], and Shimizu and Lida [10].

As elaborated in Nuñez-Antonio et al. [11], circular data possess characteristics such as high skewness or kurtosis and multimodality in many situations, for example, the data on directions of clinical vectorcardiogram (see Downs [12]), wind directions (see Fisher and Lee [13]), and animal orientation (see Oliveira et al. [14]). Such data are not well fitted by standard parametric models, and in such cases, semiparametric or nonparametric modeling may be considered more appropriate.

Fernändez-Durän [15] and Mooney et al. [16] considered modeling of circular data by semiparametric models based on mixture of circular normal and von Mises distributions whereas Hall et al. [17] and Bai et al. [18] considered nonparametric kernel-based density estimation for data on sphere and Fisher [19] and Taylor [20] considered kernel density estimation for circular data. Whereas Fisher [19] adapted the kernels used for linear data to the context of circular data, Taylor [20] used von Mises circular distribution replacing linear kernel used in the classical kernel density estimator that naturally maintains the periodicity in the resulting density estimator. More recently, Di Marzio et al. [21, 22] provided a theoretical basis for circular kernel density estimator by considering the general setting of nonparametric kernel density estimation on a -dimensional torus; the special case of provides circular kernel density estimator that is described below.

Given a random sample from the density (1), the circular kernel density estimator is given by where is a circular kernel where is a concentration parameter and is the mean direction. Note that a circular kernel is usually chosen to be a circular density, unimodal and symmetric around its mean direction which is zero, and it is characterized by a concentration parameter which governs the amount of the smoothing (see, e.g., Di Marzio et al. [21, 22] for details). Classical examples of circular kernels include the von Mises density and the densities of wrapped normal and wrapped Cauchy distributions.

In this note we consider the circular kernel density estimator using the wrapped Cauchy kernel that is given bywhereand we show that it is derived from the empirical version of Carathéodory function, used in the literature on orthogonal polynomials on the unit circle. We also show that this approach leads to Fourier series density estimation; however no truncation of the series is required.

In Section 2, some basic results from the literature on orthogonal polynomials on the unit circle are presented first and then the strategy of estimating is introduced. This in turn produces the nonparametric circular kernel density estimator given in (3). The details are in Section 3. The next section describes the Fourier series estimator of that is shown to be equivalent to the circular kernel density estimator (2) in a limiting sense when wrapped Cauchy kernel is employed.

2. Some Preliminary Results on Orthogonal Polynomials on Unit Circle

Let be the open unit disk, , in the complex plane, and let be a continuous measure defined on the boundary , that is, the circle . The point will be represented by for , and The closure of will be denoted by

Definition 1. A sequence of polynomials defined on are orthogonal with respect to a Borel measure , if they satisfywhere for and it equals 0, otherwise.

The importance of these polynomials is in approximation of bounded functions defined on the unit circle by the representation (see Cantero and Iserles [23])where Here refers to the complex conjugate of . The reader is referred to the excellent book on the subject of orthogonal polynomials on the unit circle by Simon [24]. We will use this result for approximating defined below,

The above function is a “special function,” known as Carathéodory function, that plays an important role in the study of orthogonal polynomials defined on the unit circle (see Simon [24], pp. 25). The orthogonal expansion of with respect to the basis is given by (see and of Simon [25]) whereis the trigonometric moment of the circular distribution with density . The integrand in (8) involves the functionthat is known as the complex Poisson kernel in the theory of complex analysis. Its relation with the wrapped Cauchy kernel lies in the fact thatfor and , where “” denotes the real part. The above function is known as the real Poisson kernel and it is clearly related to the wrapped Cauchy kernel since

A standard result in complex analysis, known as the Poisson representation (see ([24], p. 27)), says that if is analytic in a neighborhood of , with real, then for , This representation leads to the result (see (ii) in Section of Simon [24]) that for Lebesgue a.e. exists and that

3. Smooth Circular Density Estimator Derived from an Estimator of

Recognizing the expression for as the expectation of , its empirical version is given by Thus we may define an estimator of motivated by , the identities (15) and (16), as

Now we show that the above estimator is of the same form as the circular density estimator given in (3). Recognize that where ; then using (12), we have and thereforewhich is of the same form as in (3).

4. A Connection of the Circular Density Estimator to Fourier Series Estimator

The orthogonal series for given in (9) may be directly used to define an estimator of . This method also provides the same estimator as derived in the previous section that is demonstrated as follows. Estimating the coefficients , , by an estimator of is given by Substituting the expression for from (22), we can write which is the same as given in (19).

Equation (9) may be used to derive the Fourier series estimator. The reader may be referred to Efromovich [26] for the details about Fourier series density estimator. Truncating the infinite sum in (9) at some large index , we have an estimator of Here is chosen according to some criteria, for example, to minimize the integrated squared error. The reader may be referred to Efromovich [26]. Thus we have two contrasting situations; in one we have to choose and in the other case we have to choose a suitable concentration parameter . Numerically as well as technically, the second choice, that is, the circular kernel density estimator, may be considered more suitable.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

The author would also like to acknowledge the partial support from NSERC, Canada, through a Discovery Grant to the author.