Use PowerTOST [Power / Sample Size]

posted by d_labes  – Berlin, Germany, 2017-12-06 10:52 (2011 d 04:42 ago) – Posting: # 18037
Views: 13,474

Dear Astea!

❝ ...

❝ By the way, for extreme numbers from the paper by Endrényi L, Tóthfalusi L. we also get slightly underpowered study:

power.scABEL(CV=0.8,n=54,theta0=0.95,design="2x3x3")

❝ [1] 0.79917

power.scABEL(CV=0.8,n=72,theta0=0.9,design="2x3x3")

❝ [1] 0.79782

❝ ...


The sample sizes in the paper of the two Laszlos are based on 10 000 simulations (only).
Additionally the seed of the pseudo-random number generator for the sims was presumably different.
Try this:
power.scABEL(CV=0.8, n=54, theta0=0.95, design="2x3x3", nsims=10000)
[1] 0.7976
power.scABEL(CV=0.8, n=54, theta0=0.95, design="2x3x3", nsims=10000, setseed=F)
[1] 0.7995
power.scABEL(CV=0.8, n=54, theta0=0.95, design="2x3x3", nsims=10000, setseed=F)
[1] 0.8029


With nsims=1E6 and setseed=F on the other hand I don't succeed to obtain power>0.8 (20 tries).
One more reason to use PowerTOST's power.scABEL() / sampleN.scABEL() functions :cool:. Also w.r.t the tables in the paper of the two Laszlos.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
22,627 posts in 4,743 threads, 1,613 registered users;
26 visitors (0 registered, 26 guests [including 13 identified bots]).
Forum time: 16:34 CEST (Europe/Vienna)

Facts are meaningless. You could use facts to prove
anything that’s even remotely true!    Homer Simpson

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