Lack of ISR / only 5% ISR in study with 656 samples [Regulatives / Guidelines]

posted by nobody – 2014-10-23 12:06 (3873 d 07:06 ago) – Posting: # 13775
Views: 11,132

❝ Next do this:

for (i in 1:sims) {

  xboot = sample(data, replace=T)

  y      <- length(xboot[abs(xboot)>20]) # No. of abs. deviations > 20%

  y      <- y/ISR.no        # Fraction of repeats with > 20% deviation

  if (y > 1/3) failed <- failed + 1

}

cat(sprintf("%i",sims),"simulated studies with", size, "samples:\n",

  "Method: Bootstrapping (Mersenne-Twister) of ISR data\n",

  "Failed criterion (deviations >20% in >1/3 of repeats) seen in",

  sprintf("%.2f%% %s", 100*failed, "of studies.\n"))


Wanted to perform some kind of sensitivity analysis to see at which criterion studies start to fail. With 1/8 of samples deviating more than 20% I get 200% study failure, due to simply multiplying number of failed studies by 100 in the output, but should be divided by number of studies simulated, I guess. The code above works only as long as 0 trials fail, i guess ;-)

Οὐδείς


Edit: Full quote removed. Please delete everything from the text of the original poster which is not necessary in understanding your answer; see also this post. [Helmut]
Of course, you are correct. The line should be:
  sprintf("%.2f%% %s", 100*failed/sims, "of studies.\n"))

Kindest regards, nobody

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,681 registered users;
41 visitors (0 registered, 41 guests [including 31 identified bots]).
Forum time: 19:13 CEST (Europe/Vienna)

Nerds don’t just happen to dress informally.
They do it too consistently.
Consciously or not, they dress informally
as a prophylactic measure against stupidity.    Paul Graham

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