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

posted by d_labes  – Berlin, Germany, 2020-11-28 10:20 (99 d 10:00 ago) – Posting: # 22091
Views: 1,594

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,371 posts in 4,463 threads, 1,496 registered users;
online 6 (0 registered, 6 guests [including 3 identified bots]).
Forum time: Sunday 20:20 UTC (Europe/Vienna)

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

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