Test1aBern1p=function(K,N,P0,Seuil) # Appel de la procédure avec le nombre de succès ou de "un" K, la taille de l'échantillon N, # la valeur de P0 pour l'aternative H^{1a}_0={p <= p0} contre H^{1a}_1={p0 < p} et le seuil de signification Alpha du test. { Frequence=K/N # Proportion observée. Calpha=N+1 # Valeur initiale de la valeur critique. Palpha=0 # baleur initiale de la probabilité de H^(1a)_1. U=runif(1) # Simulation d'une loi U(]0 1[). while(Palpha < Seuil) # Condition d'exécution de la boucle pour la recherche du plus petit Calpha réalisant la condition. { Calpha=Calpha-1 # Nouvelle valeur critique. C'est la valeur critique en sortie de boucle. Palpha=pbinom(Calpha,N,P0,lower.tail=FALSE) # Calcul de la nouvelle probabilité de H^(1a)_1. } Calpha=Calpha+1 Gammaalpha=(Seuil-pbinom(Calpha,N,P0,lower.tail=FALSE))/dbinom(Calpha,N,P0) # Probabilité critique. if((Calpha