Puis2Asym2MoyeInde=function(Donnees1,Donnees2,Delta_0,Seuil,Delta) # Appel de la procédure avec le vecteur contenant les observations, Delta1 et Delta2, Delta_0 pour # l'aternative H^{(2)}_0={Delta_m = Delta_0} contre H^{(2)}_1={Delta_m non = Delta_0}, le seuil de signification du test et le point Delta 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. 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. Puissance=1-pnorm((Delta_0-Delta)/EcartTypeDifference+qnorm(1-Seuil/2))+pnorm((Delta_0-Delta)/EcartTypeDifference-qnorm(1-Seuil/2)) # Calcul de l'estimation de l'approximation de la puissance. options(digits=4) # Format d'affichage. return(Puissance) }