bear v1.0.0 - a data analytical tool for ABE in R [🇷 for BE/BA]

posted by Aceto81 – Belgium, 2008-07-15 12:07 (6550 d 01:22 ago) – Posting: # 2035
Views: 89,893

❝ ❝ more sophisticated and/or robust approaches for estimation of the terminal elimination rate are for example discussed in this post.

❝ A very good pointer and references there. Many thanks.


Yung-jin,

first of all thanks for sharing the Bear package.
For the estimation of the terminal elimination rate and the number of points to pick, I already wrote a function in the past which uses the WinNonLin approach as in the thread mentioned here before.
Maybe it is of some help (with some adjustments for your code).
dat <- data.frame(time,conc)
  n_lambda=0
  r.adj=0
  for (i in (nrow(dat)-3):which.max(dat$conc)) {
    if (r.adj - round(summary(lm(log(conc)~time,dat[i:nrow(dat),]))$adj.r.squared,4) < (0.0001))  {
      n_lambda <- nrow(dat)-i+1
      r.adj <- summary(lm(log(conc)~time,dat[i:nrow(dat),]))$adj.r.squared
    }
  }

Best regards

Ace

--
Edit: I changed the tabs in the post to serial blanks; otherwise on some displays the output gets confusing.
See also this post for updated code (preventing Cmax/tmax to be included). [Helmut]
Thread locked

Complete thread:

UA Flag
Activity
 Admin contact
23,655 posts in 4,993 threads, 1,570 registered users;
130 visitors (0 registered, 130 guests [including 9 identified bots]).
Forum time: 13:29 CEST (Europe/Vienna)

Science is simply common sense at its best that is,
rigidly accurate in observation, and
merciless to fallacy in logic.    Thomas Henry Huxley

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