Eviews PDF
Eviews PDF
Eviews PDF
Université de Paris–Dauphine
(Mise à jour janv.-12)
LOGICIEL EVIEWS
Ce mode d’emploi du logiciel EVIEWS (http:\\www.eviews.com) est un guide
d’apprentissage à partir du manuel de Régis Bourbonnais, « Économétrie : manuel et
exercices corrigés », Dunod, 8ème éd., 2011. Les données des exercices et les programmes
peuvent être téléchargés sur le site http:\\regisbourbonnais.dauphine.fr puis cliquer sur le
livre Econométrie. Les noms des fichiers correspondent au numéro de chapitre et numéro
d’exercice : par exemple C3EX2 = Exercice 2 du chapitre 3.
Merci de me signaler (sur mon E–mail : Regis.bourbonnais@dauphine.fr) toutes erreurs,
imprécisions ou anomalies.
Menu Général
Environnement Général de la Workfile d’Eviews cette barre est
toujours présente et
EVIEWS accessible
File Edit Object View Proc Quick Opt. Win Hel
p
Ligne de commande : on
peut rentrer directement une
instruction
Chapitre 2 – Exercices 1, 2 et 5
La première étape consiste à créer un espace de travail (WORKFILE) par <file> <new–
workfile>, c’est–à–dire un espace d'accueil pour les séries statistiques (mensuelles,
semestrielles, annuelles, sans dates ...) et les résultats de calcul. Dans la fenêtre proposée, il
convient, soit de choisir la périodicité des données (la date de début et de fin pour les séries
temporelles, sous le forme par exemple : 2009:08 Août 2009, 2009 si les données sont
annuelles), soit d’indiquer le nombre d’observations (séries UNDATED – pour les séries en
coupe instantanée). Dans un WORKFILE ne peut figurer que des séries de même
périodicité.
Créer des séries : <Objects> <new–objects> series puis donner un nom pour la série
(fenêtre à droite). Ensuite sélectionner, dans la workfile, le ou les noms des séries à rentrer, par
exemple REV, (éventuellement faire une sélection de plusieurs noms de séries non adjacentes
en maintenant appuyer sur Ctrl) et double cliquer et cliquer sur Open Group.
Faire un graphique
Dans la barre d’outil de la base de données cliquer sur <View> puis <Graph>.
Pour modifier les options du graphique, double cliquer sur la zone du graphique et faire
Graph option.
Coefficient de corrélation
Pour calculer le coefficient de corrélation entre REV et CONS (évidemment égal à 1)
sélectionner dans la workfile REV et CONS puis double cliquer et choisir <View>
<Correlation>.
View Procs Objects Name Freeze Estimate Forecast Stats Resids (1) t-stat:
Ho: coeff =0
Variable Coefficient St.Error t–stat Probabilty H1: coeff 0
- si Prob < 0.05
C . . . . Rejet Ho
REV . ... ... ... Coeff significatif
- si Prob > 0.05
accept Ho
R2 Mean Dependent Variable coeff non signif
R2 adjusted Std Deviation Dependent Variable (2) F-stat:
Std. Error AIC (Akaike) Ho: coeff = 0
SSE SBC (Schwartz) H1: coeff 0
Si prob < 0.05
Log Likelihood F-stat Rejet Ho, donc
DW Prob of F-stat modèle O.K
Les graphiques de y observé et estimé ainsi que les résidus sont obtenus : <Resid> dans la
fenêtre de « Estimation Output » en haut à droite, ou bien par <View> <Resid>. Pour revenir
aux résultats : <View> <Estimation Output>.
Graphe des résidus : <View> <Actual_fitted_residual>
Forecast Method
1 16 Do Graph
Forecast Evaluation
OK No
Remarque : Chaque type d’objets (équation, graphique, matrice, ...) peut être repéré dans la
workfile par un nom, ce qui permet de les conserver.
Chapitre 3 – Exercice 1
On considère que les données sont sur tableur Excel et on commence par les importer d’Excel
vers EVIEWS.
Une autre méthode, plus simple, consiste à procéder par <copier> <coller>.
On veut, par exemple, calculer les coefficients de régression linéaire, c’est–à–dire la matrice :
a ( X X ) 1 X Y
Dependent Variable: Y
Method: Least Squares
Sample: 1 14
Included observations: 14
Variable Coefficient Std. Error t-Statistic Prob.
C 32.89132 11.66331 2.820068 0.0182
X1 0.801901 0.298436 2.687012 0.0228
X2 -0.381362 0.156581 -2.435564 0.0351
X3 -0.037132 0.052023 -0.713768 0.4917
R-squared 0.702687 Mean dependent var 17.71429
Adjusted R-squared 0.613493 S.D. dependent var 4.177385
S.E. of regression 2.597069 Akaike info criterion 4.981600
Sum squared resid 67.44767 Schwarz criterion 5.164188
Log likelihood -30.87120 F-statistic 7.878181
Durbin-Watson stat 3.186886 Prob(F-statistic) 0.005452
View Procs Objects Name Freeze Estimate Forecast Stats Resids (1) t-stat:
Ho: coeff =0
Variable Coefficient St.Error t–stat Probabilty H1: coeff 0
- si prob < 0.05
C . . . . Rejet Ho
X1 . . . . Coeff significatif
X2 . . . . - si prob > 0,05
X3 . . . . accept Ho
coeff non signif
Il s’agit de l’ordre dans lequel les variables ont été citées lors de l’instruction de régression.
Vérifier par<Views> <Representation> et en double cliquant dans la workfile sur le vecteur C
Remarque : chaque type d’objets – série, équation, matrice, scalaire, ... – est repéré dans la
workfile par un symbole différent.
Omitted Variables: X2 X3
<File> <New–Programme>
Programme :
ls y c x1
scalar sce1 = @ssr/(1–@r2)*@r2
scalar n1 = @ncoef
‘ On donne le nom eq1 à l’équation de régression à trois variables explicatives
equation eq1.ls y c x1 x2 x3
scalar sce = @ssr/(1–@r2)*@r2
scalar n2 = @ncoef
scalar fe = ((sce–sce1)/(n2–n1))/(@ssr/(@regobs–@ncoef))
fe = Fisher empirique
Remarques : une ligne commençant par le symbole ‘ est une ligne de commentaire. Pour
faire apparaître la valeur d’un scalaire il faut « double cliquer » sur le nom de l’objet et sa
valeur est indiquée en bas à droite de la fenêtre.
Programme :
' Test de contrainte linéaire
SCALAR SCR = EQ1.@SSR
SCALAR N = EQ1.@REGOBS – EQ1.@NCOEF
GENR Z = Y – X1
GENR V = X2 + X3
LS Z C V
SCALAR SCR4 = @SSR
SCALAR N4 = @REGOBS – @NCOEF
SCALAR FE = (((SCR4 – SCR)/(N4 – N))/(SCR/N))
Par utilisation des commandes :Estimer le modèle complet, puis <Views CoefficientsTests>
<Wald Coefficient Restriction> et taper dans la fenêtre, C(2) = 1, C(3) = C(4)
Ce qui donne, si l’on considère des séries trimestrielles, sur deux ans :
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
d1 d 2 d 3 d 4
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
Quatre méthodes :
FOR !I = 1 TO 3
FOR !J = !I+1 TO 4
equation eq!a.ls Y C X!I X!J 'equation à deux variables
!a = !a + 1
next
next
FOR !I = 1 TO 2
FOR !J = !I+1 TO 3
Régis Bourbonnais - Logiciel Eviews - Page 11
FOR !K = !J+1 TO 4
equation eq!a.ls Y C X!I X!J X!K 'equation à trois variables
!a = !a + 1
next
next
next
equation eq!a.ls Y C X1 X2 X3 X4
ENDIF
ENDIF
NEXT
– 2) Backward elimination : régression sur le modèle complet puis élimination une par une
des variables explicatives non significatives.
vm = 3 on retient x3
t2 > t1 on retient x2
'calcul de r2yx1.x2x3
ls y c x2 x3
genr e7 = resid
ls x1 c x2 x3
genr e8 = resid
scalar cp4 = @cor(e7,e8)^2
scalar t4 = sqr(cp4/((1–cp4)/(@regobs–2)))
– 4) Stagewise regression
Première étape : Coefficient de corrélation simple le plus élevé (sélection x3), puis
ls y c x3
genr e9 = resid
for !i=1 to 3
scalar cs!i = @cor(x!i,e9)^2
scalar te!i = sqr(cs!i/((1–cs!i)/(@regobs–2)))
next
cs2 > cs1 > sc3, et cs2 significativement 0 (te2 = 2.49 > 2.447)
Régis Bourbonnais - Logiciel Eviews - Page 13
donc on retient x2.
ls y c x2 x3
genr e9 = resid
for !i = 1 to 3
scalar cs!i = @cor(x!i,e9)^2
scalar te!i = sqr(cs!i/((1–cs!i)/(@regobs–2)))
next
'estimation directe de rau partir d’une régression sur les résidus : deuxième méthode'
ls y c x1 x2 x3
genr res = resid
ls res c res(–1)
scalar rau2 = C(2)
genr dy = y–rau2*y(–1) 'On génère les quasi–diff.'
genr dx1 = x1–rau2*x1(–1)
genr dx2 = x2–rau2*x2(–1)
genr dx3 = x3–rau2*x3(–1)
equation eqm2. ls dy c dx1 dx2 dx3
scalar am2 = c(1)/(1–rau2)
– Test de Goldfeld–Quant:
Une variable doit être la cause de l’hétéroscédasticité.
Nombre d’observations doit être important
Omettre C observations centrales (C = partie entière de N/4).
Ici n = 30, alors C = 7 ou 8.
Programme
'calcul du nombre d'observations à retirer (nm)'
smpl 1 30
scalar !nobs = @obs(y)
scalar nm = @ceiling(!nobs/4)
scalar !n1 = @ceiling((!nobs–nm)/2)
smpl 1 !n1
ls y c x
scalar scr1 = @ssr
scalar ddl1 = @regobs–@ncoef
scalar !n2 = !n1 + nm + 1
smpl !n2 !nobs
ls y c x
scalar scr2 = @ssr
scalar ddl2 = @regobs–@ncoef
if scr1>scr2 then scalar fe = (scr1/ddl1)/(scr2/ddl2)
else scalar fe = (scr2/ddl2)/(scr1/ddl1)
endif
scalar test = 0
if @fdist(fe, ddl1, ddl2) < 0.05 then test = 1
endif
Conclusion: si test = 1, alors on rejette H0, hypothèse selon laquelle le modèle est
homoscédastique.
Test de Gleisjer:
Programme
smpl 1 30
equation eq.ls y c x 'Régression sur le modèle de base.
genr resa = abs(resid)
equation eq1.ls resa c x ' Régression de la valeur absolue des résidus sur la variable
genr xra = sqr(x) 'explicative supposée être la cause de l’hétéroscédasticité.
equation eq2.ls resa c xra ' Idem mais avec la racine carrée de la variable explicative.
genr xin = 1 / x
equation eq3.ls resa c xin ' Idem mais avec l’inverse de la variable explicative
scalar proba = 1
for !I = 1 to 3
scalar te =@abs( eq!I.C(2)/sqr(eq!I.@covariance(2,2)))
Régis Bourbonnais - Logiciel Eviews - Page 15
scalar ddl = eq!I.@regobs– eq!I.@ncoef
' On retient la probabilité critique la plus faible et le numéro de l’équation significative
IF @tdist(te, ddl) < proba then
proba = @tdist(te, ddl)
scalar ind = !I
endif
next
' correction de l'héteroscedasticité
genr pon = 1 / sqr(x) 'ne pas oublier que la correction affecte aussi la constante.
genr yp=y*pon
genr xp=x*pon
equation pond1.ls yp xp pon
' Correction de l’hétéroscédasticité par une commande EVIEWS : Régression pondérée
genr pon = 1/sqr(x)
equation pond.ls(w = pon) y c x
' On détruit de la workfile les objets dont on n’a plus besoin
delete eq xra xin resa ddl1 ddl2 ddl scr1 scr2 nm
Le modèle logistique y
y max
t 1 br t
'
' Question 3
'
tend = @trend(1978)
scalar somcr=9999999
for !i=680 to 900 step 10
scalar ymax=!i
genr Y=log(ymax/taux–1)
equation EQ.ls Y c tend
if @ssr<somcr then scalar somcr=@ssr
scalar ymaxmax=ymax
endif
next
genr Y=log(ymaxmax/taux–1)
equation balay.ls Y c tend
scalar b3=exp(C(1))
scalar r3=exp(C(2))
La fonction CES:
‘on affecte des valeurs initiales aux coefficients C(1) = 12 ; C(2) = –1 ; C(3) = 0.5 ; C(4) = –
0.5
param 1 12. 2 –1. 3 0.5 4 –0.5
LS q = c(1)*(c(3)*k^c(4)+(1–c(3))*l^c(4))^c(2)
Programme
'test d’autocorrélation dans modèle autorégressif'
equation eq1.ls po c po(–1) pe
scalar rau = (1–@dw/2)
scalar h = rau*sqr(@regobs/(1–@regobs*sqr(@covariance(2,2)))
if h > 1.96 then 'test bilatéral
scalar test = 1
a
else test = 0
endif 1
rau = 0.4105
Programme
'
' Estimation du modèle à équations simultanées
'
genr tend = @trend(1919)
genr sw = w + wp
' Double Moindres Carrés par programmation
equation eq1.ls p c tend wp t g p(–1) k x(–1)
forcst pa
equation eq2.ls i c pa p(–1) k
'
' D M C par utilisation de la fonction EVIEWS
'
equation eq3.tsls i c p(–1) k p @ tend wp t g p(–1) k x(–1)
En mode interactif les instructions sont – après avoir sélectionnée la série – dans le menu
View <Unit Root Test>. Dans la fenêtre affichée, Eviews propose le choix entre les trois
modèles [1], [2], et [3] et le nombre de décalages dans le cas d’un test DFA (si 0, alors on
effectue le test DF). Le choix du décalage p peut être effectué directement par application des
critères d’informations de Akaike ou Schwarz. L’interprétation du test est la suivante :
t-Statistic Prob.*
Puisque la valeur empirique (ADF Test Statistic –1.8058) est supérieure aux trois valeurs
critiques à 10, 5 et 1%, et que la probabilité critique est de 37,8 % alors on accepte
l’hypothèse H0 d’une série non stationnaire.
Pour un modèle
VAR simple on Unrestricted Vector Auto Regression
On inclue les
choisit l’option variables
“UnRestricted”, Var Specification Series (Groups Included) endogène du
l’autre est pour les UnRestricted Endogenous: système VAR,
modèles ECM Vector Error Correction Y1 Y2 ici, Y1 et Y2
Les résultats complets (estimations des coefficients, écarts types, critères AIC, SC, etc.) sont
fournis. L’option <Views> <Representation> permet d’obtenir la représentation du modèle
VAR. :
y1,t = 0,00676 * y1,t–1 – 0,6125 * y2,t–1 + 17,129
y2,t = –0,1752 * y1,t–1 + 0,2992 * y2,t–1 – 12,862
Le programme permettant de calculer la prévision assortie de son écart type est les suivant :
Expand 78:1 96:4
Var Var1.Ls 1 1 Y1 Y2
var1.makemodel(model1) @ALL F
'Calcul de la prévision
'
smpl 96:1 96:4
Model1.SOLVE
Etape n°1) Modifier <Range> et <Sample> puis dans le modèle VAR estimé <Proc> <Make
model>
Sample: 1 30
Included observations: 28
Test assumption: Linear deterministic trend in the data
Series: Y1 Y2 Y3
Lags interval: 1 to 1
Likelihood 5 Percent 1 Percent Hypothesized
Eigenvalue Ratio Critical Value Critical Value No. of CE(s)
0.603749 37.38888 29.68 35.65 None **
0.229311 11.46904 15.41 20.04 At most 1
0.138550 4.175870 3.76 6.65 At most 2 *
*(**) denotes rejection of the hypothesis at 5%(1%) significance level
L.R. test indicates 1 cointegrating equation(s) at 5% significance level
Pour estimer le modèle de type Logit, nous prenons les options suivantes sous Evieuws
Question 3) Après avoir modifié le <range> et le <sample> (1 - 61) puis renseigné les valeurs
prévues des variables explicatives, le calcul de la prévision est effectué comme pour le
modèle de régression.
Par défaut, ces statistiques font références à la dernière équation de régression utilisée ;
cependant, ces fonctions peuvent être précédées du nom de l’équation de régression :
vente@R2, renvoie la valeur du coefficient de détermination de l’équation de régression
portant le nom vente.
@R2 R2 statistic
@RBAR2 adjusted R2 statistic
@SE standard error of the regression
@SSR sum of squared residuals
@DW Durbin–Watson statistic
@F F–statistic
@LOGL value of the log–likelihood function
@REGOBS number of observations in regression
@AIC Akaike Information Criterion
@SC Schwartz Criterion
@MEANDEP mean of the dependent variable
@SDDEP standard deviation of the dependent variable
@NCOEF total number of estimated coefficients
@COVARIANCE(i,j) covariance of coefficients i and j.
@RESIDCOVA(i,j) covariance of residuals from equation i with those in equation j in a VAR or system
object. @RESIDCOVA must be preceded with a named object, i.e. VAR1.@RESIDCOVA(2,2)
Fonctions diverses
@MOVAV(X,n) n period moving average of X, where n is an integer
@MOVSUM(X,n) n period moving sum of X, where n is an integer
@TREND(d) time trend variable normalized to be zero in period d, where d is a date or
observation number
@SEAS(d) seasonal dummy equal to one when the quarter or month equals d and zero otherwise.
Lois de probabilité
@DNORM(X) standard normal density function of X
@CNORM(X) standard cumulative normal distribution function of X
@TDIST(X, d) Probability that a t–statistic exceeds X with d degrees of freedom
@FDIST(X, n, d) Probability that an F–statistic exceeds X with n numerator degrees of freedom and
d denominator degrees of freedom
@CHISQ(X, d) Probability that a Chi–squared statistic exceeds X with d degrees of freedom
[Component series to save : cocher les deux séries de suffixe _SA et _SF]
Stationnarisation :
Corrélations :
Sélectionner les deux séries à traiter, clic droit > Open as Group
Sélectionner les deux séries à traiter, clic droit > Open as Group
Modèle de régression :
Représentation :
Modélisation VAR
Jarque Bera :
View > Residual Tests > Normality Test > Cholesky of covariance