Jeudi 11 août 2011 4 11 /08 /Août /2011 09:56

 

L'analyse de la mémoire longue et plus précisement les modèles ARFIMA se sont considérablement développés ces dernières années. Certain d'entre vous seront donc intéressés de tester la présence d'integration fractionnaire contre l'alternative d'une racine unitaire. De nombreux tests ont été developpés et je présente ici celui proposé par Lobato & Velasco 2007. Notons que ce test peut également être utilisé pour explorer les relations de cointegration fractionnaire.

The following RATS code is a source file.

Procedure LVFC

*****                                                                   *****
*                                                                           *
*       ******************** Gilles de Truchis *******************          *
*          *********** D1 Rech DEFI | UnivMed (U2) ************             *
*                 ************  2010/2011  ************                     *                  *                                                                           *                         ****                                                                     ****



********
*Compute the Lobato and Velasco (2007) fractional (possibly cointegration) test
*Use the lag option to take into account autocorrelation
********



***************************************** LVFC procedure
procedure lvfc series sbeg send
type series series
type integer sbeg send

    option real d
    option integer lags

    local series zeta
    local series dzeta
    local real sum
    local real dhat
    local series dres
    local series fdres

    local vector[integer] dflags
    local vector coeffs

    dim dflags(fix(send-sbeg)) coeffs(fix(send-sbeg))

    if %defined(d)
    {
        comp dhat = d
    }
    else
    {
        disp "Syntax: @lvfc(d = estimate_of_d) series start end"
        return
    }

    comp sum = 0.0

    diff(diff=1) series sbeg+1 send dres

    ewise dflags(i) = i
    ewise coeffs(i) = %binomial(dhat,i)*(-1)**i
    filter(truncate) series sbeg+1 send fdres
    # dflags
    # coeffs

     set zeta sbeg+1 send = %if(dhat==1.0,0,(fdres - dres)/(1-dhat))

    if %defined(lags).and.lags!=0
    {
        linreg(noprint) fdres sbeg+lags+1 send
        # fdres{1 to lags}

        set dzeta sbeg+lags+1 send = sum=0.0, %do(i,1,lags, sum = sum + %beta(i)*zeta{i}), zeta{0} - sum

        linreg(noprint) dres sbeg+lags+1 send
        # dzeta dres{1 to lags}
    }
    else
    {
        linreg(noprint) dres sbeg+1 send
        # zeta
    }

    comp %lvstat = %tstats(1)
    comp %lvsignif = %ztest(%lvstat)

end

 

Sources: Lobato, I. N., Velasco, C., 2007. Efficient Wald tests for fractional unit roots. Econometrica, Vol. 75, No. 2, 575-589

Par Gilles De Truchis - Publié dans : Econométrie
Ecrire un commentaire - Voir les 0 commentaires
Retour à l'accueil

Un problème d'affichage...?

Certaines pages de ce site contiennent des animations Flash qui nécessitent la dernière version du flash player.

haut de page

Varennes EcoFin

Mon site internet migre vers une nouvelle adresse: Varennes EcoFin. Il s'agit d'un site en flash recentré sur mes travaux de recherches, ma thèse et des resources pour les étudiants (programmes, bases de données etc.). Pour accéder à la version sans flash (HTML5): Varennes EcoFin Sky.

Présentation

Recherche

Music Player

 
Créer un blog gratuit sur over-blog.com - Contact - C.G.U. - Signaler un abus - Articles les plus commentés