bootf2BCA v1.2 [Software]

posted by Helmut Homepage – Vienna, Austria, 2019-10-05 00:59 (419 d 04:12 ago) – Posting: # 20665
Views: 3,489

Dear all,

we discussed the code at BioBridges…

The problem with reproducibility can easily be solved.
In the script F2_boot.R after line 270 #main algorithm add a new line


where X is an integer (specified in the protocol and report). Use

c(-1L, +1L) * .Machine$integer.max

to check the possible range. On my machine it gives

[1] -2147483647  2147483647

However, I suggest a simple one (42?).

» Be warned: The source code is hard to read, nearly no comments at all.

The current version 1.2 is a little bit better. Still heavy stuff.

The cut-off preferred by the EMA is implemented:


Don’t get confused by Q>=85%. That’s typo. Actually it is the correct >85%.
Still missing: One value >85% for the test (FDA) or one value >85% for the comparator (WHO).

» Be further warned: If you rise the number of bootstrap samples (default is 1000, too low IMHO) the result file will contain a huge number of entries, each bootstrap sample is recorded! The same mess as in the 'Object Pascal' implementation. But could of course also changed with a little knowledge of R.

Hint: In F2_boot.R comment out the loop in lines 317–329. Then the report-file shrinks from a couple of megabytes to some kilobytes. I would keep the rest for documentation. Useful to read the bootstrapped ƒ2 values from the file and generate your own histogram (density instead of counts). I don’t like the one of [image] plotly. Way too many bins for my taste (here an example with 5,000 bootstraps). \(\small{\textrm{ceiling}(\sqrt n)}\) gives 71 bins and smells of Excel.


I prefer the [image] Freed­man–Dia­conis rule* and show additionally the median and 90% CI (green line if the lower CL ≥50).


Dif-tor heh smusma 🖖
Helmut Schütz

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

 Admin contact
21,213 posts in 4,426 threads, 1,483 registered users;
online 15 (0 registered, 15 guests [including 11 identified bots]).
Forum time: Friday 04:12 CET (Europe/Vienna)

Biostatistician. One who has neither the intellect for mathematics
nor the commitment for medicine but likes to dabble in both.    Stephen Senn

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