Bug? [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2020-11-29 17:22 (1238 d 22:49 ago) – Posting: # 22094
Views: 5,714

(edited by d_labes on 2020-11-29 17:57)

Dear Helmut,

❝ PS: Do we have a bug in scABEL()?

CV  <- c(0.3, 0.3+10^(-6:-2))

lim <- scABEL(CV)

res <- data.frame(CV = CV, L = lim[, 1], U = lim[, 2])

print(res, row.names = FALSE)


      CV         L        U

0.300000 0.8000000 1.250000

0.300001 0.8000296 1.249954

0.300010 0.8000244 1.249962

0.300100 0.7999731 1.250042

0.301000 0.7994604 1.250844

0.310000 0.7943616 1.258873


Not a bug :no:, it's a feature :cool:!
Reason: r_const=0.76.
With r_const <- log(1.25)/CV2se(0.3) = 0.7601283... we get
       CV         L        U
 0.300000 0.8000000 1.250000
 0.300001 0.7999994 1.250001
 0.300010 0.7999943 1.250009
 0.300100 0.7999430 1.250089
 0.301000 0.7994302 1.250891
 0.310000 0.7943307 1.258922

Satisfied?

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
22,988 posts in 4,825 threads, 1,649 registered users;
29 visitors (0 registered, 29 guests [including 10 identified bots]).
Forum time: 17:11 CEST (Europe/Vienna)

As soon as we abandon our own reason, and are content
to rely upon authority, there is no end to our troubles.    Bertrand Russell

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