testGaussien <- function(x,y,a){ m <- length(x); n <- length(y); sqrt(m*n/(m+n))*(mean(x)-mean(y)) > qnorm(1-a) } testMW <- function(x,y,a){ te <- wilcox.test(x,y, paired=FALSE, alternative="greater"); te$p.valuepuis], ly[lp>puis], col="red", xlim=c(0,n+1), ylim = c(0, m+1)); points(lx[lp<=puis], ly[lp<=puis], col="blue"); } m <- 50; n <- 50; a <- 0.025 # PG <- puissanceOfSize(m, n, testGaussien, a); # write.table(as.table(PG), file="PuissanceTestGaussien.html") PG <- as.matrix(read.table("PuissanceTestGaussien.html")) # PMW <- puissanceOfSize(m, n, testMW, a); # write.table(as.table(PMW), file="PuissanceTestMW.html") PMW <- as.matrix(read.table("PuissanceTestMW.html")) par(mfrow=c(1,2)); plotPuissance(PG) plotPuissance(PMW) Im <- seq(1,m,0.1) # lines(Im, 4*Im/(Im - 4), col="green", lwd=3)