# ========================================================
# Exemples d'estimations de modèle à équations simultanées
# ========================================================

# on utilise la commande "systemfit" du package "systemfit"
# ce qui nécessite qu'il soit préalablement installé

# lecture des données puis estimations du modèle de Klein
klein = read.table("data.csv", header=TRUE, sep=";", dec=",")
attach(klein)

klein # pour voir les données
fix(klein) # pour les éditer (refermer la fenêtre d'édition pour poursuivre)

library(systemfit) # pour charger le package "systemfit" supposé déjà installé

help(systemfit) # si on veut lire l'aide de la commande "systemfit"

# equations
conso =  cons ~ p + plag + wtot
wage = wpri ~ x + xlag + year
invest = inv ~ p + plag + klag
modklein = list(consommation = conso, salaires = wage, investissement = invest)

# estimations par les mco
fitols = systemfit(modklein)
summary(fitols)

# estimations par la méthode "sur"
fitsur = systemfit(modklein, "SUR")
summary(fitsur)

# estimations par les doubles moindres carrés
inst = ~ klag + plag + xlag + wpub + g + t + year
fitdmc = systemfit(modklein, "2SLS", inst)
summary(fitdmc)

# estimations par les doubles moindres carrés
# avec instruments différents
inst1 = ~ klag + plag + xlag + t + year
inst2 = ~ klag + plag + xlag + wpub + g + t + year
inst3 = ~ klag + plag + xlag + wpub + g + t
instlist = list(inst1, inst2, inst3)
fitdmc2 = systemfit(modklein, "2SLS", instlist)
summary(fitdmc2)

# estimations par les triples moindres carrés
inst = ~ klag + plag + xlag + wpub + g + t + year
fittmc = systemfit(modklein, "3SLS", inst)
summary(fittmc)