17.08.2013 Views

THÈSE Estimation, validation et identification des modèles ARMA ...

THÈSE Estimation, validation et identification des modèles ARMA ...

THÈSE Estimation, validation et identification des modèles ARMA ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapitre 4. Multivariate portmanteau test for weak structural V<strong>ARMA</strong> models 133<br />

8. Define the estimator<br />

ˆΦm = 1<br />

n<br />

9. Define the estimators<br />

n<br />

<br />

ê<br />

′<br />

t−1,...,ê ′ ′ ∂<strong>et</strong>(θ0)<br />

t−m ⊗<br />

∂θ ′<br />

<br />

θ0= ˆ .<br />

θn<br />

t=1<br />

ˆΣˆ = Γm ˆ Σγm + ˆ Φm ˆ Σˆ ˆΦ θn<br />

′ m + ˆ Φm ˆ Σˆ + θn,γm ˆ Σ ′ ˆΦ θn,γm ˆ ′ m<br />

ˆΣˆρm =<br />

<br />

Im ⊗( ˆ Se ⊗ ˆ Se) −1<br />

<br />

ˆΣˆ Im ⊗( Γm<br />

ˆ Se ⊗ ˆ Se) −1<br />

<br />

10. Compute the eigenvalues ˆ ξm = ( ˆ ξ1,d2m,..., ˆ ξd2m,d2m) ′ of the matrix<br />

<br />

ˆΣˆΓm<br />

ˆΩm =<br />

Im ⊗ ˆ Σ −1/2<br />

e0 ⊗ ˆ Σ −1/2<br />

e0<br />

Im ⊗ ˆ Σ −1/2<br />

e0 ⊗ ˆ Σ −1/2<br />

<br />

e0 .<br />

11. Compute the portmanteau statistics<br />

Pm = nˆρ ′ <br />

m Im ⊗ ˆR −1<br />

e (0)⊗ ˆ R −1<br />

e (0)<br />

<br />

ˆρm and<br />

˜Pm = n 2<br />

m 1<br />

(n−h) Tr<br />

<br />

ˆΓ ′<br />

e (h) ˆ Γ −1<br />

e (0)ˆ Γe(h) ˆ Γ −1<br />

e (0)<br />

<br />

.<br />

h=1<br />

12. Evaluate the p-values<br />

P<br />

<br />

Zm( ˆ <br />

ξm) > Pm<br />

and P<br />

<br />

Zm( ˆ ξm) > ˜ <br />

Pm , Zm( ˆ d<br />

ξm) =<br />

2 m<br />

i=1<br />

ˆξi,d 2 mZ 2 i,<br />

using the Imhof algorithm (1961). The BP test (resp. the LB test) rejects the<br />

adequacy of the weak V<strong>ARMA</strong>(p,q) model when the first (resp. the second) pvalue<br />

is less than the asymptotic level α.<br />

4.9 Numerical illustrations<br />

In this section, by means of Monte Carlo experiments, we investigate the finite<br />

sample properties of the test introduced in this paper. For illustrative purpose, we only<br />

present the results of the modified and standard versions of the LB test. The results<br />

concerning the BP test are not presented here, because they are very close to those<br />

of the LB test. The numerical illustrations of this section are made with the softwares<br />

R (see http ://cran.r-project.org/) and FORTRAN (to compute the p-values using the<br />

Imohf algorithm, 1961).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!