Research Article | Open Access

Volume 2016 |Article ID 9505794 | 8 pages | https://doi.org/10.1155/2016/9505794

# Empirical Likelihood Inference for First-Order Random Coefficient Integer-Valued Autoregressive Processes

Accepted04 Nov 2015
Published08 Feb 2016

#### Abstract

We apply the empirical likelihood method to estimate the variance of random coefficient in the first-order random coefficient integer-valued autoregressive (RCINAR(1)) processes. The empirical likelihood ratio statistic is derived and some asymptotic theory for it is presented. Furthermore, a simulation study is presented to demonstrate the performance of the proposed method.

#### 1. Introduction

Integer-valued time series data are fairly common in practice. Especially in economics and medicine, many interesting variables are integer-valued. In the last three decades, integer-valued time series have received increasing attention because of their wide applicability in many different areas, and there were many developments in the literature on it. See, for instance, Davis et al.  and MacDonald and Zucchini . For count data, so far there are two main classes of time series models that have been developed in recent years: state-space models and thinning models. For state-space models, we refer to Fukasawa and Basawa . Integer-valued autoregressive (INAR(1)) model was first defined by Steutel and Harn  through the “thinning” operator . Recall the definition of a “thinning” operator : where is an integer-valued random variable and and is an i.i.d. Bernoulli random sequence with that is independent of . Based on the “thinning” operator , the INAR(1) model is defined aswhere is a sequence of i.i.d. nonnegative integer-valued random variables.

Note that the parameter may be random and it may vary with time; Zheng et al.  introduced the following first-order random coefficient integer-valued autoregressive (RCINAR(1)) model:where is an independent identically distributed sequence with cumulative distribution function on with and ; is a sequence of i.i.d. nonnegative integer-valued random variables with and . Moreover, and are independent.

Zheng et al.  further generalized the above model to the -order cases. In recent several years, RCINAR model has been studied by many authors (see references in ). In this paper, we are concerned with estimating the variance of random coefficient in model (3). We propose an empirical log-likelihood ratio statistics for and derive its asymptotic distribution which is standard .

As a nonparametric statistical method, the empirical likelihood method was introduced by Owen . The advantages of the empirical likelihood are now widely recognized. It has sampling properties similar to the bootstrap. Many advantages of the empirical likelihood over the normal approximation-based method have also been shown in the literature. These attractive properties have motivated various authors to extend empirical likelihood methodology to other situations. Now, the empirical likelihood methods have been widely applied to the statistical inference of the time series models (see ).

The remainder of this paper is organized as follows: In Section 2, we introduce the methodology and the main results. Simulation results are given in Section 3. Section 4 provides the proofs of the main results.

Throughout the paper, we use the notations “” and “” to denote convergence in distribution and convergence in probability, respectively. Convergence “almost surely” is written as “a.s.” Furthermore, denotes the transpose matrix of the matrix , and denotes Euclidean norm of the matrix or vector.

#### 2. Methodology and Main Results

In this section, we will first discuss how to apply the empirical likelihood method to estimate the unknown parameter .

Let , and . For simplicity of notation, we write as ; parameters and will be omitted. Then, after simple algebra, we get and , where .

First we consider estimating by using the conditional least-squares method. Based on the sample , the least-squares estimator of can be obtained by minimizing with . Solving the equationfor , we haveLet , where and are given by Zheng et al. . Further let and . Then, the estimating equation of can be written aswhere .

In what follows, we apply Owen’s empirical likelihood method to make inference about . For convenience of writing, let be a probability vector with and ; also, let denote the true parameter value for . The log empirical likelihood ratio evaluated at , a candidate value of , is By using the Lagrange multiplier method, introducing a Lagrange multiplier , we havewhere satisfies

Owen’s empirical log-likelihood ratio statistic has a chi-squared limiting distribution. Similarly, we can prove that will also be asymptotically chi-squared distributed. In order to establish a theory for , we assume that the following assumptions hold:() is a strictly stationary and ergodic process.().

Remark 1. Similar conditions can be found in .

Now we can give the limiting properties of .

Theorem 2. Assume that and hold. If is the true value of , thenwhere is a chi-squared distribution with degree of freedom.

As a consequence of the theorem, confidence regions for the parameter can be constructed by (12). For , an asymptotic confidence region for is given by where is the upper -quantile of the chi-squared distribution with degrees of freedom equal to 1.

#### 3. Simulation Study

In this section, we conduct some simulation studies which show that our proposed methods perform very well.

In the first simulation study, we consider the RCINAR(1) process:where is a sequence of i.i.d. sequence with and ; . We take and and take and . Samples of size and All simulation studies are based on 1000 repetitions. The results of the simulations are presented in Table 1. The nominal confidence level is chosen to be and , and the figures in parentheses are the simulation results at the nominal level of .

 = 50 = 100 = 300 λ = 1 0.1 0.980 (0.960) 0.976 (0.949) 0.979 (0.956) 0.2 0.972 (0.946) 0.984 (0.971) 0.984 (0.966) 0.3 0.977 (0.958) 0.986 (0.959) 0.990 (0.970) 0.4 0.983 (0.943) 0.984 (0.958) 0.984 (0.968) 0.5 0.978 (0.963) 0.986 (0.966) 0.990 (0.977) λ = 2 0.1 0.980 (0.952) 0.972 (0.943) 0.970 (0.975) 0.2 0.977 (0.956) 0.980 (0.980) 0.973 (0.979) 0.3 0.989 (0.961) 0.983 (0.967) 0.969 (0.978) 0.4 0.983 (0.969) 0.983 (0.966) 0.974 (0.981) 0.5 0.973 (0.968) 0.982 (0.965) 0.970 (0.971)

From Table 1, we find that the confidence region obtained by using the empirical likelihood method has high coverage levels for different . The coverage probability has no obvious change for different and . That means that the empirical likelihood method is also robust.

In the second simulation study, we illustrate how our method can be applied to fit a set of data through a practical example. We apply model (3) to fit the number of large- and medium-sized civil Boeing 767 cargo planes over the period 1985–2013 in China. The data in Table 2 are provided by the National Bureau of Statistics of China (http://data.stats.gov.cn/easyquery.htm?cn=C01). The fitting procedure is as follows: Firstly, by using the data over the period 1985–2003, we obtain the estimator of the model parameter. Then, by using this model, we can obtain a fitting sequence over the period 2004–2013. Furthermore, in order to compare with the ordinary autoregressive (AR(1)) model, we also give the fitting results of the AR(1) model. Table 3 reports the fitting results. In Table 3, Number is the true value and RCINAR(1) and AR(1) are the fitting results obtained by the RCINAR(1) model and AR(1) model, respectively. For the simulation results of AR(1) model, we take the rounded integer values of the simulation results. From the simulation results, we can find that the RCINAR(1) model has more plausible fitting results than the AR(1) model.

 Year 1985 1986 1987 1988 1989 1990 Number 2 2 4 5 6 6 Year 1991 1992 1993 1994 1995 1996 Number 6 10 12 16 17 17 Year 1997 1998 1999 2000 2001 2002 Number 17 15 16 16 17 18 Year 2003 2004 2005 2006 2007 2008 Number 22 27 27 29 22 22 Year 2009 2010 2011 2012 2013 Number 19 18 15 13 11
 Year 2004 2005 2006 2007 2008 2009 Number 27 27 29 22 22 19 RCINAR(1) 23 22 22 21 19 18 AR(1) 22 22 23 23 23 25 Year 2010 2011 2012 2013 Number 18 15 13 11 RCINAR(1) 18 18 16 16 AR(1) 25 27 27 29

#### 4. Proofs of the Main Results

Lemma 3. Assume that and hold. Thenwhere and .

Proof. Note thatFirst, we consider . After simple algebra calculation, we haveNext, we consider . By the mean value theorem, we have where lies between and and lies between and . Therefore,Below, we prove that . For , note that By Theorem 3.1 in Zheng et al. , we know thatMoreover, by the ergodic theorem, we haveFurther note thatwhich, combined with (22) and (23), implies thatSimilarly, we can prove thatNext, we prove thatNote thatBy the ergodic theorem, we haveBy (21), we haveTherefore, by (21), we haveSimilarly, we can prove thatUsing this, together with (25), (26), and (27), we can proveFinally, we prove thatFor this, we first prove thatBy the Cramer-Wold device, it suffices to show that, for all ,Let and Then is a zero-mean, square integrable martingale array. By making use of a martingale central limit theorem , we can prove (36). Further, by (23), we know that (34) holds. Therefore, by (17), (33), and (34), we can prove Lemma 3.

Lemma 4. Assume that and hold. Then

Proof. Note thatBy (23), in order to prove Lemma 4, we have only to show thatNote thatBy the ergodic theorem, we know thatSimilar to the proof of (33), we can further prove thatThis, in conjunction with (41), yields (39). So we complete the proof of Lemma 4.

Lemma 5. Assume that and hold. Then

Proof. To prove (43), we only need to prove thatLet For , definewhere denotes the largest integer not greater than . For each , (37) implies that as . Moreover, note that For given , choose so that Therefore, for each , if , then we have So, for any ,which implies (44). So we prove (43).

Proof of Theorem 2. First, we prove thatWrite , where and Observe that This implies thatFurther, by Lemma 4, we know that By Lemma 3, we have Thus by (51) and Lemma 5, we have which implies (49).
By (49) and Lemma 5, we can prove thatExpanding (11), we haveBy (55) and Lemmas 3, 4, and 5, we know that the final term in (56) is bounded byThis, together with (41), yieldsBy the Taylor expansion, we have Below, we prove that there exists a finite number , such thatThe Taylor expansion of around yields where, as . Thus, there exists , such that for any Moreover, by (55), we have Let . Note that if , then for ,which implies that where
Moreover, by (10) and (58), we haveThis, together with Lemmas 3 and 4, implies Theorem 2.

#### Conflict of Interests

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

#### Acknowledgments

The authors acknowledge the financial supports by National Natural Science Foundation of China (nos. 11571138, 11271155, 11001105, 11071126, 10926156, and 11071269), Specialized Research Fund for the Doctoral Program of Higher Education (no. 20110061110003), Program for New Century Excellent Talents in University (NCET-08-237), Scientific Research Fund of Jilin University (no. 201100011), and Jilin Province Natural Science Foundation (nos. 20130101066JC, 20130522102JH, and 20101596).

1. R. A. Davis, T. M. Dunsmuir, and Y. Wang, “Modeling time series of count data,” in Asymptotics, Nonparametrics and Time Series, pp. 63–114, Marcel-Dekker, New York, NY, USA, 1999. View at: Google Scholar
2. I. L. MacDonald and W. Z. Zucchini, Hidden Markov and Other Models for Discrete-Valued Time Series, Chapman & Hall, London, UK, 1997.
3. T. Fukasawa and I. V. Basawa, “Estimation for a class of generalized state-space time series models,” Statistics & Probability Letters, vol. 60, no. 4, pp. 459–473, 2002. View at: Publisher Site | Google Scholar
4. F. W. Steutel and K. V. Harn, “Discrete analogues of self-decomposability and stability,” The Annals of Probability, vol. 7, no. 5, pp. 893–899, 1979. View at: Publisher Site | Google Scholar | MathSciNet
5. H. Zheng, I. V. Basawa, and S. Datta, “First-order random coefficient integer-valued autoregressive processes,” Journal of Statistical Planning and Inference, vol. 137, no. 1, pp. 212–229, 2007. View at: Publisher Site | Google Scholar | MathSciNet
6. H. Zheng, I. V. Basawa, and S. Datta, “Inference for $p$th-order random coefficient integer-valued autoregressive processes,” Journal of Time Series Analysis, vol. 27, no. 3, pp. 411–440, 2006. View at: Publisher Site | Google Scholar | MathSciNet
7. J. Kang and S. Lee, “Parameter change test for random coefficient integer-valued autoregressive processes with application to polio data analysis,” Journal of Time Series Analysis, vol. 30, no. 2, pp. 239–258, 2009.
8. H. Zhang, D. Wang, and F. Zhu, “The empirical likelihood for first-order random coefficient integer-valued autoregressive processes,” Communications in Statistics. Theory and Methods, vol. 40, no. 3, pp. 492–509, 2011.
9. H. Zhang, D. Wang, and F. Zhu, “Empirical likelihood inference for random coefficient INAR(p) process,” Journal of Time Series Analysis, vol. 32, no. 3, pp. 195–203, 2011. View at: Publisher Site | Google Scholar | MathSciNet
10. A. Roitershtein and Z. Zhong, “On random coefficient INAR(1) processes,” Science China Mathematics, vol. 56, no. 1, pp. 177–200, 2013. View at: Publisher Site | Google Scholar | MathSciNet
11. A. B. Owen, “Empirical likelihood ratio confidence intervals for a single functional,” Biometrika, vol. 75, no. 2, pp. 237–249, 1988. View at: Publisher Site | Google Scholar | MathSciNet
12. A. B. Owen, “Empirical likelihood ratio confidence region,” The Annals of Statistics, vol. 18, no. 1, pp. 90–120, 1990. View at: Publisher Site | Google Scholar | MathSciNet
13. A. B. Owen, “Empirical likelihood for linear models,” The Annals of Statistics, vol. 19, no. 4, pp. 1725–1747, 1991. View at: Publisher Site | Google Scholar | MathSciNet
14. C.-S. Chuang and N. H. Chan, “Empirical likelihood for autoregressive models, with applications to unstable time series,” Statistica Sinica, vol. 12, no. 2, pp. 387–407, 2002.
15. S. X. Chen, W. Härdle, and M. Li, “An empirical likelihood goodness-of-fit test for time series,” Journal of the Royal Statistical Society—Series B: Statistical Methodology, vol. 65, no. 3, pp. 663–678, 2003. View at: Publisher Site | Google Scholar | MathSciNet
16. N. H. Chan and S. Ling, “Empirical likelihood for GARCH models,” Econometric Theory, vol. 22, no. 3, pp. 403–428, 2006. View at: Publisher Site | Google Scholar | MathSciNet
17. S. X. Chen and J. T. Gao, “An adaptive empirical likelihood test for parametric time series regression models,” Journal of Econometrics, vol. 141, no. 2, pp. 950–972, 2007. View at: Publisher Site | Google Scholar | MathSciNet
18. P. Guggenberger and R. J. Smith, “Generalized empirical likelihood tests in time series models with potential identification failure,” Journal of Econometrics, vol. 142, no. 1, pp. 134–161, 2008. View at: Publisher Site | Google Scholar | MathSciNet
19. J. Li, W. Liang, S. He, and X. Wu, “Empirical likelihood for the smoothed LAD estimator in infinite variance autoregressive models,” Statistics & Probability Letters, vol. 80, no. 17-18, pp. 1420–1430, 2010. View at: Publisher Site | Google Scholar
20. Z.-W. Zhao and D.-H. Wang, “Empirical likelihood for an autoregressive model with explanatory variables,” Communications in Statistics. Theory and Methods, vol. 40, no. 3, pp. 559–570, 2011. View at: Publisher Site | Google Scholar | MathSciNet
21. Z.-W. Zhao, D.-H. Wang, and C.-X. Peng, “Coefficient constancy test in generalized random coefficient autoregressive model,” Applied Mathematics and Computation, vol. 219, no. 20, pp. 10283–10292, 2013. View at: Publisher Site | Google Scholar | MathSciNet
22. P. Hall and C. C. Heyde, Martingale Limit Theory and Its Application, Academic Press, New York, NY, USA, 1980. View at: MathSciNet

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions. 