library(ggplot2)

testWith <- function(n, mu, tst)
{
        pvals <- replicate(1000, tst(rnorm(n, 0, 1), rnorm(n, mu, 1))$p.value)
        data.frame(mean(pvals), mean(pvals < 0.05))
}

tbl <- data.frame()
for(mu in seq(-2, 2, 0.02)) for(n in seq(10, 100, 5)) { tbl <- rbind(tbl, data.frame(mu=mu, n=n)) }
wres <- cbind(tbl, do.call(rbind, apply(tbl, c(1), function(row) testWith(row[2], row[1], wilcox.test))))
tres <- cbind(tbl, do.call(rbind, apply(tbl, c(1), function(row) testWith(row[2], row[1], t.test))))

svg("tres_pvals.svg")
ggplot(tres, aes(x=n, y=mu)) + geom_tile(aes(fill=mean.pvals.))
dev.off()
svg("tres_power.svg")
ggplot(tres, aes(x=n, y=mu)) + geom_tile(aes(fill=mean.pvals...0.05.))
dev.off()

svg("wres_pvals.svg")
ggplot(wres, aes(x=n, y=mu)) + geom_tile(aes(fill=mean.pvals.))
dev.off()
svg("wres_power.svg")
ggplot(wres, aes(x=n, y=mu)) + geom_tile(aes(fill=mean.pvals...0.05.))
dev.off()