?***************************************
? Imperfection du marché du travail 
? et decision d'investissement
?***************************************

options memory=8;

?--------------------------------------------------------
? Lecture des donnees
?  Il y a 6 années (de 1987 à 1992) et 1471 individus
?
?  Intitulé des séries
?  invm   : Investissement
?  kapm   : Capital
?  Y      : Production
?  empr   : Endettement
?  resniv : Profit
?  tinvm  : taux d'accumulation
?  ydkm   : rapport production / capital
?  resm   : rapport profit / capital
?  totdkm : rapport endettement /capital
?  dum    : dummy
?
? Lecture des données avec toutes les variables d'un individu
? sur la même ligne
?--------------------------------------------------------
SMPL 1 1471;

READ (FILE='chapitre14.txt')

@id invm87 kapm87 Y87 empr87 resniv87 tinvm87 ydkm87 resm87 totdkm87 dum87

@id invm88 kapm88 Y88 empr88 resniv88 tinvm88 ydkm88 resm88 totdkm88 dum88

@id invm89 kapm89 Y89 empr89 resniv89 tinvm89 ydkm89 resm89 totdkm89 dum89

@id invm90 kapm90 Y90 empr90 resniv90 tinvm90 ydkm90 resm90 totdkm90 dum90

@id invm91 kapm91 Y91 empr91 resniv91 tinvm91 ydkm91 resm91 totdkm91 dum91

@id invm92 kapm92 Y92 empr92 resniv92 tinvm92 ydkm92 resm92 totdkm92 dum92

;

?**********************************************************
?Problème n°2 : Estimation d'un modèle accelérateur profit
?               Méthode des moments GMM
?**********************************************************


LIST inv inv88-inv92;

FRML inv88 tinvm88 - (beta1*tinvm87+ beta3*ydkm87+ beta7*resm87);

FRML inv89 tinvm89 - (beta1*tinvm88+beta3*ydkm88+ +beta7*resm88);

FRML inv90 tinvm90 - (beta1*tinvm89+beta3*ydkm89+ +beta7*resm89);

FRML inv91 tinvm91 - (beta1*tinvm90+beta3*ydkm90+ beta7*resm90);

FRML inv92 tinvm92 - (beta1*tinvm91+ beta3*ydkm91+ beta7*resm91);

FRML dinv89 inv89-inv88;
FRML dinv90 inv90-inv89;
FRML dinv91 inv91-inv90;
FRML dinv92 inv92-inv91;

DOT 89-92;
 EQSUB dinv. inv88-inv92;
ENDDOT;

PARAM beta1  beta3 beta7 ;

READ (nrow=13,ncol=4) selec;
1 1 1 1  
1 1 1 1 
0 1 1 1 
0 0 1 1
0 0 0 1
1 1 1 1 
0 1 1 1 
0 0 1 1
0 0 0 1
1 1 1 1 
0 1 1 1 
0 0 1 1
0 0 0 1
;

GMM(hetero,INST=(C,tinvm87-tinvm90, ydkm87-ydkm90,resm87-resm90),MASK=selec) 
   dinv89 dinv90 dinv91 dinv92;

INPUT proc_lm2test ;

LM2TEST;