Test2Asym1Moye=function(Donnees,Mu0,Seuil) # Appel de la procédure avec le vecteur contenant les observations, Mu0 pour l'aternative H_0^{(2)}={Mu = Mu0} contre H_1^{(2)}={Mu non = Mu0} et le seuil de # signification du test. { N=length(Donnees) # Nombre d'observations. Moyenne=mean(Donnees) # Moyenne des observations. EcartType=sd(Donnees) # Ecart type des observations. T=qnorm(1-Seuil/2) # Quantile d'ordre 1-seuil/2 de la loi Normale Standard. Critique1=Mu0 - (T*EcartType)/sqrt(N) # Estimation de la première valeur critique du test. Critique2=Mu0 + (T*EcartType)/sqrt(N) # Estimation de la deuxième valeur critique du test. if((Critique1 <= Moyenne)&( Moyenne <= Critique2)) # Décision. { Resultat="«H_0^{(2)}" } else{ Resultat="«H_1^{(2)}" } Pval=2*(1-pnorm(abs(Moyenne-Mu0)*sqrt(N)/EcartType)) # Estimation de la \(p-\)valeur. PuisAPoster=1-pnorm(sqrt(N)*(Mu0-Moyenne)/EcartType+T)+pnorm(sqrt(N)*(Mu0-Moyenne)/EcartType-T) # Estimation de la puissance a posteriori. options(digits=4) # Format d'affichage. cat("\nTest asymptotique de l'alternative : H_0^{(2)}={ Mu =",Mu0,"} contre H_1^{(2)}={ Mu non =",Mu0,"}.\n") # Affichage des résultats. cat("Taille de l'échantillon :",N,"; moyenne observée :",Moyenne,"; écart type observé : ",EcartType,".\n") cat("Seuil asymptotique du test :",Seuil,".\n") cat("Estimation des valeurs critiques :",Critique1," et",Critique2,".\n\n") cat("Décision :",Resultat,"est vraie».\n\n") cat("La p-valeur du test est :",Pval,".\n") cat("Estimation de la puissance a posteriori du test :",PuisAPoster,".\n") }