###############################
#    boucle_binomiale.txt     #
###############################
# Orsay, le 26/11/2008
# Exercice sur la Loi binomiale
# représentation de B(n,p)avec n et p variables
# x, VA (ensemble de definition, ensemble des possibles de la VA)
# Nous allons utiliser une boucle "for" pour faire varier n et voir les résultats 
# sur les graphes successifs (une pause sera appliquée entre chaque graphe 
# avec l'option par(ask=TRUE) en attente de la frappe d'une touche pour 
# continuer ou, mieux encore avec Sys.sleep(temps) pour introduire une pause donnée.
min=1
max=10
k=5
p=0.3
# gestion automatique à chaque graphe rencontré
# par(ask=TRUE)
for (i in min:max)
{
  n=i*k
  x=0:n
  px=dbinom(x,n,p)
  titre=paste("Loi binomiale B(",n,",",p,")","  n*p=",n*p," et n*(1-p)=",n*(1-p))
  plot(px~x,type="h",xlim=c(0,max*k),ylim=c(0,0.35),lwd=4,col="red",xlab="X=k",ylab="P(X=k)",las=1,main=titre)
  lines(x,px,col="orange",lwd=2)
  # calcul de l'espérence de X~B(5,0.5)
  m=n*p
  # calcul de la variance de X~B(5,0.5)
  s2=n*p*(1-p)
  s=s2^0.5
  result=paste("n =",n,"=> moy =",m,"=> ect =",s)
  print(result)
  # utilisation de  Sys.sleep(temps) pour introduire une pause 
  # (les temps sont en seconde)
  Sys.sleep(1.5)
#  par(ask=TRUE)/ après chaque graphe/ c'est une autre solution mais ne change rien
}

#*********************
#* Fin du traitement *
#*********************