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

posted by d_labes  – Berlin, Germany, 2020-11-28 12:20 (2037 d 14:50 ago) – Posting: # 22091
Views: 17,126

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:

UA Flag
Activity
 Admin contact
23,655 posts in 4,993 threads, 1,572 registered users;
355 visitors (0 registered, 355 guests [including 34 identified bots]).
Forum time: 04:10 CEST (Europe/Vienna)

The real struggle is not between the right and the left
but between the party of the thoughtful
and the party of the jerks.    Jimmy Wales

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