Validation of PhEq_bootstrap [Software]

posted by Shuanghe  – Spain, 2018-07-17 15:54  – Posting: # 19077
Views: 2,225

Hi Helmut!

» Now for the bad news. There is one thing stupid in the software when it comes to validation (line 535 of the source Conv_to_learn1.pas):

RandSeed:=random(random(1000000))+random(random(1000000));

That means that the seed of the (pseudo)random number generator is always a random number itself. Therefore, it is impossible to reproduce runs (try it: you will not get exactly my results).

What a coincidence! I was in a meeting recently with a regulatory agency and the issue of reproducibility came up since the agency couldn't reproduce the bootstrap f2 values we reported in the dossier. My first thought was the different seed number. I opened the source files one by one and searched words like random, seed, to get the line. :-D

» IMHO, that’s bad coding practice.

Indeed. End user should be able to specify the seed number to be able to reproduce the analysis. I always leave a variable for seed number in my SAS code for bootstrap of BE data.

All the best,
Shuanghe

Complete thread:

Activity
 Mix view
Bioequivalence and Bioavailability Forum |  Admin contact
19,817 posts in 4,200 threads, 1,361 registered users;
online 9 (0 registered, 9 guests [including 8 identified bots]).
Forum time (Europe/Vienna): 16:22 CEST

A big computer, a complex algorithm and a long time
does not equal science.    Robert Gentleman

The BIOEQUIVALENCE / BIOAVAILABILITY FORUM is hosted by
BEBAC Ing. Helmut Schütz
HTML5