# Orsay le 28/11/2008
# exercice séries de manips répétées 5 fois
# 1 - importation des observations enregistrèes dans le fichier manips.txt 
#     (distributions observées)
distribs=read.table("./manips.txt",h=TRUE)
distribs
attach(distribs)
names(distribs)
noserie
nrepetitions
# 2 - selection d'une distribution observée particulière
# i, indice de la distribution
i=2
distrib_obs=c(f0[i],f1[i],f2[i],f3[i],f4[i],f5[i])
distrib_obs
# 3 - graphe de la distribution observée
#     comparée à la distribution theorique
# division de la fenetre graphique en 2 parties grace à l'option mfrow
par(mfrow=c(1,2))
etendue=0:5
ymax=0.35
titre1=paste("distribution observé No :",i)
plot(distrib_obs~etendue,type="h",col="blue",main=titre1,xlab="X=k",ylab="fobs",lwd=3,ylim=c(0,ymax))
# ditribution binomiale de référence
distrib_th=dbinom(etendue,5,0.5)
titre2="binomiale B(5,0.5) de référence"
plot(distrib_th~etendue,type="h",col="red",main=titre2,xlab="X=k",ylab="P(X=k)",lwd=3,ylim=c(0,ymax))
# Fin de traitement