**AUTEUR : Isabelle CADORET-DAVID**

ALLOCATE 20
OPEN DATA surcout.wks
DATA(FORMAT=wks,ORG=obs) /
PRINT /
*************estimation du modèle à plusieurs équation apparemment non reliées********
SET logc = log(cost)
SET logy = log(output)
SET logpk = log(pk)
SET logpl = log(pl)
SET logpf = log(pf)
EQUATION cout logc
# CONSTANT logy logpk logpl logpf

EQUATION partk sk
# CONSTANT

EQUATION partl sl
# CONSTANT


EQUATION partf sf
# CONSTANT
SUR(iters=100) 4
# cout
# partk
# partl
# partf

RESTRICT(REPLACE) 3
# 3 6
# 1 -1 0
# 4 7
# 1 -1 0
# 5 8
# 1 -1 0
SUR(CREATE) 4
# cout
# partk
# partl
# partf
**cette estimation ne peut être réalisée car la matrice de variance-covariance est singulière**


************estimation en enlevant une équation******************************
SET logcpm = log(cost/Pl)
SET logy = log(output)
SET logpkpl = log(pk/pl)
SET logpfpl = log(pf/pl)

EQUATION cout logcpm
# CONSTANT logy logpkpl logpfpl

EQUATION partk sk
# CONSTANT

EQUATION partf sf
# CONSTANT

DISP '**estimation avec les MCO**'
LINREG logcpm
# CONSTANT logy logpkpl logpfpl

LINREG sk
# CONSTANT

LINREG sf
# CONSTANT


DISP '**estimation avec la méthode de Zellner itérative avec prise en compte des contraintes interéquations**'
SUR(iters=100) 3
# cout
# partk
# partf

RESTRICT(REPLACE) 2
# 3 5
# 1 -1 0
# 4 6
# 1 -1 0
SUR(CREATE) 3
# cout
# partk
# partf

COMPUTE bl = 1 - %beta(5) - %beta(6)
DISP 'bl' bl

DISP '**calcul de l"élasticité d"échelle **'
** cf Chambers : "Applied production Analysis" sur ce point**

COMPUTE echelle = 1/%beta(2)
DISP 'echelle' echelle