Abstract and Applied Analysis

Volume 2014 (2014), Article ID 635917, 14 pages

http://dx.doi.org/10.1155/2014/635917

Research Article

Nonparametric Regression with Subfractional Brownian Motion via Malliavin Calculus

Yuquan Cang, Junfeng Liu, and Yan Zhang

School of Mathematics and Statistics, Nanjing Audit University, 86 West Yushan Road, Pukou, Nanjing 211815, China

Received 27 November 2013; Accepted 13 December 2013; Published 6 February 2014

Academic Editor: Litan Yan

Copyright © 2014 Yuquan Cang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We study the asymptotic behavior of the sequence as tends to infinity, where and are two independent subfractional Brownian motions with indices and , respectively. is a kernel function and the bandwidth parameter satisfies some hypotheses in terms of and . Its limiting distribution is a mixed normal law involving the local time of the sub-fractional Brownian motion . We mainly use the techniques of Malliavin calculus with respect to sub-fractional Brownian motion.

1. Introduction

The asymptotic theory of nonlinear regression plays a central role in econometrics, underlying models as diverse as simultaneous equations systems and discrete choice. Two examples in econometrics are GMM estimation and nonlinear cointegration. GMM limit theory was originally developed for ergodic and strictly stationary time series (Hansen [1]) for which all measurable functions are stationary and ergodic, so that applications of strong laws and central limit theory are straightforward. Although some attempts have been made to extend the theory to models with deterministically trending data (e.g., Andrewa and McDermott [2] and Wooldridge [3]), traditional CLT approaches have still been used and no significant progress has been made. Nonlinear cointegrating models also seem important in a range of applications (e.g., Granger [4]) and models with nonlinear attractor sets have been popular in economics for many years.

The motivation of this paper comes from the econometric theory. Consider a nonlinear structural model of cointegration

where is a stationary equilibrium error process, is a jointly dependent nonstationary regressor, and is an unknown function to be estimated with the observed data . Let be a nonnegative real function and set , where the bandwidth parameter . The conventional kernel estimate of in model (1) is given by

The limit behavior of has recently been investigated in Karlsen et al. [5] in the situation where is a recurrent Markov chain. The main theorem in Karlsen et al. ([5], Theorem 3.1) relies on the asymptotic theory developed in Karlsen et al. [5] involving the conditions on the invariant measure associated with a recurrent Markov chain. These conditions are not always easy to check in practice and do not include some cases of econometric interest such as fractional processes. Wang and Phillips [6] provided an alternative approach to nonparametric cointegration in developing the asymptotics. In particular, instead of the recurrent Markov chain in Karlsen et al. [5]), they worked with partial sum representations of the type , where is a Gaussian process or a general linear process, and obtained the limit behavior for kernel functions of this process. This specification corresponds to the conventional formulation of unit root and cointegration models, and the limit theory has links to traditional nonparametric asymptotics for stationary models even though rates of convergence are different. This approach also allows them to work with cases where the regressor is a non-stationary long memory time series. An important assumption in the main part of the above references is the fact that the error process is a martingale difference sequence and there is no contemporaneous correlation between and .

The estimation error in the kernel estimator (2) has the usual decomposition

The second term of (3) affects bias, and, at least when this is of smaller order, it is the first term that determines the asymptotic distribution. The asymptotic behavior of the estimator is usually related to the behavior of the sequence

The limit in distribution as of the sequence has been widely studied in the literature in various situations. Inspired by Bourguin and Tudor [7], we will consider the following situation: we assume that the regressor is a subfractional Brownian motion with Hurst index and the error is , where is a sub-fractional Brownian motion with which is independent of . In this case, our error process is no semimartingale property. We will set with . A supplementary assumption on will be imposed later in terms of the Hurst parameters and . The sequence can be now written as

Our purpose is to give an approach based on stochastic calculus for this asymptotic theory. Recently, the stochastic integration with respect to the sub-fractional Brownian motion has been widely studied (see Yan et al. [8, 9] and references therein). Various types of stochastic integrals, based on Malliavin calculus, has been introduced and change of variables formulas have been derived. We will use all these different techniques in our work. The general idea is as follows. Suppose that . We will first observe that the asymptotic behavior of the sequence will be given by the sum

This is easy to understand since the conditional distribution of given is given by , where is a standard normal random variable. The double sum can be decomposed into two parts: a “diagonal” part and a “nondiagonal” part. We will restrict ourselves to the situation where the diagonal part is dominant (in a sense that will be defined later) with respect to the non-diagonal part. This will imply a certain assumption on the bandwidth parameter in terms of and . We will therefore need to study the asymptotic behavior of

In the case this is actually the bracket of which is a martingale; this motivates our choice of notation. We will assume that the kernel is the standard Gaussian kernel

This choice is motivated by the fact that can be decomposed into an orthogonal sum of multiple Wiener-Itô integrals and the Malliavin calculus can be used to treat the convergence of (7). Its limit in distribution will be after normalization of the local time of the sub-fractional Brownian motion denoted by , where is positive constant. Consequently, we will find that the (renormalized) sequence converges in law to a mixed normal random variable , where is a Brownian motion independent of and is a positive constant.

We also mention that, although the error process does not appear in the limit of (5), it governs the behavior of this sequence. Indeed, the parameter is involved in the renormalization of (5) and the stochastic calculus with respect to is crucial in the proof of our main results.

We have organized our paper as follows. Section 2 contains the notations, definitions, and results from the stochastic calculus that will be needed throughout our paper. In Section 3 we will find the renormalization order of the sequence (5), while Section 4 contains the result of the convergence of the “bracket” (7). In Section 5 we will prove the limit theorem in distribution for .

2. Preliminaries and Notations

In this section we describe the elements from stochastic calculus that we will need in the paper. Let be a standard subfractional Brownian motion (subfBm for short) with parameter . It is well known that this process is a centered Gaussian process with the following covariance function:

for all . For , coincides with the standard Brownian motion . is neither a semimartingale nor a Markov process unless , so many of the powerful techniques from stochastic analysis are not available when dealing with . The subfBm has properties analogous to those of fBm (self-similarity, long-range dependence, and Hölder paths) and satisfies the following estimates:

But its increments are not stationary; more works for sub-fractional Brownian motion can be found in Bojdecki et al. [10–12], Liu et al. [13, 14], Tudor [15], and Yan et al. [8, 9].

Consider a real separable Hilbert space and an isonormal Gaussian process defined on a complete probability space , that is, centered Gaussian family of random variables such that . In this case for sub-fractional Brownian motion, the space is the canonical Hilbert space of the sub-fractional Brownian motion which is defined as the closure of the linear space generated by the indicator functions with respect to the scalar product

Denote by the multiple Itô stochastic integral of a symmetric kernel , with respect to . This is actually an isometry between the Hilbert (symmetric tensor product) equipped with the norm and the th Wiener chaos which is defined as the closed linear span of the random variables , where , and is the Hermite polynomial of degree defined for by

and . The isometry of multiple integrals can be written as follows for positive integers ,

It also holds that

where denotes the symmetrization of defined by

Recall that the Wiener chaos expansion of a square integrable Brownian random variable is given by

where are symmetric functions and .

Let be the Ornstein-Uhlenbeck operator

where is given by (16). If and , we define the Sobolev-Watanabe spaces as the closure of the set of polynomial random variables with respect to the norm

where stands for the identity mapping.

Let us denote by the set of smooth functionals of the form

where and . The Malliavin derivative of a functional as above is given by

The derivative operator is then a closable operator from into . We denote by the closure of with respect to the norm

The divergence integral is the adjoint operator of . That is, we say that a random variable in belongs to the domain of the divergence operator , denoted by , if

for every . In this case is defined by the duality relationship

for any . We have . We will use the notation

to express the Skorohod integral of a process with respect to .

Let be a stochastic process having the chaotic decomposition , where for every . One can prove that if and only if for every , and converges in . In this case,

We denote by the derivative operator, defined on multiple integrals as

This operator is continuous from into . It is known that a random variable belongs to , if and only if its chaotic decomposition satisfies

Set . The Stroock formula that gives the Wiener chaos decomposition of a functional is

For a complete survey of these materials we refer the reader to the book by Nualart [16].

3. Renormalization of the Sequence

We will assume throughout the paper that in (5); then

We compute in this part the -norm of in order to renormalize it. We have

The summand will be called the “diagonal” term while the summand will be called the “non-diagonal” term. We will analyze each of them separately.

Let be the covariance function of the subfBm , and set

In the following we will denote that and , . The following lemma is due to (68) in Yan and Shen [8].

Lemma 1. For all , , and , then one has

where .

Lemma 2. Let ; one has

for .

Proof. Make the change of variable with . One can easily show that the function

is nonpositive by convexity. This completes the proof.

The following lemma, which will be needed in the sequel, studies the properties of the function and the behavior of when goes to infinity.

Lemma 3. There exists a constant , independent of , such that

for all integers , where is positive constant depending only on .

Combining Lemmas 2 and 3, we can easily check the following estimate:

Concerning the term in (30) we have the following.

Proposition 4. As tends to infinity,

Proof. Due to the independence of and , we have that

We denote that for convenience. Then

where is a standard normal random variable. Recall that, if is a standard normal random variable and if ,

Consequently

Then

As tends to infinity, the term behaves as such

where the sign “” means that the left-hand side and the right-hand side have the same limit as tends to infinity. We will use this sign throughout this paper. Hence

Since

so

Thus, according to expression (47), the proof is complete.

We now compute the term .

Proposition 5. Suppose

Then, as tends to infinity,

Proof. Using again the independence of and ,

where

We need to evaluate the expectation . Let be the covariance matrix of given by

We have and

The density of is then

We obtain

Thus

Suppose ; we use Lemma 1 to bound from below; therefore

By virtue of classical inequality with and ,

For the function , according to [17], we can rewrite it as follows:

and thus

So

and consequently

It follows that, under condition (47), converges to zero as tends to infinity.

As a consequence of Propositions 4 and 5 we obtain the following -norm estimate for .

Theorem 6. Suppose condition (47) holds; then, as tends to infinity,

4. Limit Distribution of

Theorem 6 implies that the diagonal part of is dominant in relation to the non-diagonal part, in the sense that this diagonal part is responsible for the renormalization order of which is . As a consequence we need to investigate the limit distribution of

Using the self-similarity of sub-fractional Brownian motion, we have

The limit of sequence (64) is linked to the local time of the sub-fractional Brownian motion . From Geman-Horwitz one can find that the sub-fractional Brownian motion has a local time continuous in which satisfies the occupation formula

for every continuous and bounded function and such that

where denotes the Lebesgue measure and is the Dirac delta function.

Next we will give an important convergence result that will be necessary in proving the main result of this section.

Proposition 7. The following convergence in distribution result holds:

Proof. For fixed , denote by the Gaussian kernel with variance given by . Note that for every

Using (68) we can rewrite the left-hand side of (67) as

We will show that each of , converges to zero (in some sense). Let us first consider the term . We have

It follows from Yan et al. [9] that

in and almost surely, where is the local time of sub-fractional Brownian motion. Then tends to zero as and . For the second term one can express it as follows:

and for every it converges almost surely to zero as tends to infinity using the convergence of Riemann sum convergence. Let us now handle the term given by

We will treat by using the chaos decomposition of the Gaussian kernel applied to random variables in the first Wiener chaos. Recall that for every (which is the canonical Hilbert space associated with the sub-fractional Brownian motion )

where . Using the chaos decomposition we can write almost surely as follows:

where

We will show that converges to zero in as tends to infinity and tends to zero. From (73) we can easily see that the diagonal part of converges to zero. We can also see, from the expression of , that the summands with vanish. Then, by using the orthogonality of multiple stochastic integral, we obtain

We can also write

where

We now claim that, for every fixed ,

In fact, for every , we get

Now for every we have for every and

as tends to infinity because of .

Furthermore, we have that

converges to

as tends to infinity. Since the quantity is finite, it implies (80).

We will now prove

Equalities (80) and (85) will imply the convergence of to zero in . Furthermore one can easily check

By bounding from above the terms and by 1 in , we obtain that

Let us focus on the case where first. For the function defined by (31), due to Yan and Shen [8], one can check . It follows that

Given that, by using Stirling’s formula, the coefficient behaves as , we obtain that the above sum is finite. Thus we obtain the convergence of to zero in for .

Let us now treat the case . We can easily check that the function defined by (31) is increasing on . Since , it holds that . Then

We know that ; by adapting Lemma  2 in [18], we can prove that

with depending only on . As a consequence

The Stirling formula implies again that the above series is convergent.

Theorem 8. Let be defined by (63); then as tends to infinity, one has the convergence in distribution of

where is the local time of sub-fractional Brownian motion .

Proof. Using Proposition 7, it suffices to check that

as tends to infinity. Using the occupation time formula, one can obtain

which converges to as tends to infinity by using the continuity of the local time .

5. Limit Distribution of

In this section, we will investigate the limit in distribution of

Let us consider the Gaussian vector

From the above definition, one can obtain

Theorem 9. Let be given by (95) and assume that

Then one has the convergence in law

where is the local time of and is a Brownian motion independent of .

Proof. We will study the characteristic function of . Let be the imaginary unit and given by

Because of the independence of the two sub-fractional Brownian motions and computing the conditional expectation of given we obtain

It follows that, with ,

Consider the process defined by

and the function defined by . Note that is a bounded variation process (its quadratic variation process is zero). Furthermore

Using the change of variables formula for bounded variation process it follows that

which means that

Therefore

We will focus on the term and show that

From

we obtain

Recall that the following holds:

This can be seen for as an integer at the beginning of the proof and also it can be checked for any . We will use this property to compute the following upper bound for :

This is because all the terms of the double integral are measurable with respect to the filtration generated by . At this point we use (111) to write

Suppose that , which means that and ; consequently

This term is exactly the nondiagonal term of the -norm of and we know that, under condition (98), it converges to zero as tends to infinity. Finally we have

Suppose . It follows that , , which gives us

We get the following result again which is now valid for :

Concerning the term we note that

and the result follows from Theorem 8.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors want to thank the editor and anonymous referees whose valuable remarks and suggestions greatly improved the presentation of this paper. The project is sponsored by NSFC (no. 11226198 and no. 11171062), Innovation Program of Shanghai Municipal Education Commission (12ZZ063), and Priority Academic Program Development of Jiangsu Higher Education Institutions.

References

L. P. Hansen, “Large sample properties of generalized method of moments estimators,” Econometrica, vol. 50, pp. 1029–1054, 1982.

D. W. K. Andrewa and C. J. McDermott, “Nonlinear econometric models with deterministically trending variables,” Reiew of Economic Studies, vol. 62, pp. 343–360, 1995.

J. M. Wooldridge, “Chapter 45 Estimation and inference for dependent processes,” Handbook of Econometrics, vol. 4, pp. 2639–2738, 1994.

C. W. J. Granger, “Nonlinear relationships between extended-memory variables,” Econometrica, vol. 63, no. 2, pp. 265–280, 1995.

H. A. Karlsen, T. Myklebust, and D. Tjøstheim, “Nonparametric estimation in a nonlinear cointegration type model,” Annals of Statistics, vol. 35, no. 1, pp. 252–299, 2007.

Q. Wang and P. C. B. Phillips, “Asymptotic theory for local time density estimation and nonparametric cointegrating regression,” Econometric Theory, vol. 25, no. 3, pp. 710–738, 2009.

S. Bourguin and C. A. Tudor, “Asymptotic theory for fractional regression models via malliavin calculus,” Journal of Theoretical Probability, vol. 25, no. 2, pp. 536–564, 2012.

L. Yan and G. Shen, “On the collision local time of sub-fractional Brownian motions,” Statistics and Probability Letters, vol. 80, no. 5-6, pp. 296–308, 2010.

L. Yan, G. Shen, and K. He, “Itô's formula for a sub-fractional Brownian motion,” Communications on Stochastic Analysis, vol. 5, no. 1, pp. 135–159, 2011.

T. Bojdecki, L. G. Gorostiza, and A. Talarczyk, “Fractional Brownian density process and its self-intersection local time of order k,” Journal of Theoretical Probability, vol. 17, no. 3, pp. 717–739, 2004.

T. Bojdecki, L. G. Gorostiza, and A. Talarczyk, “Limit theorems for occupation time fluctuations of branching systems I: long-range dependence,” Stochastic Processes and their Applications, vol. 116, no. 1, pp. 1–18, 2006.

T. Bojdecki, L. G. Gorostiza, and A. Talarczyk, “Some extensions of fractional brownian motion and sub-fractional browian motion related to particle systems,” Electronic Communications in Probability, vol. 12, pp. 161–172, 2007.

J. Liu and L. Yan, “Remarks on asymptotic behavior of weighted quadratic variation of subfractional Brownian motion,” Journal of the Korean Statistical Society, vol. 41, no. 2, pp. 177–187, 2012.

J. Liu, L. Yan, Z. Peng, and D. Wang, “Remarks on confidence intervals for self-similarity parameter of a subfractional brownian motion,” Abstract and Applied Analysis, vol. 2012, Article ID 804942, 14 pages, 2012.

C. Tudor, “Inner product spaces of integrands associated to subfractional Brownian motion,” Statistics and Probability Letters, vol. 78, no. 14, pp. 2201–2209, 2008.

D. Nualart, Malliavin Calculus and Related Topics, Springer, New York, NY, USA, 2nd edition, 2006.

C. Tudor, “Some properties of the sub-fractional Brownian motion,” Stochastics, vol. 79, no. 5, pp. 431–448, 2007.

M. Eddahbi, R. Lacayo, J. L. Solé, J. Vives, and C. A. Tudor, “Regularity of the local time for the d fractional Brownian motion with N,” Stochastic Analysis and Applications, vol. 23, no. 2, pp. 383–400, 2005.