setwd("H:/MES DOCUMENTS/Mes sites Web/monsiteweb2/LivreEconometrie/programmes/R")
getwd()

## lecture des données
chap2<-read.table("chapitre02R.csv",head=TRUE,sep=";")
head(chap2)## visualise les premières observations

## problème 2.1
reg7098=lm(txi7098~txs7098,data=chap2)
summary(reg7098) ##  estimation MCO

aovreg7098<-aov(txi7098~txs7098,data=chap2)
aovreg7098   ## 	analyse de la variance

plot(txi7098~txs7098,data=chap2)
abline(reg7098)# trace la droite de regression

chap2_15=chap2[chap2$pays<16,]# selection de 15 pays
reg7098_15=lm(txi7098~txs7098,data=chap2_15)
summary(reg7098_15) ##  estimation MCO sur 15 pays

ttab<-qt(0.975,17)# intervalle de confiance
var1<-vcov(reg7098)
c(reg7098$coef[2]-ttab*sqrt(var1[2,2]),reg7098$coef[2]+ttab*sqrt(var1[2,2])) ## intervalle de confiance

confint(reg7098,level=0.95)# intervalle de confiance pour les 2 coefficients (plus rapide)

reg7079=lm(txi7079~txs7079,data=chap2)
summary(reg7079)
confint(reg7079,level=0.95)

reg8089=lm(txi8089~txs8089,data=chap2)
summary(reg8089)
confint(reg8089,level=0.95)


reg9098=lm(txi9098~txs9098,data=chap2)
summary(reg9098)
confint(reg9098,level=0.95)

## problème 2.2
reg1=lm(txi7098~txsprive7098+txspub7098,data=chap2)
summary(reg1)
reg11=lm(txi7098~I(txsprive7098+txspub7098),data=chap2)
summary(reg11)
anova(reg1,reg11)

reg2=lm(txi7079~txsprive7079+txspub7079,data=chap2)
summary(reg2)
reg22=lm(txi7079~I(txsprive7079+txspub7079),data=chap2)
summary(reg22)
anova(reg2,reg22)

reg3=lm(txi8089~txsprive8089+txspub8089,data=chap2)
summary(reg3)
reg33=lm(txi8089~I(txsprive8089+txspub8089),data=chap2)
summary(reg33)
anova(reg3,reg33)

reg4=lm(txi9098~txsprive9098+txspub9098,data=chap2)
summary(reg4)

reg44=lm(txi9098~I(txsprive9098+txspub9098),data=chap2)
summary(reg44)
anova(reg4,reg44)

reg42=lm(txi9098~txsprive9098,data=chap2)
res42<-residuals(reg42)
summary(reg42)
reg43=lm(txspub9098~txsprive9098,data=chap2)
res43<-residuals(reg43)
summary(reg43)

reg45=lm(res42~res43)
summary(reg45)


## problème 2.3
regDC=lm(dc7098~txs7098,data=chap2)
summary(regDC)

regDC1=lm(dc7098~txsprive7098+txspub7098,data=chap2)
summary(regDC1)
regDC2=lm(dc7098~I(txsprive7098+txspub7098),data=chap2)
summary(regDC2)
anova(regDC1,regDC2)