Bug? [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2020-11-29 16:22 (295 d 23:17 ago) – Posting: # 22094
Views: 2,816

(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,691 posts in 4,534 threads, 1,541 registered users;
online 15 (0 registered, 15 guests [including 4 identified bots]).
Forum time: Tuesday 16:40 CEST (Europe/Vienna)

“Data! Data! Data!” he cried impatiently.
“I can’t make bricks without clay!”    Arthur Conan Doyle

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