BinomInterPred=function(N,P,Alpha) # appel de la procédure avec les paramètres N et P de la loi B(N ; P) et le seuil de signification de l'intervalle de prédiction { Proba=1 # initialisation de la variable qui contiendra la probabilité de l'intervalle K1=0 # initialisation de la borne inférieure de l'intervalle - modalité à gauche K2=N # initialisation de la borne supérieure de l'intervalle - modalité à droite Arret=FALSE # indicateur de fin de procédure while(Arret==FALSE) { Pr1=dbinom(K1,N,P) # probabilité de la modalité à gauche Pr2=dbinom(K2,N,P) # probabilité de la modalité à droite if (Pr1>=Pr2) # comparaison des probabilités { # celle de la modalité à droite est la plus petite I=2 # indicatrice de cette modalité Proba=Proba-Pr2 # la probabilité de l'intervalle est diminué de la plus petite des deux probabilités } else { # celle de la modalité à gauche est la plus petite I=1 # indicatrice de cette modalité Proba=Proba-Pr1 # la probabilité de l'intervalle est diminué de la plus petite des deux probabilités } if (Proba>1-Alpha) { # le seuil de confiance n'est pas atteint if (I==1) { # la dernière probabilité retranchée est celle de la modalité à gauche K1=K1+1 # changement de cette modalité à gauche } else { # la dernière probabilité retranchée est celle de la modalité à droite K2=K2-1 # changement de cette modalité à droite } } else { # le seuil de confiance est atteint Arret=TRUE } } options(digits=7) # format d'affichage cat("Intervalle de prédiction de la loi B(",N,";",P,"),\n") # affichage des résultats cat("Seuil de l'intervalle :",100*Alpha,"%,\n") cat("Intervalle : [",K1,";",K2,"].\n") }