Codes/Figures                        Données                         Compléments                     Erratum





Comp11. Production de lait



# 11.1. Exploration
# Lag plot

lait2=read.table(system.file("/import/collecteLait.txt",package="caschrono"),
 header=FALSE, sep=";",colClasses=c("character",rep("numeric",3)),dec=".",
col.names=c("mois","an","evol","coll.v","cum.v","coll.m","cum.m"))
lait= ts(lait2$coll.v/1000,start=c(1979,1), frequency = 12)
log.lait=log(lait)
laitdeb = window(lait, end=c(1983,12))
laitfin = window(lait, start=c(2004,1))
lag.plot(rev(laitdeb),set=c(1:12),pch="+",col="black",do.lines=FALSE)

Fig11-1b

lag.plot(rev(laitfin),set=c(1:12),pch="+",col="black",do.lines=FALSE)

Fig11-1c

# 11.3 Essai de modélisation après l'introduction des quota
apres95 = window(log.lait, start = c(1995,1))
plot(apres95)

Fig11-2b