Bug? [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2020-11-29 17:22 (1268 d 17:52 ago) – Posting: # 22094
Views: 5,894

(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
23,034 posts in 4,835 threads, 1,647 registered users;
35 visitors (0 registered, 35 guests [including 6 identified bots]).
Forum time: 12:14 CEST (Europe/Vienna)

Give me a fruitful error any time, full of seeds, bursting with its own corrections.
You can keep your sterile truth for yourself.    Vilfredo Pareto

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