Puis1aAsym2MoyeAppa=function(Donnees1,Donnees2,Delta0,Seuil,Delta_D) # Appel de la procédure avec les objets Donnees1 et Donnees2 contenant les observations, # Delta0 pour l'aternative H^{(1a)}_0={Delta_D =< Delta0} contre H^{(1a)}_1={Delta0 < Delta_D}, le seuil de signification du test et Delta_D où calculer la puissance. { 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(1-Seuil,lower.tail=TRUE) # Quantile d'ordre 1-alpha de la loi N(0 ; 1). Puissance=1-pnorm(Quant+sqrt(N1)*(Delta0-Delta_D)/EcartType) # Estimation de l'approximation de la puissance du test au point Delta_D. return(Puissance) # Retour du résultat. } }