# NOT RUN {
sim.1 <- pwrFDR(effect.size=0.8, n.sample=60, lambda=0.90, r.1=0.05, N.tests=450, alpha=0.15,
method="sim", FDP.control.method="BHCLT")
sim.2 <- pwrFDR(effect.size=0.8, n.sample=60, lambda=0.90, r.1=0.05, N.tests=450, alpha=0.15,
method="sim", FDP.control.method="both", control=list(sim.level=2))
with(detail(sim.1)$reps, cbind(R.st/100, nna(R.R)/450))
with(detail(sim.2)$reps, cbind(R.st/100, nna(R.R)/450))
# }
Run the code above in your browser using DataLab