Puis4Asym1Prop=function(N,P1,P2,Seuil,PP) # Appel de procédure non interactive avec la taille de l'échantillon N, les valeurs de P1 et P2 # pour l'alternative H_0^{(4)}={P <= P1, ou P2 <= P} contre H_1^{(4)}={P1 < P < P2}, le seuil de signification du test et le PP # où calculer une estimation de la puissance. { Quant=qnorm(1-Seuil,lower.tail=TRUE) # Quantile d'ordre 1-alpha de la loi N(0 ; 1). Puis=pnorm(sqrt(N)*(P2-PP)/sqrt(PP*(1-PP))-Quant*sqrt((P2*(1-P2))/(PP*(1-PP))))-pnorm(sqrt(N)*(P1-PP)/sqrt(PP*(1-PP))+Quant*sqrt((P1*(1-P1))/(PP*(1-PP)))) # Calcul de l'approximation de la puissance asymptotique au point PP. options(digits=4) # Format d'affichage des résultats. return(Puis) # Renvoi de l'approximation de la puissance asymptotique. }