Bioequivalence and Bioavailability Forum

Main page Policy/Terms of Use Abbreviations Latest Posts

 Log-in |  Register |  Search

Back to the forum  Query: 2018-02-22 21:56 CET (UTC+1h)

PowerTOST sample size EMA full replicate [R for BE/BA]

posted by Shuanghe - Spain, 2016-03-07 18:04  - Posting: # 16070
Views: 6,423

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. :-P

Edit: Category changed. [Helmut]

All the best,

Complete thread:

Back to the forum Activity
 Mix view
Bioequivalence and Bioavailability Forum | Admin contact
17,931 Posts in 3,818 Threads, 1,121 registered users;
17 users online (0 registered, 17 guests).

When people learn no tools of judgment
and merely follow their hopes,
the seeds of political manipulation are sown.    Stephen Jay Gould

BEBAC Ing. Helmut Schütz