Test2Asym1Prop_Int=function() # Appel de la procédure interactive pour l'aternative H_0^{(2)}={P = P0} contre H_1^{(2)}={P non = P0}. { K=as.integer(readline("Saisir le nombre «un» : ")) # Saisie du nombre de succès. 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 le seuil du test : ")) # Saisie du seuil du test. Frequence=K/N # Fréquence observée. Quant=qnorm(1-Seuil/2) # quantile d'ordre 1-alpha/2 de la loi N(0 ; 1). Condition=as.integer(1+max(c(7,Quant*sqrt(P0/(1-P0)),Quant*sqrt((1-P0)/P0))))**2 # Calcul de la condition de cohérence. Critique1=P0-Quant*sqrt(P0*(1-P0)/N) # Approximation de la première valeur critique du test. Critique2=P0+Quant*sqrt(P0*(1-P0)/N) # Approximation de la deuxième valeur critique du test. Pval=2*(1-pnorm(sqrt(N)*abs(Frequence-P0)/sqrt(P0*(1-P0)))) # p-valeur du test. if((Critique1 <= Frequence)&(Frequence <= Critique2)) # Décision. { Resultat="«H_0^{(2)}={p =" } else{ Resultat="«H_1^{(2)}={p not =" } PuisAPoster=1-pnorm(sqrt(N)*(P0-Frequence)/sqrt(Frequence*(1-Frequence))+Quant*sqrt((P0*(1-P0))/(Frequence*(1-Frequence))))+ pnorm(sqrt(N)*(P0-Frequence)/sqrt(Frequence*(1-Frequence))-Quant*sqrt((P0*(1-P0))/(Frequence*(1-Frequence)))) # Estimation de la puissance a posteriori. options(digits=4) # Format d'affichage. cat("\nTest asymptotique de l'alternative : H_0^{(2)}={ p =",P0,"} contre H_1^{(2)}={ p non =",P0,"}.\n") # Affichage des résultats. cat("Condition de cohérence du test : ",Condition,".\n") cat("Taille de l'échantillon :",N,"; nombre de «un» :",K,"; proportion observée :",Frequence,".\n") cat("Seuil asymptotique du test :",Seuil,"; valeurs critiques :",Critique1," et",Critique2,".\n") cat("P-valeur du test 1a :",Pval,".\n\n") cat("Décision :",Resultat,P0,"} est vraie».\n\n") cat("Approximation de la puissance a posteriori du test :",PuisAPoster,".\n") }