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

posted by d_labes  – Berlin, Germany, 2020-11-28 12:20 (1238 d 23:43 ago) – Posting: # 22091
Views: 5,809

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
22,988 posts in 4,825 threads, 1,657 registered users;
102 visitors (0 registered, 102 guests [including 9 identified bots]).
Forum time: 13:04 CEST (Europe/Vienna)

The whole purpose of education is
to turn mirrors into windows.    Sydney J. Harris

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