Test1aBern1p_Int=function() # Appel de la procédure interactive pour le test de l'aternative H^{1a}_0={p <= p0} contre H^{1a}_1={p0 < p}. { K=as.integer(readline("Saisir le nombre de «un» : ")) # Saisie du nombre de «un». N=as.integer(readline("Saisir la taille de l'échantillon : ")) # Saisie de la taille de l'échantillon. P0=as.numeric(readline("Saisir la proportion P0 : ")) # Saisie de la proportion P0 définissant les hypothèses de l'alternative. Seuil=as.numeric(readline("Saisir seuil du test : ")) # Saisie du seuil 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