Bioequivalence and Bioavailability Forum

❝ Can't enough!

❝ One small correction: The graphics under "Sample size for RSABE" is apparent not correct here.

library(PowerTOST)
library(magicaxis)
runs  <- 5
nsims <- as.integer(10^(seq(3, 6, 0.1)))
res   <- data.frame(nsims = nsims)
for (k in seq_along(nsims)) {
  for (j in 1:runs) {
    ifelse (j == 1, setseed <- TRUE, setseed <- FALSE)
    res[k, j+1] <- sampleN.RSABE(theta0 = 0.9, CV = 0.7, nsims = nsims[k],
                                 setseed = setseed, design = "2x2x4",
                                 targetpower = 0.9, print = FALSE,
                                 details = FALSE)[["Sample size"]]
  }
}
ylim <- range(res[, 2:(runs+1)], na.rm = TRUE)
dev.new(width = 5, height = 5)
op <- par(no.readonly = TRUE)
par(pty = "s")
plot(res$nsims, res[, 1], type = "n", log = "x", ylim = ylim, axes = FALSE,
     xlab = "Number of simulations", ylab = "Sample size")
grid(); box()
magaxis(1)
magaxis(2, minorn = 2, las = 1)
for (j in 1:runs) {
  ifelse (j == 1,  pch <- 16, pch <- 1)
  points(res$nsims, res[, j+1], pch = pch, col = "blue", cex = 1.15)
}
points(1e5, 48, pch = 16, col = "#00AA00", cex = 1.15) # default
points(1e4, 46, pch = 16, col = "red", cex = 1.15)     # Lászlós
par(op)
print(res, row.names = FALSE)

   nsims V2 V3 V4 V5 V6
    1000 50 48 54 52 50
    1258 50 48 52 50 52
    1584 50 50 52 50 56
    1995 50 48 50 50 50
    2511 48 48 48 48 52
    3162 48 50 50 50 50
    3981 46 48 50 50 48
    5011 46 50 50 50 48
    6309 46 48 48 48 50
    7943 46 50 48 48 50
   10000 46 50 48 50 50
   12589 46 50 48 50 48
   15848 48 48 48 48 48
   19952 48 48 50 50 48
   25118 48 48 50 48 48
   31622 48 48 48 48 48
   39810 48 48 50 48 50
   50118 48 48 48 48 48
   63095 48 48 48 48 48
   79432 48 48 48 48 48
  100000 48 48 48 48 48
  125892 48 48 48 48 48
  158489 48 48 48 48 48
  199526 48 48 48 48 48
  251188 48 48 48 48 48
  316227 48 48 48 48 48
  398107 48 48 48 48 48
  501187 48 48 48 48 48
  630957 48 48 48 48 48
  794328 48 48 48 48 48
 1000000 48 48 48 48 48