Bioequivalence and Bioavailability Forum

Hi All,

I was trying to produce sample size plots of full replicate BE for EMA regulation for my colleague, similar to those in this post, that would cover ISCV from 0.1 to 1, GMR from 0.9 to 1, using Helmut's code.

I noticed that at certain combinations PowerTOST would produce errors. e.g., with the following values
sampleN.scABEL(alpha = 0.05, targetpower = 0.8, theta0 = 0.91,
                    CV = 0.1, design = "2x2x4", regulator = "EMA")

would produce error
Error in if (pwr <= targetpower) cat(n, " ", formatC(pwr, digits = pd,  :
  missing value where TRUE/FALSE needed

These error only appear when CV is extremely low so in practice it doesn't matter since no one will use 2x2x4 replicate study for product with ISCV of 10%. I guess that's the reason I never noticed it before. So just for education purpose, can anyone tell me what's going on here? Maybe it's apparent for R gurus but as R novice I have no idea what the message says. According to Detlew, it's not a bug for sure.

---

Edit: Category changed. [Helmut]