BE-proff ● 2016-08-30 10:35 (3129 d 03:09 ago) Posting: # 16612 Views: 7,044 |
|
Hi All, How significant is a use of upper/low limits for calculation of sample size? Will I have a serious error using only CV without limits? ![]() |
ElMaestro ★★★ Denmark, 2016-08-30 11:04 (3129 d 02:39 ago) @ BE-proff Posting: # 16613 Views: 5,893 |
|
Hi BE-proff, ❝ How significant is a use of upper/low limits for calculation of sample size? The use of confidence limits on the CV not in particularly wide practical use, but it is advocated by some. If you go with an upper limit (as compared to the CV itself), that is going to be conservative, sample size will increase, power will go up. If you go with the lower limit, then you may be risking a too low sample size, the risk of failure increases. If you use limit then you need a level like 80% or 90% or whatever. It is really not set in stone and guidelines do not give an answer. ❝ Will I have a serious error using only CV without limits? Error is a tricky word, and it has a negative feel doesn't. You will always have some uncertainty on the CV before you plan the trial. Whether you want to translate this uncertainty into a confidence limit is of course up to you. If I were you and having some uncertainty I would do the following: a. Try and generate a curve of power versus N for some value of CV and GMR. b. Now adjust the GMR so that it is a little farther from 1.0 than the GRM under a. and generate the curve. c. Now try and increase CV a generate the curve. You will quickly be able to derive a few important, but mostly non-quantitative relationships, from such experiements. It will give you some ideas about the "what if"-scenario that you are alluding to, one way or another. At the end of the day all this stuff about power, confidence limits on GMR, confidence limits on CV, type I errors etc only works if the residual in the linear model is normal. While it is mandatory to use that model (normal residual) for evaluation and assessment, usually/often/sometimes the residual is really not normal. We never test for it, and that in itself translates directly into additioanal uncertainty on whatever estimate you are working on, regardless of how fancy your formulae otherwise are. — Pass or fail! ElMaestro |
BE-proff ● 2016-08-30 15:53 (3128 d 21:51 ago) @ ElMaestro Posting: # 16614 Views: 5,793 |
|
Hi ElMaestro, Thank you for clarification. Do you have an example of code for PowerTOST to plot a curve for power assessment? |
ElMaestro ★★★ Denmark, 2016-08-30 16:07 (3128 d 21:36 ago) @ BE-proff Posting: # 16615 Views: 5,790 |
|
Hi BE-proff, ❝ Do you have an example of code for PowerTOST to plot a curve for power assessment? Something like this for a 222BE study: GetMeTheCurve=function(CV, GMR) — Pass or fail! ElMaestro |
mittyri ★★ Russia, 2016-08-30 16:25 (3128 d 21:19 ago) @ ElMaestro Posting: # 16616 Views: 5,864 |
|
Hi ElMaestro and BE-proff, or even simpler (thanks to Detlew and Helmut!): library("PowerTOST") ![]() — Kind regards, Mittyri |
d_labes ★★★ Berlin, Germany, 2016-08-30 18:39 (3128 d 19:04 ago) @ mittyri Posting: # 16617 Views: 5,745 |
|
Dear mittyri, ❝ or even simpler (thanks to Detlew and Helmut!): The Academy award of the best actor goes to ... Helmut ![]() ![]() BTW: @All Don't forget the note under 'Details' of the man page of pa.ABE() and similar functions pa.scABE() or pa.NTIDFDA() :"It should be kept in mind that this is not a substitute for the "Sensitivity Analysis" recommended in ICH-E9. In a real study a combination of all effects occurs simultaneously. It is up to you to decide on reasonable combinations and analyze their respective power." — Regards, Detlew |