Test1bAsym2MoyeAppa=function(Donnees1,Donnees2,Delta0,Seuil) # Appel de la procédure avec les objets Donnees1 et Donnees2 contenant les observations, # la valeur de Delta0 pour l'aternative H^{(1b)}_0={Delta0 <= Delta_D} contre H^{(1b)}_1={DeltaD < Delta_0} et le seuil de signification du test. { N1=length(Donnees1) # Nombre d'observations du premier échantillon. N2=length(Donnees2) # Nombre d'observations du deuxième échantillon. if(N1!=N2) # Contrôle sur la taille des deux échantillons. { cat("\n Erreur sur les données : les échantillons n'ont pas la même taille\n") } else { Donnees=Donnees1-Donnees2 # Objet contenant la différence des observations appariées. Moyenne1=mean(Donnees1) # Moyenne des observations du premier échantillon. Moyenne2=mean(Donnees2) # Moyenne des observations du deuxième échantillon. Moyenne=mean(Donnees) # Moyenne des différences. EcartType1=sd(Donnees1) # Ecart type des observations du premier échantillon. EcartType2=sd(Donnees2) # Ecart type des observations du deuxième échantillon. EcartType=sd(Donnees) # Ecart type des différences. Quant=qnorm(Seuil,lower.tail=TRUE) # Quantile d'ordre alpha de la loi N(0 ; 1). Critique=Delta0+Quant*EcartType/sqrt(N1) # Estimation de la valeur critique du test. Pval=pnorm(sqrt(N1)*(Moyenne-Delta0)/EcartType) # Estimation de la p-valeur du test. if(Moyenne < Critique) # Résultats du test. { Resultat="«H^{(1b)}_1={ Delta_D < " } else { Resultat="«H^{(1b)}_0={ Delta_D >= " } options(digits=4) # Format d'affichage des résultats. cat("\nTest asymptotique de l'alternative : H^{(1b)}_0={",Delta0," <= Delta_D } contre H^{(1b)}_1={ Delta_D < ",Delta0,"}.\n\n") # Affichage des résultats. cat("Premier échantillon, taille :",N1,"; moyenne :",Moyenne1,"; écart type : ",EcartType1,".\n") # Caractéristiques du premier échantillon. cat("Deuxième échantillon, taille :",N2,"; moyenne :",Moyenne2,"; écart type : ",EcartType2,".\n") # Caractéristiques du deuxième échantillon cat("Moyenne observée des différences :",Moyenne,"; écart type observé des différences: ",EcartType,".\n\n") # Caractéristiques des différences. cat("Seuil asymptotique du test :",Seuil,"; estimation de la valeur critique :",Critique,".\n") cat("Estimation de la p-valeur :",Pval,".\n\n") cat("Décision :",Resultat,Delta0,"} est vraie».\n\n") if((Moyenne < Delta0)&(Moyenne>=Critique)) # L'estimation de l'approximation de la puissance asymptotique a posteriori n'est affichée que si la Moyenne observée # est dans H^{(1b)}_1 et la décision est «H^{(1b)}_0 vraie». { PuisAPoster=pnorm(Quant+sqrt(N)*(Delta0-Moyenne)/EcartType) # Estimation de l'approximation de la puissance asymptotique a posteriori. cat("Estimation de l'approximation de la puissance asymptotique a posteriori :",PuisAPoster,".\n") } } }