bootf2bca v1.2 [Software]

posted by Helmut Homepage – Vienna, Austria, 2019-10-05 00:59  – Posting: # 20665
Views: 639

Dear all,

we discussed the code at BioBridges…

The problem with reproducibility can easily be solved.
In F2_boot.R after line 270 #main algorithm add e.g., set.seed(123456).

» 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:

[image]

Don’t get confused by Q>=85%. Actually it is the correct >85%. Still missing: one value >85% for the test (FDA) and one value >85% for the reference (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 plotly. Way too many bins for my taste (here an example with 5,000 bootstraps).

[image]


I prefer [image] Freed­man–Dia­conis* and show additionally the median and 90% CI.

[image]




Cheers,
Helmut Schütz
[image]

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

Complete thread:

Activity
 Admin contact
20,140 posts in 4,247 threads, 1,386 registered users;
online 7 (0 registered, 7 guests [including 3 identified bots]).
Forum time (Europe/Vienna): 04:35 CET

The scientific spirit is of more value than its products,
and irrationally held truths may be more harmful
than reasoned errors.    Thomas Henry Huxley

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