load("/Users/carraro/Documents/Laurent/enseignement/axe_MSA/06-07/cartons/cartons.RData")
#création du data frame dechets
plot(dechets)
#pb avec les variables qualitatives prestataire et categorie
pairs(cartons~surface+CA, data=dechets,labels=c("cartons en tonnes","surface en m^2","CA en keuros"))
boxplot(cartons~categorie+prestataire, data=dechets,main="tonnage cartons")
#Premières modélisations
abscisses<-min(dechets$CA)+1:500*(max(dechets$CA)-min(dechets$CA))/500
mod1<-lm(cartons~CA,data=dechets)
p1<-predict(mod1,newdata=data.frame(CA=abscisses))
mod2<-lm(cartons~CA+I(CA^2),data=dechets)
p2<-predict(mod2,newdata=data.frame(CA=abscisses))
mod3<-lm(cartons~CA+I(CA^2)+I(CA^3),data=dechets)
p3<-predict(mod3,newdata=data.frame(CA=abscisses))
plot(cartons~CA,data=dechets,xlab="CA en keuros",ylab="cartons en tonnes")
lines(abscisses,p1,col=2)
lines(abscisses,p2,col=3)
lines(abscisses,p3,col=4)
legend("topleft",lty=1,c("degré 1","degré 2","degré 3"),col=c(2,3,4))
mod_1<-lm(cartons~CA+surface+categorie+prestataire,data=dechets)
summary(mod_1)
plot(mod_1$residuals[prestataire=="A"]~CA[prestataire=="A"],data=dechets,main="résidus tonnage cartons",col="red",xlab="CA",ylab="résidus")
points(mod_1$residuals[prestataire=="B"]~CA[prestataire=="B"],data=dechets,main="résidus tonnage cartons",col="blue",pch="*",cex=1.5)
legend("topleft",pch=c("o","*"),c("prestataire A","prestataire B"),col=c("red","blue"),cex=1.5)

#étude pour prestataire A
dechetsA<-dechets[dechets$prestataire=="A",]
pairs(cartons~surface+CA, data=dechetsA,labels=c("cartons en tonnes","surface en m^2","CA en keuros"))
boxplot(cartons~categorie, data=dechetsA,main="tonnage cartons")
#on passe au log
pairs(log10(cartons)~log10(surface)+log10(CA), data=dechetsA,labels=c("log(cartons) en tonnes","log
(surface) en m^2","log(CA) en keuros"))
boxplot(log10(cartons)~categorie, data=dechetsA,main="tonnage cartons")
mod_log_A<-lm(log10(cartons)~log10(CA)+log10(surface)+categorie,data=dechetsA)
summary(mod_log_A)
plot(mod_log_A$residuals~CA,data=dechetsA,main="résidus tonnage cartons",xlab="CA",ylab="résidus")

#étude pour prestataire B
dechetsB<-dechets[dechets$prestataire=="B",]
pairs(cartons~surface+CA, data=dechetsB,labels=c("cartons en tonnes","surface en m^2","CA en keuros"))
boxplot(cartons~categorie, data=dechetsB,main="tonnage cartons")
#on passe au log
pairs(log10(cartons)~log10(surface)+log10(CA), data=dechetsB,labels=c("log(cartons) en tonnes","log
(surface) en m^2","log(CA) en keuros"))
boxplot(log10(cartons)~categorie, data=dechetsB,main="tonnage cartons")
mod_log_B<-lm(log10(cartons)~log10(CA)+log10(surface)+categorie,data=dechetsB)
summary(mod_log_B)
plot(mod_log_B$residuals~CA,data=dechetsB,main="résidus tonnage cartons",xlab="CA",ylab="résidus")