PORTMANTEAU TESTS FOR ARMA MODELS WITH INFINITE VARIANCE
|
|
- Bartholomew Lyons
- 6 years ago
- Views:
Transcription
1 doi: /j x PORTMANTEAU TESTS FOR ARMA MODELS WITH INFINITE VARIANCE By J.-W. Lin and A. I. McLeod The University of Western Ontario First Version received Septeber 2006 Abstract. Autoregressive and oving-average (ARMA) odels with stable Paretian errors are soe of the ost studied odels for tie series with infinite variance. Estiation ethods for these odels have been studied by any researchers but the proble of diagnostic checking of fitted odels has not been addressed. In this article, we develop portanteau tests for checking the randoness of a tie series with infinite variance and for ARMA diagnostic checking when the innovations have infinite variance. It is assued that least squares or an asyptotically equivalent estiation ethod, such as Gaussian axiu likelihood, is used. It is also assued that the distribution of the innovations is identically and independently distributed (i.i.d.) stable Paretian. It is seen via siulation that the proposed portanteau tests do not converge well to the corresponding liiting distributions for practical series length so a Monte Carlo test is suggested. Siulation experients show that the proposed Monte Carlo test procedure works effectively. Two illustrative applications to actual data are provided to deonstrate that an incorrect conclusion ay result if the usual portanteau test based on the finite variance assuption is used. Keywords. ARMA odel diagnostic check; portanteau test; residual autocorrelation function; stable paretian distribution; testing for randoness. 1. INTRODUCTION Tie-series odels with stable Paretian errors have been studied by any researchers. Adler et al. (1998) discussed any aspects of how to apply standard Box Jenkins techniques to stable autoregressive oving-average (ARMA) processes. Adler et al. (1998) concluded that, in principle, the standard Box Jenkins techniques do carry over to the stable setting but a great deal of care needs to be exercised. In Section 2 we briefly review the stable Paretian distribution and in Section 3 we develop portanteau tests for whiteness or randoness for an identically and independently distributed (i.i.d.) series. The whiteness test is illustrated with a brief application to daily returns on the S&P 500 stock index. In Section 4 we develop portanteau diagnostic checks for residuals of an AR odel fitted by least squares assuing that the true innovations are i.i.d.-stable Paretiandistributed. This is extended to the ARMA odel using the equality of residuals in AR and ARMA odels. An illustrative exaple shows the differences in inferences that ay result between the finite-variance and infinite-variance portanteau tests /08/ Journal copilation Ó 2008 Blackwell Publishing Ltd., 9600 Garsington Road, Oxford OX4 2DQ, UK and 350 Main Street, Malden, MA 02148, USA.
2 PORTMANTEAU TESTS THE STABLE PARETIAN DISTRIBUTION A stable distribution is usually defined through its characteristic function. A rando variable Z, or Z a (r,b,l), is said to have a stable distribution if its characteristic function has the following for: Eðe itz Þ¼ exp rjtja 1 ib sgnðtþ tan pa 2 þ ilt if a 6¼ 1 exp rjtj 1 þ ib 2 ð1þ p sgnðtþ log jtj þ ilt if a ¼ 1, where i 2 ¼ 1, t the paraeter of the characteristic function, a the index of stability or the characteristic exponent, satisfying 0 < a 2, r > 0 the scale paraeter, b the skewness satisfying 1 b 1, l 2 R 1 is the location paraeter, and ( 1 if t > 0 sgnðtþ ¼ 0 if t ¼ 0 1 if t < PORTMANTEAU TESTS FOR RANDOMNESS OF STABLE PARETIAN TIME SERIES In this section, we study portanteau tests for checking randoness of a sequence of stable Paretian rando variables. We consider the stable analogues of portanteau tests of Box and Pierce (1970) as well as Penˇ a and Rodriguez (2002), denoted by Q BP and ^D, respectively. For this, we require soe iportant properties of saple autocorrelation functions (ACF) and saple partial autocorrelation functions (PACF) of stable Paretian ARMA processes (Brockwell and Davis, 1991, Ch. 13; Saorodnitsky and Taqqu, 1994; Adler et al., 1998) Asyptotic distribution of saple ACF Let fz t : t ¼ 0, ±1, ±2,...g be an (i.i.d.) sequence of stable Paretian rando variables and X t be a strictly stationary process defined by X t ¼ X1 j¼ 1 w j Z t j ; t ¼ 1;...; n; ð2þ where X 1 j¼ 1 jjw j d j < 1; for soe d 2ð0; aþ\½0; 1Š: ð3þ The stable analogue of the autocorrelation function at lag k is defined as Journal copilation Ó 2008 Blackwell Publishing Ltd.
3 602 J.-W. LIN AND A. I. MCLEOD q k ¼ P j w jw jþk P ; j w2 j k ¼ 1; 2;...: ð4þ Equation (4) can be estiated by the saple autocorrelation function as follows: P n k t¼1 r k ¼ P X tx tþk n t¼1 X t 2 ; k ¼ 1; 2;...; ð5þ for a > 0. According to Davis and Resnick (1986), for any positive integer k, the liiting distribution of saple autocorrelation functions is given by 1 n a ð r1 q 1 ;...; r k q k Þ T! ðy 1 ;...; Y k Þ T ; ð6þ where Ô!Õ denotes convergence in distribution and Y h ¼ X1 j¼1 S j q kþj þ q k j 2q j q k ; h ¼ 1;...; k; ð7þ S 0 where S 0, S 1,... are independent stable variables; S 0 is positive with S 0 Z a=2 ðc 2=a a=2 ; 1; 0Þ; and S j are Z a ðca 1=a ; 0; 0Þ, where ð1 aþ=ðcð2 aþ cosðpa=2þþ if a 6¼ 1 C a ¼ 2=p if a ¼ 1. Under the null hypothesis that X t are a sequence of i.i.d. stable Paretian rando variables, we have q 0 ¼ 1 and q k ¼ 0 for k 1, so the liiting distribution of saple ACFs can be further siplified as follows: where W h are given by 1 n a ð r1 ;...; r k Þ T! ðw 1 ;...; W k Þ T ; ð8þ W h ¼ S h ; h ¼ 1;...; k: ð9þ S 0 Note that, for a > 1, we ay also use the ean-corrected saple ACF at lag k, denoted as ~r k, which is given by ~r k ¼ P n k t¼1 ðx t X ÞðX tþk X Þ P n t¼1 ðx t X Þ 2 ; k ¼ 1; 2;...: ð10þ Davis and Resnick (1986) indicated that the liiting distribution of ~r k is the sae as that of r k. Journal copilation Ó 2008 Blackwell Publishing Ltd.
4 PORTMANTEAU TESTS Asyptotic distribution of saple PACF The saple PACF at lag k is defined as the saple estiate of the k-th eleent of the Yule Walker solution to an AR process. The saple PACF ay also be calculated using the Durbin Levison algorith. Let p k be the saple PACF at lag k, and p () ¼ (p 1,..., p ) T. By the Durbin Levison algorith, the vector p () can be expressed as a function of r (), p () ¼ w(r () ), with the k-th eleent given by p k ¼ wðr ðkþ Þ¼ r k r T ðk 1Þ R 1 ðk 1Þ r ðk 1Þ 1 r T ðk 1Þ R 1 ðk 1Þ r ; ð11þ ðk 1Þ where r (p) ¼ (r 1,..., r p ) T is the p 1 vector of saple ACFs, R (p) ¼ (r ji jj ) pp is the p p saple autocorrelation atrix and r ðkþ ¼ðr k;...; r 1 Þ T. Following the proof in Monti (1994), we can derive the asyptotic distribution of saple PACFs. Under the null hypothesis that X t are independent, the ACFs are all zero, and according to Brockwell and Davis (1991, Ch. 13),! n 1=a r h ¼ O p ; h ¼ 1; 2;...: Therefore, R ðkþ ¼ 1 k þ O p! n 1=a ; where 1 k is a k k identity atrix. By eqn (1),! n 2=a p ðþ ¼ r ðþ þ O p : ð12þ Using eqn (8), we have 1 n a ð p1 ;...; p Þ T! ðw 1 ;...; W Þ T : ð13þ 3.3. Asyptotic distributions of Q BP and ^D tests Under the assuption 1 < a < 2, Runde (1997) derived the liiting distribution of Q BP, based on the ean-corrected saple ACFs. His result is given by n 2=aX Q BP ðþ ¼ ~r j 2! W 1 2 þþw 2 ; ð14þ j¼1 where fw k : k ¼ 1,..., g are defined in eqn (9). Note that if 0 < a 1, the liiting distribution of eqn (14) reains the sae if ~r k are replaced by r k. Journal copilation Ó 2008 Blackwell Publishing Ltd.
5 604 J.-W. LIN AND A. I. MCLEOD Consider the ^D test of Penˇ a and Rodriguez (2002) next. In the stable case we ay define the test statistic, n 2=a ^DðÞ ¼ 1 jr ðþj 1= : ð15þ Using the results in Sections 3.1 and 3.2, and following the arguents of Penˇ a and Rodriguez (2002), we ay have the asyptotic distribution of eqn (15) in Theore 1 Theore 1. ^DðÞ in eqn (15) is asyptotically distributed as X i¼1 þ 1 i Wi 2 ; where fw i : i ¼ 1,..., g are as defined in eqn (9). The proof of this theore is given in Appendix A. Reark 1. The liiting distributions of the Q BP and ^D tests can be coputed by aking use of the change variable technique and soe nuerical algoriths of calculating the probability density function of stable rando variables, such as Mittnik et al. (1999). This approach requires, however, intensive nuerical coputations. Reark 2. Another approach to obtaining the asyptotic distributions of the Q BP and ^D tests is to siulate the aforeentioned tests based on their asyptotic distributions. For exaple, ^D is siulated as defined in Theore 1. This approach also requires lengthy coputations but it is uch less intensive coputationally than the approach entioned in Reark 1. This approach will be adopted in the subsequent analysis based on 10 4 siulations. Reark 3. Penˇ a and Rodriguez (2006) consider a slightly different noralization for the ^D statistic and an iproved finite-saple approxiation to its liiting distribution is provided. In the stable case, the new test statistic ay be written, ^D n 2=a ðþ ¼ jr ðþj 1= : ð16þ It ay be shown that the liiting distribution of ^D ðþ is as stated in Theore 1. Moreover, identical results are produced if ^D ðþ is used instead of ^DðÞ in the Monte Carlo test given in Appendix B Siulation experients Based on 250 siulations, the 5, 10, 30, 50, 70, 90, 95, 97.5, 99 (%) epirical quantiles of both Q BP () and ^DðÞ tests with lag ¼ 5 were Journal copilation Ó 2008 Blackwell Publishing Ltd.
6 PORTMANTEAU TESTS 605 calculated and plotted against the corresponding asyptotic distributions as stated in Section 3.3. It is seen in Figures 1 and 2 that the epirical and asyptotic quantiles do not agree very well unless n is very large. The practical applications of both tests would be severely ipeded by the slow convergence. A solution to this proble is to use the Monte Carlo test or paraetric bootstrap. Consider the siulation experients: i.i.d. rando sequence of Z a (1,0,0) with series length n ¼ 250 and a ¼ 1.9,1.7,1.5,1.3,1.1 was siulated. The epirical sizes of both tests were calculated based on N ¼ 10 4 siulations and each Monte Carlo test was siulated based on B ¼ 10 3 siulations see Appendix B for algorith details. The results are tabulated in Table I. It is seen that the epirical sizes of both tests are very close to the 5% noinal level even with n ¼ 250. n = 1000 epirical asyptotic n = 2000 epirical asyptotic D^ D^ n = 5000 epirical asyptotic n = 10,000 epirical asyptotic D^ D^ Figure 1. plot coparing the epirical and asyptotic distributions of ^DðÞ; ¼ 5 defined in eqn (15). Rando sequences of series length n ¼ 1000,2000,5000,10,000 were siulated fro S a (1,0,0), a ¼ 1.2 and epirical quantiles corresponding to 5, 10, 30, 50, 70, 90, 95, 97.5, 99 (%) are shown. Journal copilation Ó 2008 Blackwell Publishing Ltd.
7 606 J.-W. LIN AND A. I. MCLEOD n = 1000 epirical asyptotic n = 2000 epirical asyptotic 0 50 Q BP Q BP 150 n = 5000 epirical asyptotic n = 10,000 epirical asyptotic Q BP Q BP Figure 2. plot coparing the epirical and asyptotic distributions of Q BP (), ¼ 5 defined in eqn (14). Rando sequences of series length n ¼ 1000,2000,5000,10,000 were siulated fro S a (1,0,0),a ¼ 1.2 and epirical quantiles corresponding to 5, 10, 30, 50, 70, 90, 95, 97.5, 99 (%) are shown. TABLE I Epirical Sizes (%) of Monte Carlo Tests for a Noinal 5% Test for Randoness of Stable White Noise Innovations with Index a ^D a ¼ 5 ¼ 10 ¼ 15 ¼ 5 ¼ 10 ¼ The epirical size for each test was calculated based on N ¼ 10 4 siulations. Each Monte Carlo test used B ¼ 10 3 siulations. Series length n ¼ 250 and lags ¼ 5,10,15 were investigated Illustrative exaple The daily S & P 500 stock index fro 2 January, 1999 to 29 Deceber, 2006 was obtained fro Wharton Data Research Services. This results in a series with Journal copilation Ó 2008 Blackwell Publishing Ltd. Q BP
8 PORTMANTEAU TESTS 607 TABLE II P-values Using Monte Carlo Portanteau Tests and Ljung Box v 2 Test, Q LB, for Testing Randoness of Returns of S&P 500 Stock Index Test Statistic Method ¼ 5 ¼ 10 ¼ 20 ¼ 50 ^D Monte Carlo/stable ^D Monte Carlo/noral Q BP Monte Carlo/stable Q BP Monte Carlo/noral Q LB v length n ¼ The returns, log(z tþ1 /z t ), were coputed and tested for randoness. The consistent estiators of McCulloch (1986) were used to estiate a and b for the returns. We obtained ^a M ¼ 1:587 and ^b M ¼ 0:081. It is seen that ^b M is close to zero so the series is not highly skewed. Since ^a M is uch less than 2, the Monte Carlo test using the stable distribution should be used (Appendix B). For coparison, the Monte Carlo tests were also perfored using noral rando variables as well. The Ljung Box tests based on the v 2 approxiation are shown in Table II. The Ljung Box test (Ljung and Box, 1978) rejects at level a when Q LB > v 2 ð1 aþ, where Q LB ðþ ¼nðn þ 2Þ X k¼1 r 2 k ðn kþ ; ð17þ and v 2 ð1 aþ denotes the 1 a quantile of a v2 distribution with degrees of freedo. As expected the Ljung Box test agrees very well with the Monte Carlo test using the Box Pierce statistic Q BP for the Monte Carlo test using noral rando variables. There is a striking difference in the stable/noral Monte Carlo tests when ¼ 50. Under the norality assuption both ^D and Q BP indicate that the randoness assuption is strongly rejected whereas the P-value is only about 5% when the stable assuption is ade. Reark 4. Portanteau tests based on the nonparaetric bootstrap procedure could also be used but it is expected that they would be less powerful since less inforation is used. In addition, the proposed paraetric bootstrap test procedure can be used to test nonlinear tie series driven by stable Paretian innovations. 4. DIAGNOSTIC CHECK FOR MODEL ADEQUACY OF AR(p) MODELS WITH STABLE PARETIAN ERRORS 4.1. Soe asyptotic results In this section, we consider Q BP and ^D tests for diagnostic checks in odel adequacy of AR(p) odels with stable Paretian errors. Journal copilation Ó 2008 Blackwell Publishing Ltd.
9 608 J.-W. LIN AND A. I. MCLEOD Define the general AR(p) process as follows: /ðbþx t ¼ Z t ; ð18þ where fz t : t ¼ 0, ±1, ±2,...g is an i.i.d. sequence of stable Paretian rando variables, B denotes the backward operator, and /(B) ¼ 1 / 1 B / p B p. Let ^/ ðpþ ¼ð^/ 1 ;...; ^/ p Þ T denote the estiates of autoregressive coefficients. The residuals of the fitted odel are given as follows: ^Z t ¼ Z t ð^/ ðpþ Þ¼X t ^/ 1 X t 1 ^/ p X t p ¼ ^/ðbþx t ; ð19þ and the corresponding residual autocorrelation at lag k is given by P ^Z t^z t k ^r k ¼ P : ^Z t 2 Consider the estiators of ^/ ðpþ satisfying! n 1=a ^/ ðpþ ¼ / ðpþ þ O p : Fro Lin and McLeod (2007, Appendix B) the residual autocorrelation at lag k, ^r k, can be approxiated by the first-order Taylor expansion about error ACFs, r k. Specifically, the approxiation is ^r k ¼ r k þ Xp j¼1 ð/ j ^/ j Þw k j þ O p n where w j is the ipulse response coefficient at lag j and P Zt Z t k r k ¼ P Z 2 t! 2=a ; ð20þ is the error autocorrelation at lag k. Equation (20) can also be written in atrix for, to order O p ([n/ log(n)] 2/a ), ^r ðpþ ¼ r ðpþ þ Xð/ ðpþ ^/ ðpþ Þ; ð21þ where w X ¼ : ð22þ w 1 w 2 w p By aking use of eqn (20) or (21) as well as following the proof in Penˇ a and Rodriguez (2002), we ay derive the asyptotic distributions of the Journal copilation Ó 2008 Blackwell Publishing Ltd.
10 PORTMANTEAU TESTS 609 aforeentioned portanteau tests for diagnostic checks in AR(p) odels. This distribution, however, is usually very coplicated and ay not be traceable unless the AR(p) odels of interest are fitted by least squares (LS). For siplicity, therefore, we only consider the case that eqn (18) is estiated using least squares in the subsequent analysis. According to Section 4 in Davis (1996), if the ARMA paraeters, b, are estiated using least squares, we have ½n=Š 1=a ð^b LS bþ converge in distribution, where ^b LS denotes the LS estiates of b. Hence, in ters of our notation, we have ^/ ðpþ / ðpþ ¼ O p ð½n=š 1=a Þ. Then, by Box and Pierce (1970), p f^z t g in eqn (19) satisfy the orthogonality conditions and, to order O p ð1= ffiffi n ½n=Š 1=a Þ, ^r T ðpþ X ¼ 0: ð23þ If we now ultiply eqn (21) on both sizes by then using eqn (23) we have Q ¼ XðX T XÞ 1 X T ; ^r ðpþ ¼ð1 QÞr ðpþ ð24þ approxiately, where 1 is an identity atrix and Q ¼ X(X T X) 1 X T.It was shown by Box and Pierce (1970) that 1 Q is idepotent of rank p. Hence, the asyptotic distribution of the Q BP test is given by n 2=aX ^r k 2 log n! WT ð1 QÞW ; ð25þ 1 where W ¼ (W 1,..., W ) T and fw i : i ¼ 1,..., g are defined in eqn (9). Consider the asyptotic distributions of residual partial autocorrelations next. Let ^p ðþ be the vector of the first residual partial autocorrelations and p () the vector of error partial autocorrelations. The Taylor expansion of wð^r ðþ Þ around r () yields ^p ðþ ¼ p ðþ ðþ ð^r ðþ r ðþ ÞþO p! n 2=a : ð26þ log n By eqns (11) and (12), eqn (26) becoes! n 2=a ^p ðþ ¼ ^r ðþ þ O p : ð27þ log n Consider the Penˇ a Rodriguez test as the for of n 2=a ^D ¼ ð1 j^r ðþj 1= Þ; ð28þ log n Journal copilation Ó 2008 Blackwell Publishing Ltd.
11 610 J.-W. LIN AND A. I. MCLEOD where ^R ðþ ¼ð^r ji jj Þ ; is the residual autocorrelation atrix. By eqn (27) and following the proof in Penˇ a and Rodriguez (2002), the liiting distribution of eqn (28) is W T A W, where A ¼ (1 Q) T W, (1 Q) and W, is an diagonal atrix with (i, i)-th eleent equal to ( i þ 1)/ for i ¼ 1,...,. Reark 5. It is shown in Lin and McLeod (2007, Appendix B.4) that the residuals in a fitted ARMA odel are asyptotically equivalent to those in a particular AR odel. Hence the asyptotic results for the AR ay be extended to the ARMA case Soe size and power calculations As in Section 3.4, the slow convergence of Q BP and ^D tests to their asyptotic distributions is also present at the residual autocorrelations. The first-order autoregressive process X t ¼ 0.5X t 1 þ Z t with Z t Z 1.2 (1,0,0) was siulated and AR(1) odels were fitted to the data. Then the 5, 10, 30, 50, 70, 90, 95, 97.5, 99 (%) epirical quantiles of ^r 1 were copared with the corresponding asyptotic quantiles. It is seen fro Figure 3 that the epirical quantiles of ^r 1 get closer to the asyptotic quantiles as n increases. The slow convergence of residual autocorrelations to their asyptotic distribution ay cause difficulties in using portanteau tests in practice. Therefore, as in Section 3.4, we suggested using the Monte Carlo test to iprove the effectiveness of portanteau tests. We now investigate the effectiveness of Q BP and ^D tests for diagnostic check in fitted AR odels with stable Paretian errors. The epirical sizes of ^D and Q BP tests for a 5% significance test were first calculated via siulation. In this experient, AR(1) odels, X t ¼ / 1 X t 1 þ Z t, were siulated, where Z t Z 1.5 (1,0,0) and / 1 ¼ 0, ±0.1, ±0.3, ±0.5, ±0.7, ±0.9 and AR(1) odels were fitted to the siulated data by the Burg algorith. The epirical size for each test was calculated based on N ¼ 10 4 siulations and each Monte Carlo test used 10 3 siulations. Series length n ¼ 100 and lags ¼ 5,10,20 were investigated. It is seen in Table III that the epirical sizes of both tests are very close to their noinal level. The epirical powers of ^D and Q BP tests as diagnostic tools were also investigated via siulation. Twelve ARMA (2,2) odels of series length n ¼ 100 in Table IV of Penˇ a and Rodriguez (2002) were siulated and AR(1) odels were fitted to the siulated data using the Burg algorith. Both tests with lags ¼ 5,10,20 were calculated using the paraetric bootstrap procedure. The epirical powers were calculated based on N ¼ 10 3 siulations and each Monte Carlo test used 10 3 siulations. It is seen in Table IV that the epirical powers of both tests are reasonably good for ost odels. Soe of the are even better than the powers listed in Penˇ a and Rodriguez (2002). In addition, increasing the series length can also iprove the effectiveness of the proposed test procedure. For exaple, with odel 3 in Table II, if the series length was increased to n ¼ 250, Journal copilation Ó 2008 Blackwell Publishing Ltd.
12 PORTMANTEAU TESTS 611 n = 100 epirical asyptotic n = 500 epirical asyptotic ^ r ^ r 1 n = 1000 epirical asyptotic n = 10,000 epirical asyptotic ^ r ^ r 1 Figure 3. plot ^r 1 and its asyptotic approxiation. The AR(1) process, X t ¼ 0.5X t 1 þ Z t, of series length n ¼ 100,500,10,000 was siulated with stable innovations, fz t gs 1.2 (1,0,0). For each n,10 4 siulations were perfored and AR(1) odels were then fitted to siulated series. The 5, 10, 30, 50, 70, 90, 95, 97.5, 99 (%) epirical quantiles of ^r 1 are plotted and the corresponding asyptotic distribution is shown. the epirical powers of the ^D test at lags ¼ 5,10,20 were increased significantly fro 23.37%, 20.10% and 17.61% to 58.27%, 43.71% and 35.52%, respectively. Siilar iproveent was also found in the Q BP test. Finally, as in Penˇ aand Rodriguez (2002), our siulation experients show that ^D is ore powerful than Q BP as a diagnostic tool. Reark 6. It is well known that the Burg estiate of / 1 is close to the LS estiate. The advantage of using Burg estiate is that it is always in the stationary region and this is needed for the Monte Carlo test Illustrative application Consider the onthly siple returns of CRSP value-weighted index fro January 1926 to Deceber 1997, n ¼ 864 (Tsay, 2002, Ch. 2). After fitting an AR(5), we Journal copilation Ó 2008 Blackwell Publishing Ltd.
13 612 J.-W. LIN AND A. I. MCLEOD TABLE III Epirical Sizes (%) of^d and Q BP for a 5% Significance Test / 1 ¼ 5 ¼ 10 ¼ 20 ^D Q BP ^D Q BP ^D Q BP ^D and Q BP tests for checking odel adequacy of AR(1) odels fitted by the Burg algorith. Both tests were ipleented by the paraetric bootstrap procedure. The epirical size for each test was calculated based on N ¼ 10 4 siulations. Each Monte Carlo test also used B ¼ 10 3 siulations. Series length n ¼ 100 and lags ¼ 5,10,20 were investigated. TABLE IV Epirical Powers (%) of^d and Q BP for a 5% Significance Test ¼ 5 ¼ 10 ¼ 20 Model ^D Q BP ^D Q BP ^D Q BP ^D and Q BP tests for checking odel adequacy of twelve ARMA(2,2) odels in Table 3 of Peňa and Rodriguez (2002) fitted by AR (1) using the Burg algorith. Both tests were ipleented based on the paraetric bootstrap procedure. The epirical power for each test was calculated based on N ¼ 10 4 siulations. Each Monte Carlo test also used B ¼ 10 3 siulations. Series length n ¼ 100 and lags ¼ 5,10,20 were investigated. obtained ^a ¼ 1:635 using the ethod of McCulloch (1986). So the infinite variance hypothesis is plausible for this data. In Table V we copare the P-values obtained for the Monte Carlo tests with ^D and Q BP for both the stable and noral cases with ¼ 10,20,30. The P-value for the Q LB test using the v 2 approxiation is also shown. As expected the Monte Carlo test using Q BP agrees well with the Q LB test using the v 2 approxiation since n is quite large. It is interesting that when ¼ 10 all tests have siilar P-values but when ¼ 20,30 the P-values of Journal copilation Ó 2008 Blackwell Publishing Ltd.
14 PORTMANTEAU TESTS 613 TABLE V An Exaple Using the Monthly Siple Return of CRSP Value-Weighted Index Data fro Tsay (2002) Test statistic Method ¼ 10 ¼ 20 ¼ 30 ^D Monte Carlo/stable Q BP Monte Carlo/stable ^D Monte Carlo/noral Q BP Monte Carlo/noral Q LB v 2 test The data were fitted by an AR(5) odel. The entries in the first two coluns are the P-values, in percent, of ^D S and ^Q S in Section 4 based on the Monte Carlo test; those in the third and fourth coluns are the P-values, in percent, of the Monte Carlo test for ^D N and ^Q N based on the noral distribution and in the last colun the P-value for the standard asyptotic Ljung Box test based on the noral assuption is given. the noral-based tests are uch saller, by a factor of about 10, than the ore correct test based on the stable distribution. In general it appears that using tests based on the assuption of norally distributed innovations ay produce P-values which are too sall when the innovations are generated by a stable distribution. In other words, the test based on the noral distribution would err in the opposite direction as a conservative test and lead to an inflated type I error rate. We ay conclude that in fitting ARMA tie series with stable innovations, it is iportant to use the Monte Carlo diagnostic test presented in this article. 5. CONCLUDING REMARKS McLeod and Li (1983) considered using a portanteau diagnostic test based on the squared residuals. Li (2004, Ch. 5) shows that this test is also useful for testing for the presence of conditional heteroscedasticity in the residuals of fitted ARMA odels. The Monte Carlo test presented in this article is generalized to the squared residuals case in Appendix B and ipleented in our R package (Lin and McLeod, 2007). APPENDIX A. PROOF OF THEOREM 1 First, by decoposing the deterinant of the saple autocorrelation atrix R (), Pena and Rodriguez (2002) showed that jr () j 1/ is a weighted function of the first partial autocorrelations. Specifically, jr ðþ j 1= ¼ Y ð1 p 2 i Þðþ1 iþ= : i¼1 ð29þ Journal copilation Ó 2008 Blackwell Publishing Ltd.
15 614 J.-W. LIN AND A. I. MCLEOD Suppose that under the null hypothesis, ^D is asyptotic distributed as X. By applying the d-ethod to g(x) ¼ log(1 x), it follows that (n/log(n)) 2/a log (jr () j 1/ ) is asyptotically distributed as X. Fro eqn (29), we can have n 2=a logðjr j 1= n Þ¼ 2=aX i¼1 i þ 1 logð1 p 2 i Þ: Next suppose that n 2=a p 2 1 ; p2 2 ;...; T! p2 Y ; ð31þ and apply the ultivariate d-ethod to it follows that gðp 2 1 ; p2 2 ;...; p2 Þ¼ X i¼1 X i¼1 i þ 1 logð1 p 2 i Þ; i þ 1 logð1 p 2 1 i Þ! 1; ;...; 1 Y : Fro the Craer Wold theore, it follows that 1; 1 ;...; 1 n 2=a p 2 1 ;...; n 2=a p 2! 1; 1 ;...; 1 Y ð33þ Under the null hypothesis that X t are a sequence of i.i.d. stable Paretian rando variables, we have q 0 ¼ 1 and q k ¼ 0 for k 1 so the liiting distribution of saple ACFs can be further siplified as follows: 1 n a ð r1 ;...; r k Þ T! ðw 1 ;...; W k Þ T ; ð34þ where W h are given by W h ¼ S h ; h ¼ 1;...; k: ð35þ S 0 Note that, for a > 1, we ay also use the ean-corrected saple autocorrelation function at lag k, denoted as ~r k, which is given by P n k t¼1 ~r k ¼ ðx t X ÞðX tþk X Þ P n t¼1 ðx t X Þ 2 ; k ¼ 1; 2;...: ð36þ Using eqn (8), we have 1 n a ð p1 ;...; p Þ T! ðw 1 ;...; W Þ T : ð37þ it follows that! T ð30þ ð32þ Journal copilation Ó 2008 Blackwell Publishing Ltd.
16 PORTMANTEAU TESTS 615 1; 1 ;...; 1 n 2=a p 2 1 ;...; n 2=a p 2! W1 2 þ 1 W 2 2 þþ1 W 2 ; ð38þ Finally, fro eqns (33) and (38), and fro eqn (31), we have 1; 1 ;...; 1 Y! X ^D! X i¼1 i¼1 þ 1 i Wi 2 : þ 1 i Wi 2 ;! T B. MONTE CARLO TEST FOR RANDOMNESS AND GOODNESS-OF-FIT OF ARMA MODELS The Monte Carlo test procedure for diagnostic checking of AR and ARMA odels with stable Paretian errors can be suarized below. Step 1. Fit an AR odel to data using least squares or the Burg algorith or for ARMA, an approxiate Gaussian axiu likelihood algorith is used. Calculate residuals f^z t g and the portanteau test of interest, say ^D. Step 2. Estiate the paraeters of the stable distribution fro residuals f^z t g in step 1. The quantile estiator given by McCulloch (1986) ay be used or the MLE ipleented in the R package fbasics. Step 3. Select the nuber of Monte Carlo siulations, B. Typically 100 B Step 4. Siulate the fitted odel using the estiated AR or ARMA paraeters in step 1 and ^a in step 2. Obtain ^D after estiating the paraeters in the siulated series. Step 5. Repeat step 4 B ties counting the nuber of ties k that a value of ^D greater than or equal to that in step 1 has been obtained. Step 6. The P-value for the test is (k þ 1)/(B þ 1). Step 7. Reject the null hypothesis if the P-value is saller than a predeterined significance level. This algorith is easily odified to handle the proble of testing a tie series for randoness as well as for Monte Carlo testing when the norality assuption is ade. Siulation of ARMA odels with stable innovations The stable paraeters ay be estiated using the ethod of McCulloch (1986). Independent and identically distributed stable rando variables ay be siulated using the fbasics package. The ARMA(p,q) odel ay be approxiated as a high-order ovingaverage odel, Journal copilation Ó 2008 Blackwell Publishing Ltd.
17 616 J.-W. LIN AND A. I. MCLEOD Y t ¼ : Z t þ w 1 Z t 1 þþw Q Z t Q : ð39þ Initial values Y 1,..., Y r,r ax(p, q) are obtained fro eqn (39) using the fast Fourier transfor to copute the convolution. The reaining portion of the tie series ay be coputed recursively. For speed, these recursions are ipleented in C and interfaced to our R functions. The function siulateara ipleents these ideas in our R package. Nonlinear test In steps 1 and 4 of the algorith in Appendix B, square the residuals, copute the eancorrected saple autocorrelations and use these autocorrelations to copute the test statistic ^D. R package The R package prtest that ipleents all ethods described above is available fro our online suppleent (Lin and McLeod, 2007). Scripts are also given (Lin and McLeod, 2007) for generating all figures and tables which were given in our paper. ACKNOWLEDGEMENTS A. I. McLeod acknowledges with thanks a Discovery Grant Award fro NSERC. The authors would also like to thank the referee for helpful suggestions. NOTE Corresponding author: A.I. McLeod, Departent of Statistical and Actuarial Sciences, The University of Western Ontario, London, Ontario N6A 5B7, Canada. E-ail: aicleod@uwo.ca REFERENCES Adler, R. J., Feldan, R. E. and Gallagher, C. (1998) Analysing stable tie series. In A Practical Guide to Heavy Tails: Statistical Techniques and Applications (eds R. J. Adler, R. E. Feldan and M. S. Taqqu). Boston, MA: Birkhäuser, Box, G. E. P. and Pierce, D. A. (1970) Distribution of residual autocorrelation in autoregressiveintegrated oving average tie series odels. Journal of Aerican Statistical Association 65, Brockwell, P. J. and Davis, R. A. (1991) Tie Series: Theory and Methods. New York: Springer. Journal copilation Ó 2008 Blackwell Publishing Ltd.
18 PORTMANTEAU TESTS 617 Davis, R. A. (1996) Gauss Newton and M-estiation for ARMA processes. Stochastic Processes and their Applications 63, Davis, R. A. and Resnick, S. (1986) Liit theory for the saple covariance and correlation functions of oving averages. The Annals of Statistics 14, Li, W. K. (2004) Diagnostic Checks in Tie Series. New York: Chapan and Hall/CRC. Lin, J.-W. and McLeod, A. I. (2007) Online suppleents to Portanteau Tests for ARMA odels with Infinite Variance (available at Ljung, G. M. and Box, G. E. P. (1978) On a easure of lack of fit in tie series odels. Bioetrika 65, McCulloch, J. H. (1986) Siple consistent estiators of stable distribution paraeters. Counication in Statistics Coputation and Siulation 15, McLeod, A. I. and Li, W. K. (1983) Diagnostic checking ARMA tie series odels using squaredresidual autocorrelations. Journal of Tie Series Analysis 4, Mittnik, S., Rachev, S. T., Doganoglu, T. and Chenyao, D. (1999) Maxiu likelihood estiation of stable Paretian odels. Matheatical and Coputer Modelling 29, Monti, A. C. (1994) A Proposal for residual autocorrelation test in linear odels. Bioetrika 81, Peňa, D. and Rodriguez, J. (2002) A powerful portanteau test of lack of fit for tie series. Journal of Aerican Statistical Association 97, Peňa, D. and Rodriguez, J. (2006) The log of the autocorrelation atrix for testing goodness of fit in tie series. Journal of Statistical Inference and Planning 136, Runde, R. (1997) The asyptotic null distribution of the Box Pierce Q-statistic for rando variable with infinite variance: an application to Geran stock returns. Journal of Econoetrics 78, Saorodnitsky, G. and Taqqu, M. (1994) Stable-Non-Gaussian Rando Processes. New York: Chapan-Hall. Tsay, R. S. (2002) Analysis of Financial Tie Series. New York: Wiley. Journal copilation Ó 2008 Blackwell Publishing Ltd.
Testing the lag length of vector autoregressive models: A power comparison between portmanteau and Lagrange multiplier tests
Working Papers 2017-03 Testing the lag length of vector autoregressive odels: A power coparison between portanteau and Lagrange ultiplier tests Raja Ben Hajria National Engineering School, University of
More informationThe proofs of Theorem 1-3 are along the lines of Wied and Galeano (2013).
A Appendix: Proofs The proofs of Theore 1-3 are along the lines of Wied and Galeano (2013) Proof of Theore 1 Let D[d 1, d 2 ] be the space of càdlàg functions on the interval [d 1, d 2 ] equipped with
More informationKeywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution
Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality
More informationA Powerful Portmanteau Test of Lack of Fit for Time Series
A Powerful Portanteau Test of Lack of Fit for Tie Series Daniel Peña and Julio Rodríguez A new portanteau test for tie series, ore powerful than the tests of Ljung and Box and Monti, is proposed. The test
More informationTesting equality of variances for multiple univariate normal populations
University of Wollongong Research Online Centre for Statistical & Survey Methodology Working Paper Series Faculty of Engineering and Inforation Sciences 0 esting equality of variances for ultiple univariate
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More informationDEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS
ISSN 1440-771X AUSTRALIA DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS An Iproved Method for Bandwidth Selection When Estiating ROC Curves Peter G Hall and Rob J Hyndan Working Paper 11/00 An iproved
More informationMULTISTEP YULE-WALKER ESTIMATION OF AUTOREGRESSIVE MODELS YOU TINGYAN
MULTISTEP YULE-WALKER ESTIMATION OF AUTOREGRESSIVE MODELS YOU TINGYAN NATIONAL UNIVERSITY OF SINGAPORE 2010 MULTISTEP YULE-WALKER ESTIMATION OF AUTOREGRESSIVE MODELS YOU TINGYAN (B.Sc. Nanjing Noral University)
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationEstimation of the Mean of the Exponential Distribution Using Maximum Ranked Set Sampling with Unequal Samples
Open Journal of Statistics, 4, 4, 64-649 Published Online Septeber 4 in SciRes http//wwwscirporg/ournal/os http//ddoiorg/436/os4486 Estiation of the Mean of the Eponential Distribution Using Maiu Ranked
More informationBiostatistics Department Technical Report
Biostatistics Departent Technical Report BST006-00 Estiation of Prevalence by Pool Screening With Equal Sized Pools and a egative Binoial Sapling Model Charles R. Katholi, Ph.D. Eeritus Professor Departent
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationBootstrapping Dependent Data
Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly
More informationTEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES
TEST OF HOMOGENEITY OF PARALLEL SAMPLES FROM LOGNORMAL POPULATIONS WITH UNEQUAL VARIANCES S. E. Ahed, R. J. Tokins and A. I. Volodin Departent of Matheatics and Statistics University of Regina Regina,
More informationA Jackknife Correction to a Test for Cointegration Rank
Econoetrics 205, 3, 355-375; doi:0.3390/econoetrics3020355 OPEN ACCESS econoetrics ISSN 2225-46 www.dpi.co/journal/econoetrics Article A Jackknife Correction to a Test for Cointegration Rank Marcus J.
More informationarxiv: v1 [stat.ot] 7 Jul 2010
Hotelling s test for highly correlated data P. Bubeliny e-ail: bubeliny@karlin.ff.cuni.cz Charles University, Faculty of Matheatics and Physics, KPMS, Sokolovska 83, Prague, Czech Republic, 8675. arxiv:007.094v
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationAn Improved Particle Filter with Applications in Ballistic Target Tracking
Sensors & ransducers Vol. 72 Issue 6 June 204 pp. 96-20 Sensors & ransducers 204 by IFSA Publishing S. L. http://www.sensorsportal.co An Iproved Particle Filter with Applications in Ballistic arget racing
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationStatistical Logic Cell Delay Analysis Using a Current-based Model
Statistical Logic Cell Delay Analysis Using a Current-based Model Hanif Fatei Shahin Nazarian Massoud Pedra Dept. of EE-Systes, University of Southern California, Los Angeles, CA 90089 {fatei, shahin,
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationCS Lecture 13. More Maximum Likelihood
CS 6347 Lecture 13 More Maxiu Likelihood Recap Last tie: Introduction to axiu likelihood estiation MLE for Bayesian networks Optial CPTs correspond to epirical counts Today: MLE for CRFs 2 Maxiu Likelihood
More information1 Proof of learning bounds
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #4 Scribe: Akshay Mittal February 13, 2013 1 Proof of learning bounds For intuition of the following theore, suppose there exists a
More informationA Simple Regression Problem
A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where
More informationA Note on the Applied Use of MDL Approximations
A Note on the Applied Use of MDL Approxiations Daniel J. Navarro Departent of Psychology Ohio State University Abstract An applied proble is discussed in which two nested psychological odels of retention
More informationTail Estimation of the Spectral Density under Fixed-Domain Asymptotics
Tail Estiation of the Spectral Density under Fixed-Doain Asyptotics Wei-Ying Wu, Chae Young Li and Yiin Xiao Wei-Ying Wu, Departent of Statistics & Probability Michigan State University, East Lansing,
More informationSharp Time Data Tradeoffs for Linear Inverse Problems
Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used
More informationIn this chapter, we consider several graph-theoretic and probabilistic models
THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions
More informationMachine Learning Basics: Estimators, Bias and Variance
Machine Learning Basics: Estiators, Bias and Variance Sargur N. srihari@cedar.buffalo.edu This is part of lecture slides on Deep Learning: http://www.cedar.buffalo.edu/~srihari/cse676 1 Topics in Basics
More informationUfuk Demirci* and Feza Kerestecioglu**
1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationAN OPTIMAL SHRINKAGE FACTOR IN PREDICTION OF ORDERED RANDOM EFFECTS
Statistica Sinica 6 016, 1709-178 doi:http://dx.doi.org/10.5705/ss.0014.0034 AN OPTIMAL SHRINKAGE FACTOR IN PREDICTION OF ORDERED RANDOM EFFECTS Nilabja Guha 1, Anindya Roy, Yaakov Malinovsky and Gauri
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationInternational Journal of Pure and Applied Mathematics Volume 37 No , IMPROVED DATA DRIVEN CONTROL CHARTS
International Journal of Pure and Applied Matheatics Volue 37 No. 3 2007, 423-438 IMPROVED DATA DRIVEN CONTROL CHARTS Wille Albers 1, Wilbert C.M. Kallenberg 2 1,2 Departent of Applied Matheatics Faculty
More informationSupport recovery in compressed sensing: An estimation theoretic approach
Support recovery in copressed sensing: An estiation theoretic approach Ain Karbasi, Ali Horati, Soheil Mohajer, Martin Vetterli School of Coputer and Counication Sciences École Polytechnique Fédérale de
More informationSPECTRUM sensing is a core concept of cognitive radio
World Acadey of Science, Engineering and Technology International Journal of Electronics and Counication Engineering Vol:6, o:2, 202 Efficient Detection Using Sequential Probability Ratio Test in Mobile
More informationAn Adaptive UKF Algorithm for the State and Parameter Estimations of a Mobile Robot
Vol. 34, No. 1 ACTA AUTOMATICA SINICA January, 2008 An Adaptive UKF Algorith for the State and Paraeter Estiations of a Mobile Robot SONG Qi 1, 2 HAN Jian-Da 1 Abstract For iproving the estiation accuracy
More informationA method to determine relative stroke detection efficiencies from multiplicity distributions
A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,
More informationBayesian Approach for Fatigue Life Prediction from Field Inspection
Bayesian Approach for Fatigue Life Prediction fro Field Inspection Dawn An and Jooho Choi School of Aerospace & Mechanical Engineering, Korea Aerospace University, Goyang, Seoul, Korea Srira Pattabhiraan
More informationN-Point. DFTs of Two Length-N Real Sequences
Coputation of the DFT of In ost practical applications, sequences of interest are real In such cases, the syetry properties of the DFT given in Table 5. can be exploited to ake the DFT coputations ore
More informationInference in the Presence of Likelihood Monotonicity for Polytomous and Logistic Regression
Advances in Pure Matheatics, 206, 6, 33-34 Published Online April 206 in SciRes. http://www.scirp.org/journal/ap http://dx.doi.org/0.4236/ap.206.65024 Inference in the Presence of Likelihood Monotonicity
More informationThe Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters
journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn
More informationGeneralized Augmentation for Control of the k-familywise Error Rate
International Journal of Statistics in Medical Research, 2012, 1, 113-119 113 Generalized Augentation for Control of the k-failywise Error Rate Alessio Farcoeni* Departent of Public Health and Infectious
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationEXTENDED NEYMAN SMOOTH GOODNESS-OF-FIT TESTS, APPLIED TO COMPETING HEAVY-TAILED DISTRIBUTIONS. J. Huston McCulloch and E. Richard Percy, Jr.
EXTENDED NEYMAN SMOOTH GOODNESS-OF-FIT TESTS, APPLIED TO COMPETING HEAVY-TAILED DISTRIBUTIONS J. Huston McCulloch and E. Richard Percy, Jr.* J. Econoetrics Special Issue Subission: Sept. 4, 200 First Revision:
More informationCompression and Predictive Distributions for Large Alphabet i.i.d and Markov models
2014 IEEE International Syposiu on Inforation Theory Copression and Predictive Distributions for Large Alphabet i.i.d and Markov odels Xiao Yang Departent of Statistics Yale University New Haven, CT, 06511
More informationIN modern society that various systems have become more
Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationE0 370 Statistical Learning Theory Lecture 6 (Aug 30, 2011) Margin Analysis
E0 370 tatistical Learning Theory Lecture 6 (Aug 30, 20) Margin Analysis Lecturer: hivani Agarwal cribe: Narasihan R Introduction In the last few lectures we have seen how to obtain high confidence bounds
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationA proposal for a First-Citation-Speed-Index Link Peer-reviewed author version
A proposal for a First-Citation-Speed-Index Link Peer-reviewed author version Made available by Hasselt University Library in Docuent Server@UHasselt Reference (Published version): EGGHE, Leo; Bornann,
More informationRecursive Algebraic Frisch Scheme: a Particle-Based Approach
Recursive Algebraic Frisch Schee: a Particle-Based Approach Stefano Massaroli Renato Myagusuku Federico Califano Claudio Melchiorri Atsushi Yaashita Hajie Asaa Departent of Precision Engineering, The University
More informationPseudo-marginal Metropolis-Hastings: a simple explanation and (partial) review of theory
Pseudo-arginal Metropolis-Hastings: a siple explanation and (partial) review of theory Chris Sherlock Motivation Iagine a stochastic process V which arises fro soe distribution with density p(v θ ). Iagine
More informationBest Procedures For Sample-Free Item Analysis
Best Procedures For Saple-Free Ite Analysis Benjain D. Wright University of Chicago Graha A. Douglas University of Western Australia Wright s (1969) widely used "unconditional" procedure for Rasch saple-free
More informationA Poisson process reparameterisation for Bayesian inference for extremes
A Poisson process reparaeterisation for Bayesian inference for extrees Paul Sharkey and Jonathan A. Tawn Eail: p.sharkey1@lancs.ac.uk, j.tawn@lancs.ac.uk STOR-i Centre for Doctoral Training, Departent
More informationUvA-DARE (Digital Academic Repository) Recursive unsupervised learning of finite mixture models Zivkovic, Z.; van der Heijden, F.
UvA-DARE (Digital Acadeic Repository) Recursive unsupervised learning of finite ixture odels Zivkovic, Z.; van der Heijden, F. Published in: IEEE Transactions on Pattern Analysis and Machine Intelligence
More informationUsing EM To Estimate A Probablity Density With A Mixture Of Gaussians
Using EM To Estiate A Probablity Density With A Mixture Of Gaussians Aaron A. D Souza adsouza@usc.edu Introduction The proble we are trying to address in this note is siple. Given a set of data points
More informationSymbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm
Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat
More informationStatistical properties of contact maps
PHYSICAL REVIEW E VOLUME 59, NUMBER 1 JANUARY 1999 Statistical properties of contact aps Michele Vendruscolo, 1 Balakrishna Subraanian, 2 Ido Kanter, 3 Eytan Doany, 1 and Joel Lebowitz 2 1 Departent of
More informationOptimal Jackknife for Discrete Time and Continuous Time Unit Root Models
Optial Jackknife for Discrete Tie and Continuous Tie Unit Root Models Ye Chen and Jun Yu Singapore Manageent University January 6, Abstract Maxiu likelihood estiation of the persistence paraeter in the
More informationA MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION
A eshsize boosting algorith in kernel density estiation A MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION C.C. Ishiekwene, S.M. Ogbonwan and J.E. Osewenkhae Departent of Matheatics, University
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationA Comparative Study of Parametric and Nonparametric Regressions
Iranian Econoic Review, Vol.16, No.30, Fall 011 A Coparative Study of Paraetric and Nonparaetric Regressions Shahra Fattahi Received: 010/08/8 Accepted: 011/04/4 Abstract his paper evaluates inflation
More informationTracking using CONDENSATION: Conditional Density Propagation
Tracking using CONDENSATION: Conditional Density Propagation Goal Model-based visual tracking in dense clutter at near video frae rates M. Isard and A. Blake, CONDENSATION Conditional density propagation
More informationResearch Article Integer-Valued Moving Average Models with Structural Changes
Matheatical Probles in Engineering, Article ID 21592, 5 pages http://dx.doi.org/10.1155/2014/21592 Research Article Integer-Valued Moving Average Models with Structural Changes Kaizhi Yu, 1 Hong Zou, 2
More informationSequence Analysis, WS 14/15, D. Huson & R. Neher (this part by D. Huson) February 5,
Sequence Analysis, WS 14/15, D. Huson & R. Neher (this part by D. Huson) February 5, 2015 31 11 Motif Finding Sources for this section: Rouchka, 1997, A Brief Overview of Gibbs Sapling. J. Buhler, M. Topa:
More informationGeneral Properties of Radiation Detectors Supplements
Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent
More informationA remark on a success rate model for DPA and CPA
A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationLost-Sales Problems with Stochastic Lead Times: Convexity Results for Base-Stock Policies
OPERATIONS RESEARCH Vol. 52, No. 5, Septeber October 2004, pp. 795 803 issn 0030-364X eissn 1526-5463 04 5205 0795 infors doi 10.1287/opre.1040.0130 2004 INFORMS TECHNICAL NOTE Lost-Sales Probles with
More informationare equal to zero, where, q = p 1. For each gene j, the pairwise null and alternative hypotheses are,
Page of 8 Suppleentary Materials: A ultiple testing procedure for ulti-diensional pairwise coparisons with application to gene expression studies Anjana Grandhi, Wenge Guo, Shyaal D. Peddada S Notations
More informationOn Constant Power Water-filling
On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives
More informationEmpirical phi-divergence test-statistics for the equality of means of two populations
Epirical phi-divergence test-statistics for the equality of eans of two populations N. Balakrishnan, N. Martín and L. Pardo 3 Departent of Matheatics and Statistics, McMaster University, Hailton, Canada
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationDiagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations
Diagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations Farhat Iqbal Department of Statistics, University of Balochistan Quetta-Pakistan farhatiqb@gmail.com Abstract In this paper
More informationMeasuring Temperature with a Silicon Diode
Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis
More informationSome Proofs: This section provides proofs of some theoretical results in section 3.
Testing Jups via False Discovery Rate Control Yu-Min Yen. Institute of Econoics, Acadeia Sinica, Taipei, Taiwan. E-ail: YMYEN@econ.sinica.edu.tw. SUPPLEMENTARY MATERIALS Suppleentary Materials contain
More informationSoft Computing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis
Soft Coputing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis Beverly Rivera 1,2, Irbis Gallegos 1, and Vladik Kreinovich 2 1 Regional Cyber and Energy Security Center RCES
More informationC na (1) a=l. c = CO + Clm + CZ TWO-STAGE SAMPLE DESIGN WITH SMALL CLUSTERS. 1. Introduction
TWO-STGE SMPLE DESIGN WITH SMLL CLUSTERS Robert G. Clark and David G. Steel School of Matheatics and pplied Statistics, University of Wollongong, NSW 5 ustralia. (robert.clark@abs.gov.au) Key Words: saple
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationKernel-Based Nonparametric Anomaly Detection
Kernel-Based Nonparaetric Anoaly Detection Shaofeng Zou Dept of EECS Syracuse University Eail: szou@syr.edu Yingbin Liang Dept of EECS Syracuse University Eail: yliang6@syr.edu H. Vincent Poor Dept of
More informationEE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng
EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel
More informationTHE KALMAN FILTER: A LOOK BEHIND THE SCENE
HE KALMA FILER: A LOOK BEHID HE SCEE R.E. Deain School of Matheatical and Geospatial Sciences, RMI University eail: rod.deain@rit.edu.au Presented at the Victorian Regional Survey Conference, Mildura,
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 9, September ISSN
International Journal of Scientific & Engineering Research, Volue 4, Issue 9, Septeber-3 44 ISSN 9-558 he unscented Kalan Filter for the Estiation the States of he Boiler-urbin Model Halieh Noorohaadi,
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationMoments of the product and ratio of two correlated chi-square variables
Stat Papers 009 50:581 59 DOI 10.1007/s0036-007-0105-0 REGULAR ARTICLE Moents of the product and ratio of two correlated chi-square variables Anwar H. Joarder Received: June 006 / Revised: 8 October 007
More informationInspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information
Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub
More informationInflation Forecasts: An Empirical Re-examination. Swarna B. Dutt University of West Georgia. Dipak Ghosh Emporia State University
Southwest Business and Econoics Journal/2006-2007 Inflation Forecasts: An Epirical Re-exaination Swarna B. Dutt University of West Georgia Dipak Ghosh Eporia State University Abstract Inflation forecasts
More informationNonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy
Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 22 Oct 1998
arxiv:cond-at/9810285v1 [cond-at.stat-ech] 22 Oct 1998 Statistical Properties of Contact Maps Michele Vendruscolo (1), Balakrishna Subraanian (2), Ido Kanter (3), Eytan Doany (1) and Joel Lebowitz (2)
More informationarxiv: v2 [math.st] 11 Dec 2018
esting for high-diensional network paraeters in auto-regressive odels arxiv:803659v [aths] Dec 08 Lili Zheng and Garvesh Raskutti Abstract High-diensional auto-regressive odels provide a natural way to
More informationCOMPUTER ALGEBRA DERIVATION OF THE BIAS OF LINEAR ESTIMATORS OF AUTOREGRESSIVE MODELS
doi:10.1111/j.1467-9892.2005.00459.x COMPUTER ALGEBRA DERIVATION OF THE BIA OF LINEAR ETIMATOR OF AUTOREGREIVE MODEL By Y. Zhang A. I. McLeod Acadia University The University of Western Ontario First Version
More informationScale Invariant Conditional Dependence Measures
Sashank J. Reddi sjakkar@cs.cu.edu Machine Learning Departent, School of Coputer Science, Carnegie Mellon University Barnabás Póczos bapoczos@cs.cu.edu Machine Learning Departent, School of Coputer Science,
More information