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

posted by Aceto81 – Belgium, 2008-07-15 12:07 (6131 d 17:38 ago) – Posting: # 2035
Views: 60,796

❝ ❝ 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,424 posts in 4,927 threads, 1,672 registered users;
27 visitors (0 registered, 27 guests [including 3 identified bots]).
Forum time: 05:46 CEST (Europe/Vienna)

The whole purpose of education is
to turn mirrors into windows.    Sydney J. Harris

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