Anu ★ India, 2013-02-01 08:01 (4469 d 18:53 ago) Posting: # 9938 Views: 12,093 |
|
Hi All, Greetings. Kindly help me out with, how can we calculate sample size in case of RSABE according to US FDA Guidelines? |
d_labes ★★★ Berlin, Germany, 2013-02-01 10:46 (4469 d 16:08 ago) @ Anu Posting: # 9940 Views: 11,290 |
|
Dear Anu, there is no simple (or complicated) analytical way to calculate power or estimate the sample size for RSABE due to the complex conditions and requirements for the described methods (FDA or EMA) ![]() The only way known at moment is via simulations. Have a look into the paper: Laszlo Tothfalusi and Laszlo Endrenyi "Sample Sizes for Designing Bioequivalence Studies for Highly Variable Drugs" J Pharm Pharmaceut Sci (www.cspsCanada.org) 15(1) 73 - 84, 2011 Can be accessed online. See also this thread here in the forum. Hope this helps. — Regards, Detlew |
Anu ★ India, 2013-02-07 07:58 (4463 d 18:56 ago) @ d_labes Posting: # 9982 Views: 10,544 |
|
Hi Detlew, Thankyou so much for your reply. I have used the R-code as well that you have provided. In my data (CV=40%, Ratio=1.16), through R-code sample size came 39 and with ref to the Table A3, Sample-size is between (47-104). |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-07 13:32 (4463 d 13:22 ago) @ Anu Posting: # 9985 Views: 10,590 |
|
Hi Anu, ❝ I have used the R-code as well that you have provided. In my data (CV=40%, Ratio=1.16), through R-code sample size came 39 and with ref to the Table A3, Sample-size is between (47-104). As Detlew already pointed out there is no analytical solution for RSABE due to the GMR-restriction (+ the 50% CV cap for EMA). Either you use the tables or you set up your own simulations. require(PowerTOST) Table A3 (FDA method, 3-period, 80% power) for CV 40% gives 47 (GMR 1.15) and 104 (GMR 1.20) subjects. Sample sizes are substantially larger (especially for GMR 1.20) due to the [0.80–1.25] restriction.Good news: For GMR 1.16 your 39 would be correct …Bad news: … but only by ignoring the restriction.— Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
d_labes ★★★ Berlin, Germany, 2013-02-08 12:57 (4462 d 13:57 ago) @ Helmut Posting: # 9992 Views: 10,602 |
|
Dear Helmut, dear Anu, ❝ As Detlew already pointed out there is no analytical solution for RSABE due to the GMR-restriction (+ the 50% CV cap for EMA). Either you use the tables or you set up your own simulations. No need to set up your own simulations. Just let PowerTOST do the job ![]() sampleN.scABEL(CV=0.4, theta0=1.16, design="2x3x3", regulator="FDA") And the best of all: computational time for each step approx. 1-2 sec for 1E+06 sims, believe it or not ![]() But ... Cave 1: Only functions with PowerTOST V1.1-00. Distribution via CRAN may take some time. Cave 2: PowerTOST is doing the sample size estimation only for balanced designs since the break down of the total subject number in case of unbalanced sequence groups is not unique. Moreover the formulas used are only for balanced designs. That may give different sample sizes compared to the two Laszlos. Cave3: In case of regulator="FDA" the sample size is only approximate since the BE decision method via widened BE limits is not exactly what is expected by the FDA. But the two Laszlos state that the scABEL method should be 'operational' equivalent to the FDA method. Thus the sample size should be comparable.— Regards, Detlew |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-08 13:19 (4462 d 13:35 ago) @ d_labes Posting: # 9993 Views: 10,663 |
|
Dear Detlew! ❝ No need to set up your own simulations. Just let PowerTOST do the job ❝ And the best of all: computational time for each step approx. 1-2 sec for 1E+06 sims, believe it or not Do you want to marry me? The community will love you for that! Opens even the door for two-stage simulations where the intermediate power/sample size was the actual show-stopper. ❝ In case of Agree. Remember this one? Somewhere in electronic nirvana I have presentations by Donald Schuirmann showing scaled limits ‘for demonstrational purposes’ as well. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Anu ★ India, 2013-02-09 10:20 (4461 d 16:34 ago) @ Helmut Posting: # 9994 Views: 10,483 |
|
Hi Helmut & Detlew, Many thanks for your valueable replies. Please help me with one thing that while using the code sampleN.scABEL(CV=0.4, theta0=1.16, design="2x3x3", regulator="FDA") inspite of PowerTOST & mvtnorm, do we need to install someother package as well. The following error came while execution: Error: could not find function "sampleN.scABEL" Thanks & Regards |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-09 10:35 (4461 d 16:19 ago) @ Anu Posting: # 9995 Views: 10,533 |
|
Hi Anu, see Detlew’s first warning above. You need V1.1-00 of PowerTOST which is not available at CRAN yet. Edit 2013-02-10: Go ahead – it’s there. ![]() — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
d_labes ★★★ Berlin, Germany, 2013-02-09 19:11 (4461 d 07:43 ago) @ Helmut Posting: # 9996 Views: 10,486 |
|
Dear Helmut! ❝ Do you want to marry me? The community will love you for that! They would call me "Bigamist" and send me to prison in many countries. And I think my spouse wouldn't be very amused if I convert to a Mussulmán in my old days ![]() But anyhow: Thanks for the flowers. ❝ Opens even the door for two-stage simulations where the intermediate power/sample size was the actual show-stopper. I'm not quite sure if you are right here. In simulating 2-stage adaptive designs the power/sample size steps may contribute considerably to the computation times. Thus they may stay the show-stopper even if we can push on the calculations of stage 1/stage 2 BE by the ideas behind power.scABEL() or behind the code inthis post. But lets see ... 2-stage design combined with RSABE will have another drawback: Where is the point for simulations which represents alpha ![]() Eventually power.scABEL() will help to explore. — Regards, Detlew |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-10 19:35 (4460 d 07:19 ago) @ d_labes Posting: # 9997 Views: 10,815 |
|
Dear Detlew! ❝ They would call me "Bigamist" and send me to prison in many countries. And I think my spouse wouldn't be very amused if I convert to a Mussulmán in my old days I know. Was in the tradition of this post. ![]() ❝ 2-stage design combined with RSABE will have another drawback: Where is the point for simulations which represents alpha Duno. Of course the limits are flexible, but the GMR-restriction cuts in anyway. Have to think about it. In the following results of some testing: The sample size converges pretty fast if the ratio is not extreme. ![]() In the common range nsims=1e+05 should be sufficient.In some cases the sample size is lower than Lászlós’. For CV 0.5, ratio 0.85, 90% power, 3-period (EMA) I got n=150 even with 100 mio sim’s whereas in Table A1 I find n=158. Likely their number of simulations (10,000) was too low. EMA, 3-period, 80% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 EMA, 3-period, 90% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 EMA, 4-period, 80% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 EMA, 4-period, 90% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 PT: PowerTOST (higher due to balanced sequences, higher, lower), TE: Tóthfalusi/Endrényi — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-11 02:00 (4460 d 00:54 ago) @ d_labes Posting: # 9998 Views: 10,613 |
|
Continuing… FDA, 3-period, 80% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 FDA, 3-period, 90% power
0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 FDA, 4-period, 80% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 FDA, 4-period, 90% power 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 I will ask them whether they have used 0.89 like in their paper from 2009 instead of log(1.25)/0.25 ≈ 0.8925742… I like the option to specify unequal CVs of test and reference as CV=c(T,R) . Example for EMA, 3-period, T/R 0.95, 80% power (pooled CVintra 0.3 in all cases): CVWT CVWR n It will pay off to have a ‘better’ formulation. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
d_labes ★★★ Berlin, Germany, 2013-02-11 13:03 (4459 d 13:51 ago) @ Helmut Posting: # 10001 Views: 10,532 |
|
Dear Helmut! Thank you so much for this comprehensive testing ![]() Saves me a lot of work. Can you give me a statement on the validity of power.scABEL() and sampleN.scABEL() ?❝ I will ask them whether they have used 0.89 like in their paper from 2009 instead of log(1.25)/0.25 ≈ 0.8925742… Maybe some of the differences in case of regulator="FDA" are due to the widened ABEL approach instead of the approach via the linearized SABE criterion and its 95% upper confidence limit <0. Again: Coming soon ... or later ... ![]() But if you ask them: Ask in this respect if they have used "The unknown x" used in the SAS code of the progesterone guidance in their implementation of the linearized SABE criterion. I haven't found any hint in their papers on HVD/HVDP regarding such a 'correction' for the upper CI part of the point estimator. — Regards, Detlew |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-11 15:40 (4459 d 11:14 ago) @ d_labes Posting: # 10002 Views: 10,537 |
|
Dear Detlew! ❝ Can you give me a statement on the validity of Cross-validating two simulations? ![]() ❝ Maybe some of the differences in case of Yes, may be another cause. ❝ Again: Coming soon ... or later ... Relax. ❝ But if you ask them: Ask in this respect if they have used "The unknown x" used in the SAS code of the progesterone guidance in their implementation of the linearized SABE criterion. I haven't found any hint in their papers on HVD/HVDP regarding such a 'correction' for the upper CI part of the point estimator. Right guess, IMHO. Remember this thread? I’m short in time right now; if you want I can give you their e-mail addresses and you do it yourself. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
d_labes ★★★ Berlin, Germany, 2013-02-15 09:47 (4455 d 17:07 ago) @ Helmut Posting: # 10027 Views: 10,401 |
|
Dear Helmut, a follow up to your post: ❝ I like the option to specify unequal CVs of test and reference as ❝ ❝ ..... ❝ It will pay off to have a ‘better’ formulation. This is what one would expect. I least I think I had read something like that in SABE papers. [nitpicking] But I doubt that your chosen CVs give a pooled CV of 0.3. As a one-liner using the PowerTOST helper functions I got: mse2CV(mean(CV2mse(c(0.25,0.3458)))) .[/nitpicking] What did the one-liner do? Calculate the intra-subject variances for T and R, calculate the mean of both and back-calculate from that mean the CV ![]() Or did I miss somefink Welch here? — Regards, Detlew |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2013-02-15 14:01 (4455 d 12:53 ago) @ d_labes Posting: # 10032 Views: 10,478 |
|
Dear Detlew, ❝ As a one-liner using the PowerTOST helper functions I got: ❝ ❝ ❝ ❝ What did the one-liner do? Calculate the intra-subject variances for T and R, calculate the mean of both and back-calculate from that mean the CV Sure. ![]() ❝ Or did I miss somefink Welch here? Shame on me; was a typo – should read 0.3438 …mse2CV(mean(CV2mse(c(0.25,0.3438)))) — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |