Bioequivalence and Bioavailability Forum

Main page Policy/Terms of Use Abbreviations Latest Posts

 Log-in |  Register |  Search

Back to the forum  Query: 2017-11-24 13:52 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: 5,794

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,495 Posts in 3,746 Threads, 1,085 registered users;
27 users online (0 registered, 27 guests).

The belief that there is only one truth and
that oneself is in possession of it,
seems to me the deepest root of all
that is evil in the world.    Max Born

BEBAC Ing. Helmut Schütz