GM ★ India, 2018-11-02 07:25 (1994 d 20:26 ago) Posting: # 19522 Views: 4,206 |
|
Dear Dlabes, I am trying to calculate the sample size for replicate design using R software with sampleN.RSABE. With keen interest to know about the functionality of sampleN.RSABE, I have searched in the web. I got some backend code of this function from Github (see link). I am a beginner of R software, but I tried to understand this. There are some intermediate functions like .sampleN0, .sampleN0.2, .power.RSABE are involved. I got the informationn of first two functions from the same github. but not the third one. As you are the auther of this function, am directly requesting you. Kindly explain or suggest any link about this function. Thank you in advance. — Best Regards, GM |
Helmut ★★★ Vienna, Austria, 2018-11-02 11:50 (1994 d 16:02 ago) @ GM Posting: # 19523 Views: 3,521 |
|
Hi GM, not Detlew but answering anyway. ❝ I am trying to calculate the sample size for replicate design using R software with sampleN.RSABE. So why are you not simply using it?
?sampleN.RSABE in the R-console / R-Studio’s console or consult the online manual.❝ I am a beginner of R software If you want to dive into the code, be aware of R’s steep learning curve. ❝ […] but I tried to understand this. OK, this is a laudable. ❝ There are some intermediate functions like .sampleN0, .sampleN0.2, .power.RSABE are involved. ❝ I got the informationn of first two functions from the same github. but not the third one. .power.RSABE is contained in power_RSABE2L_isc.R .— Dif-tor heh smusma 🖖🏼 Довге життя Україна! Helmut Schütz The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
GM ★ India, 2018-11-03 06:09 (1993 d 21:43 ago) @ Helmut Posting: # 19525 Views: 3,503 |
|
Hi Helmut, ❝ not Detlew but answering anyway. Thank you for quick reply at every time. ❝ So why are you not simply using it? ❝
❝ sampleN.RSABE(CV=0.4, theta0=0.9, targetpower=0.8, design="2x2x4") If you want further information about the arguments and their defaults, type?sampleN.RSABE in the R-console / R-Studio’s console or consult the online manual.As this is a replicated study, there are only few articles for calculating the sample size. There is no clear information about this. So that, I suppose to understand concept through R-program. I have few more doubts, regarding the terms CVswitch, r_const and pe_constr. I thought CVswitch=0.30 for both FDA and EMA, r_const=0.25 and 0.760 for FDA and EMA respectively and don't know about the pe_constr. Kindly correct me if I am wrong... Thank you in advance. — Best Regards, GM |
Helmut ★★★ Vienna, Austria, 2018-11-03 14:04 (1993 d 13:47 ago) @ GM Posting: # 19526 Views: 3,667 |
|
Hi GM, ❝ As this is a replicated study, there are only few articles for calculating the sample size. If reference-scaling (RSABE, ABEL) is concerned, do you know any article except the two Lászlós’?1 We need simulations because the methods are implicitly sequential. Like in any framework (following a decision scheme) an analytical solution for power does not exist.
However, the sample size estimation for ABE in replicate designs is straightforward (i.e., does not require simulations). Use sampleN.TOST() instead.❝ There is no clear information about this. So that, I suppose to understand concept through R-program. Tóthfalusi et. al.2 is a good starting point. Other references are given in the man-pages of PowerTOST .❝ I thought CVswitch=0.30 for both FDA and EMA, … Correct. ❝ … r_const=0.25 and 0.760 for FDA and EMA respectively … Not quite. The switching condition θs (aka regulatory constant) is based on the regulatory standardized variation σw0. For the FDA σw0 = 0.25 and for the EMA based on CVwR 30% as \({\sigma _{w0}} = \sqrt {\log \left( {{{0.30}^2} + 1} \right)} = 0.2935604 \ldots \) Then \({\theta _s} = \log (1.25)/\sigma _{w0} = 0.8925742 \ldots \) (FDA) and \(0.7601283 \ldots \) (EMA). Note that σw0 0.25 is explicitly given by the FDA and therefore, the regulatory constant in PowerTOST is used in full precision.* On the other hand, the EMA requires the rounded 0.760 (termed k in the GL) and an upper cap for scaling at CVwR 50% (EL 69.84–143.19%). Check the conditions:
❝ … and don't know about the pe_constr. It was introduced by all agencies following suggestions by Les Benet (see this post). It is “political” and leads to statistical troubles (essentially we are truncating the distribution and the entire concept is built on sand).
reg_const("HC") .
— Dif-tor heh smusma 🖖🏼 Довге життя Україна! Helmut Schütz The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |