Puis1bAsym2MoyeAppa=function(Donnees1,Donnees2,Delta0,Seuil,Delta_D) # 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={Delta_D < Delta0}, le seuil de signification du test et le point # Delta_D où calculer une estimation de l'approximation de la puissance asymptotique. { 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 { options(digits=4) # Format d'affichage des résultats. Donnees=Donnees1-Donnees2 # Objet contenant la différence des observations appariées. EcartType=sd(Donnees) # Ecart type des différences. Quant=qnorm(Seuil,lower.tail=TRUE) # Quantile d'ordre alpha de la loi N(0 ; 1). Puissance=pnorm(Quant+sqrt(N1)*(Delta0-Delta_D)/EcartType) # Estimation de l'approximation de la puissance asymptotique du test au point Delta_D. return(Puissance) } }