Puis1aAsym2MoyeInde=function(Donnees1,Donnees2,Delta_0,Seuil,Delta) # Appel de la procédure avec les vecteurs contenant les observations, la valeur de Delta_0 pour l'aternative. # H^{(1a)}_0={Delta_m <= Delta_0} contre HH^{(1a)}_01={Delta_0 < Delta_m}, le Seuil de signification du test et Delta le point où estimer la puissance. { N1=length(Donnees1) # Nombre d'observations du premier échantillon. N2=length(Donnees2) # Nombre d'observations du deuxième échantillon. Moyenne1=mean(Donnees1) # Moyenne des observations du premier échantillon. Moyenne2=mean(Donnees2) # Moyenne des observations du deuxième échantillon. Variance1=var(Donnees1) # Variance corrigée des observations du premier échantillon. Variance2=var(Donnees2) # Variance corrigée des observations du deuxième échantillon. EcartType1=sd(Donnees2) # Ecart type corrigé des observations du premier échantillon. EcartType2=sd(Donnees2) # Ecart type corrigé des observations du deuxième échantillon. Difference=Moyenne1-Moyenne2 # Différence des deux moyennes observées. EcartTypeDifference=sqrt(Variance1/N1+Variance2/N2) # Ecart type de la différence des moyennes observées. Quant=qnorm(1-Seuil,lower.tail=TRUE) # Quantile d'ordre 1-alpha de la loi N(0 ; 1). Puissance=1-pnorm(Quant + (Delta_0-Delta)/EcartTypeDifference) # Estimation de l'approximation de la puissance asymptotique du test au point Delta. options(digits=4) # Format d'affichage des résultats. return(Puissance) }