# NOT RUN {
powChapman(nullN=500, trueN=400, n1=100, n2=100, nsim=1000)
Ntotry <- seq(from=250, to=450, by=25)
pows <- sapply(Ntotry, function(x)
powChapman(nullN=500, trueN=x, n1=100, n2=100, nsim=1000))
plot(Ntotry, pows)
# }
Run the code above in your browser using DataLab