Données <- read.table("https://husson.github.io/livre_exo/chap7_ex7.csv",header=TRUE,sep=";", encoding = "latin1")

lambda <- sum(Données[,1]*Données[,2]) / sum(Données[,2])

Theo <- dpois(0:8,lambda)*sum(Données[,2])
Données <- cbind.data.frame(Données,Theo=Theo)

## regroupement des dernières classes
Données[5,2] <- sum(Données[5:9,2])
Données[5,3] <- sum(Données[5:9,3])
Données <- Données[1:5,]

(contribution <- (Données[,2]-Données[,3])^2/Données[,3])
(chi <- sum(contribution))
pchisq(chi,4,lower.tail=FALSE)