**question 1**
end 1
CALENDAR(PANELOBS=28) 1970
ALLOCATE 0 34//1997:1


OPEN DATA convergence.xls
DATA(FORMAT=xls,ORG=obs) /

PRINT /
*impression de toutes les séries*


disp '******************'
disp '** question 2 **'
disp '******************'
SET va = VAB80/EFTOT
SET fbcf = FBCF80/VAB80

disp '******************'
disp '** question 3 **'
disp '******************'
COMPUTE sec = 34
COMPUTE per = 28
DEC VECTOR[SERIES] D(sec)
DO i=1,sec
SET D(i) = %INDIV(T)==i
END DO I
*creation des variables indicatrices pour chaque secteur*

SET logva = log(va)
SET logfbcf = log(fbcf)
SET dlogva = logva-logva{1}

**équation 1**
DISP '***estimation secteur par secteur*****'
DO i=1,sec
LINREG(SMPL=d(i)) dlogva
# logva{1} logfbcf constant
END DO I

**équation 2**
DISP '***estimation avec un effet fixe sectoriel*****'
PREGRESS(EFFECTS=indiv,METHOD=fixed) dlogva
# logva{1} logfbcf

**équation 3**
DISP '***estimation avec un effet fixe temporel*****'
PREGRESS(EFFECTS=time,METHOD=fixed) dlogva
# logva{1} logfbcf

**équation 4**
DISP '***estimation avec un effet sectoriel aléatoire*****'
PREGRESS(EFFECTS=indiv,METHOD=random) dlogva
# logva{1} logfbcf constant


disp '******************'
disp '** question 6 **'
disp '******************'

DISP '**test de Breush et Pagan**'
*compare le modèle avec effet individuel aléatoire et le pooled*

LINREG(NOPRINT) dlogva / resids
# logva{1} logfbcf constant
COMPUTE SCRT = %RSS
SET residmco = resids
PANEL(ENTRY=0.0,INDIV=per) residmco / residm
SET period = %period(T)==1
SMPL(SERIES=period)
SET scrim = residm**2
STATISTICS(noprint) scrim
COMPUTE scrm=sec*%mean
COMPUTE LM = ((sec*(per-1))/(2*((per-1)-1)))*(((scrm/scrt)-1)**2)
CDF CHISQUARED LM 1
smpl

disp '******************'
disp '** question 7 **'
disp '******************'

DISP ' **décomposition de la variance de la série de productivité**'
STAT(NOPRINT) va
COMPUTE vartot = %variance
PANEL(ENTRY=1.0,INDIV=-1.0) va / vaw
STAT(NOPRINT) vaw
COMPUTE varw = %variance
PANEL(ENTRY=0.0,INDIV=1.0) va / vab
STAT(NOPRINT) vab
COMPUTE varb = %variance
DISP 'within' varw/vartot
DISP 'between' varb/vartot