bootf2bca v1.2 [Software]

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

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]




Dif-tor heh smusma 🖖
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,810 posts in 4,355 threads, 1,447 registered users;
online 11 (0 registered, 11 guests [including 10 identified bots]).
Forum time: 18:51 CEST (Europe/Vienna)

Medical statistician: One who will not accept that Columbus discovered America…
because he said he was looking for India in the trial plan.    Stephen Senn

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