Power of the GCC framework and power of PowerTOST [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2020-11-28 11:20 (298 d 00:59 ago) – Posting: # 22091
Views: 2,893

Dear Helmut,

» Stolen the code of subject data sims and scaled ABE evaluation in the working horse function .pwr.SABE.sds().
» Have a look into the code of the new but undocumented function power.fwl.sds() in the GitHub repository.

there is no need to define and program a new function :cool:!
Try the following:
# define a new regulator object
regGCC <- reg_const("user", r_const= log(1/0.75)/CV2se(0.3), CVswitch = 0.3, CVcap = 0.3,
                    pe_const=TRUE)
regGCC
USER defined regulatory settings
- CVswitch            = 0.3
- cap on scABEL if CVw(R) > 0.3
- regulatory constant = 0.9799758
- pe constraint applied
# verify if the limits are what they should be
scABEL(CV=0.2, regulator = regGCC)
lower upper
 0.80  1.25
scABEL(CV=0.3, regulator = regGCC)
lower upper
 0.80  1.25
scABEL(CV=0.31, regulator = regGCC)
   lower    upper
0.750000 1.333333
scABEL(CV=0.35, regulator = regGCC)
   lower    upper
0.750000 1.333333
# use the function already in the package
power.scABEL.sds(CV=0.3, n=40, design="2x2x4", regulator = regGCC, theta0 = 1.25)
[1] 0.20463

Result matches the previous calculations.

Regards,

Detlew

Complete thread:

Activity
 Admin contact
21,691 posts in 4,534 threads, 1,541 registered users;
online 6 (0 registered, 6 guests [including 4 identified bots]).
Forum time: Wednesday 13:19 CEST (Europe/Vienna)

With four parameters I can fit an elephant,
and with five I can make him wiggle his trunk.    John von Neumann

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5