EstimatAsympDiff2MoyennesAppar=function(Donnees1,Donnees2,Seuil) # appel de procédure avec les Objets Donnees1 et Donnees2 contenant les observations et le seuil de signification de l'intervalle { N1=length(Donnees1) # nombre d'observations du premier échantillon N2=length(Donnees2) # nombre d'observations du deuxième échantillon options(digits=7) # format d'affichage if(N1==N2) { Donnees=Donnees1-Donnees2 # objet contenantles différences des observation 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) # écart type des observations du premier échantillon EcartType2=sd(Donnees2) # écart type des observations du deuxième échantillon EcartType=sd(Donnees) # écart type des différences calpha=qnorm(1-(Seuil/2),lower.tail=TRUE) # quantile d'ordre 1-alpha/2 de la loi N(0 ; 1) BorneInf=Moyenne-calpha*EcartType/sqrt(N1) # borne inférieure de l'intervalle asymptotique de confiance BorneSup=Moyenne+calpha*EcartType/sqrt(N1) # borne supérieure de l'intervalle asymptotique de confiance cat("\nEstimation asymptotique de la différence de deux moyennnes théoriques de v. a. appariées.\n") # affichage des résultats cat("\n\nPremier é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("\n\nEstimation de la moyenne théorique des différrences :",Moyenne,".\n") cat("Estimation de l'écart type théorique des différences :",EcartType,".\n") cat("\n\nSeuil asymptotique de l'intervalle de confiance :",Seuil,".\n") cat("L'intervalle asymptotique de confiance est : [",BorneInf,"; ",BorneSup,"].\n") } else { cat("\n Erreur sur les données : les échantillons n'ont pas la même taille\n") } }