Bug? [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2020-11-29 16:22 (47 d 14:11 ago) – Posting: # 22094
Views: 1,232

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

Activity
 Admin contact
21,303 posts in 4,441 threads, 1,488 registered users;
online 6 (0 registered, 6 guests [including 6 identified bots]).
Forum time: Saturday 06:33 CET (Europe/Vienna)

If you can’t solve a problem, then there is an easier problem
you can solve: find it.    George Pólya

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