pvalores <- NA
pvaloresW <- NA
podert <- NA
poderw <- NA
tamanho <- c(6,10,20,30,60,100,200)
for (i in 1:length(tamanho)) {
for (simulacao in 1:1000) {
  grupo1 <- round(rnorm(n = tamanho[i], mean = 26, sd = 5), digits = 2)
  grupo2 <- round(rnorm(n = tamanho[i], mean = 28, sd = 5), digits = 2)
  a <- t.test(x = grupo1, y = grupo2)
  b <- wilcox.test(x = grupo1, y = grupo2)
  pvalores[simulacao] <- a$p.value
  pvaloresW[simulacao] <- b$p.value
}
podert[i] <- sum(pvalores <= 0.05)/1000
poderw[i] <- sum(pvaloresW <= 0.05)/1000
}

plot(x = tamanho, y = podert)
lines(x = tamanho, y = podert, add = TRUE)
lines(x = tamanho, y = poderw, col = 'red')