42! SASophylistic power oracle [🇷 for BE/BA]

posted by d_labes  – Berlin, Germany, 2009-09-29 15:37 (5295 d 01:03 ago) – Posting: # 4267
Views: 38,767

Dear guys,

let me take part in your hairsplitting party :-D .
Hope you have already waited for me.

Here are the results of the SAS Proc Power.
Code can be found here, of course appropriately adapted.

                 The POWER Procedure
        Equivalence Test for Paired Mean Ratio

               Fixed Scenario Elements

       Distribution                   Lognormal
       Method                             Exact
       Lower Equivalence Bound              0.8
       Upper Equivalence Bound             1.25
       Alpha                               0.05
       Geometric Mean Ratio                   1
       Coefficient of Variation           0.075
       Correlation                            0


              Computed Power
                        N
            Index    Pairs    Power

               1        3    0.604
               2        4    0.867
               3        5    0.968
               4        6    0.993

The power is given with 3 decimals only by default. No simple way to change.

Seems the power is calculated in this case totally wrong!
This is astonishing enough, because the same Proc Power code can reproduce most of the sample sizes given in Diletti et.al., at least for power 80%.
Ok, Proc Power is called experimental by SAS, usually means the user shall do the testing and report bugs. Others call such software beta.

My own rolled code (relying on Owen's Q-functions, undocumented but available in SAS) gives:
        Coeff. of
        variation   'true'  lower BE  upper BE        achieved   required
 alpha     (%)      µ(T/R)    bound    bound     N     power      power

  0.05     7.5        1       0.8      1.25      4     0.729014    0.7
  0.05     7.5        1       0.8      1.25      6     0.986992    0.7
  0.05     7.5        1       0.8      1.25      8     0.999577    0.7

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
113 visitors (0 registered, 113 guests [including 7 identified bots]).
Forum time: 15:40 CET (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