Bioequivalence and Bioavailability Forum

Main page Policy/Terms of Use Abbreviations Latest Posts

 Log-in |  Register |  Search

Back to the forum  Query: 2017-03-29 01:31 CEST (UTC+2h)
 

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

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

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,
Shuanghe

Complete thread:

Back to the forum Activity
 Mix view
Bioequivalence and Bioavailability Forum | Admin contact
16,703 Posts in 3,586 Threads, 1,029 registered users;
14 users online (0 registered, 14 guests).

Most scientists today are devoid of ideas, full of fear, intent on
producing some paltry result so that they can add to the flood
of inane papers that now constitutes “scientific progress”
in many areas.    Paul Feyerabend

The BIOEQUIVALENCE / BIOAVAILABILITY FORUM is hosted by
BEBAC Ing. Helmut Schütz
XHTML/CSS RSS Feed