Approximate Power [🇷 for BE/BA]
Dear all,
Dieter Hauschke informed me that Jones’/Kenward’s code uses Satterthwaite’s approximation for degrees of freedom (for heterogenous variances), which was corrected in a recent book.
In power calculations degrees of freedom are fixed by the number of subjects (and sequences).
If you want to play around again, please replace code line #22
Power (CV=20%) by the (old) new codes for n=24 are:
Power >80% (CV=20%) by the (old) new codes are:
Power is >80% for all combinations tested.
According to Dieter Hauschke Owen’s exact method is implemented in nQuery Advisor.
Since there is still some confusion about different methods for sample size estimation in bioequivalence (software, tables, approximations) an entire chapter of a new book
Dieter Hauschke informed me that Jones’/Kenward’s code uses Satterthwaite’s approximation for degrees of freedom (for heterogenous variances), which was corrected in a recent book.
[1]
In power calculations degrees of freedom are fixed by the number of subjects (and sequences).
If you want to play around again, please replace code line #22
df <- sqrt(n2)*(nc1-nc2)/(2*t1)
df <- n2
Power (CV=20%) by the (old) new codes for n=24 are:
+======+=======+=======+
| GMR | old | new |
+------+-------+-------+
| 0.85 | 28.16 | 26.86 |
| 0.90 | 64.08 | 63.73 |
| 0.95 | 89.19 | 89.60 |
| 1.00 | 96.28 | 96.72 |
| 1.05 | 89.90 | 90.32 |
| 1.10 | 70.02 | 69.89 |
| 1.15 | 41.92 | 40.85 |
| 1.20 | 18.30 | 17.04 |
+======+=======+=======+
Power >80% (CV=20%) by the (old) new codes are:
+======+=====+===========+===========+
| GMR | n | old | new |
+------+-----+-----------+-----------+
| 0.85 | 134 | 0.8014178 | 0.8017723 |
| 0.90 | 38 | 0.8140704 | 0.8154940 |
| 0.95 | 20 | 0.8300156 | 0.8346802 |
| 1.00 | 16 | 0.8214263 | 0.8331980 |
| 1.05 | 18 | 0.7950343 | 0.8001854 |
| 1.10 | 32 | 0.8084890 | 0.8100682 |
| 1.15 | 72 | 0.8035456 | 0.8042181 |
| 1.20 | 294 | 0.8017617 | 0.8019247 |
+======+=====+===========+===========+
Power is >80% for all combinations tested.
❝ Maybe somebody of you has access to SAS or software specialized in power analysis (e.g., PASS or NQuery) and would like to check these results?
According to Dieter Hauschke Owen’s exact method is implemented in nQuery Advisor.
Since there is still some confusion about different methods for sample size estimation in bioequivalence (software, tables, approximations) an entire chapter of a new book
[2]
will be devoted to this issue.- S Patterson and B Jones
Bioequivalence and Statistics in Clinical Pharmacology
Chapman & Hall/CRC, Boca Raton, pp 230 (2006)
- D Hauschke, VW Steinijans and I Pigeot
Bioequivalence Studies in Drug Development: Methods and Applications
Wiley, New York (2007)
—
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
Helmut Schütz
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
Helmut Schütz
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Complete thread:
- R-Code for Power in 2×2 Cross-over Helmut 2006-12-28 21:33 [🇷 for BE/BA]
- Power in Diletti's sample size table Helmut 2007-01-01 23:59
- Approximate PowerHelmut 2007-01-06 15:17
- Power Calculation for replicate crossover design mathews 2008-04-03 13:19
- Power Calculation for replicate crossover design Helmut 2008-04-03 13:49
- Power Calculation for replicate crossover design mathews 2008-04-04 13:21
- Power Calculation for replicate crossover design Helmut 2008-04-03 13:49
- Exact Power Helmut 2009-09-25 19:20
- Exact Power ElMaestro 2009-09-27 23:27
- Software validation Helmut 2009-09-28 15:21
- All the important answers are here! ElMaestro 2009-09-28 19:12
- All the important answers are here! Helmut 2009-09-28 19:52
- All the important answers are here! ElMaestro 2009-09-28 20:18
- All the important answers are here! yjlee168 2009-09-28 21:30
- All the important answers are here! Helmut 2009-09-28 23:18
- 42! SASophylistic power oracle d_labes 2009-09-29 13:37
- Al Gore Rhythms Helmut 2009-09-29 14:52
- Al Gore Rhythms ElMaestro 2009-09-29 18:14
- Al Gore Rhythms yjlee168 2009-09-29 23:55
- Al Gore Rhythms Helmut 2009-09-29 14:52
- 42! SASophylistic power oracle d_labes 2009-09-29 13:37
- All the important answers are here! yjlee168 2009-09-29 23:49
- All the important answers are here! Helmut 2009-09-28 23:18
- All the important answers are here! Helmut 2009-09-28 19:52
- All the important answers are here! ElMaestro 2009-09-28 19:12
- Software validation Helmut 2009-09-28 15:21
- Exact Power ElMaestro 2009-09-27 23:27
- Power Calculation for replicate crossover design mathews 2008-04-03 13:19
- Approximate PowerHelmut 2007-01-06 15:17
- R-Code for Power in 2×2 Cross-over KDA 2019-01-23 04:24
- R-Code for Power in 2×2 Cross-over Ohlbe 2019-01-27 18:24
- R-Code for Power in 2×2 Cross-over Brus 2021-12-02 15:00
- R-Code for Power in 2×2 Cross-over ElMaestro 2021-12-02 15:15
- R-Code for Power in 2×2 Cross-over Helmut 2021-12-02 17:30
- R-Code for Power in 2×2 Cross-over Helmut 2021-12-02 16:57
- R-Code for Power in 2×2 Cross-over ElMaestro 2021-12-02 15:15
- Power in Diletti's sample size table Helmut 2007-01-01 23:59