CECIF
☆    

Colombia,
2017-11-02 23:45
(2337 d 20:22 ago)

Posting: # 17964
Views: 5,262
 

 Vd/F calculation in bear [🇷 for BE/BA]

Hello everybody

I have been using bear for bioequivalence calculation since version 6.2. I consider it very useful and approachable and I am really grateful with yjlee168 for it.

For regulations in my country I have to make my own validation of bear's estimates. therefore I have been using a study of my own and I comapre my result with the ones I got from bear with the same study. Most of the values came out quite well, on exception of some of the AUC0inf and AUC0t mostly on the last sampling times which came out a little higher in my calculations (I can only suposse it is because of the decimal numbers).

However, the one parameter I have been struggling to get right in any way was Vd/F. Therefore I really need to ask what method bear uses to calculate it. I have been trying different equations to get it, but I failed in every one.
is it Vd/F= D*Ka/(I(Ka-Ke))? with ka from the residuals method or the wagner-nelson method?
or is it simply Vd/F=D/(AUC0inf*Ke)?
Vd/F= Apparent volume of distribution
D=dosage
I=intercept of the concentration vs time graph (x=0) when extrapolating the excretion line
ke=elimination constant (lambda z)
ka=absortion constant

Feel free to correct any missconception from my equations too.

I would really appreciate any light you could shed on this matter.
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2017-11-03 00:46
(2337 d 19:20 ago)

@ CECIF
Posting: # 17965
Views: 4,532
 

 Vd/F calculation in bear

Hi CECIF,

❝ […] the one parameter I have been struggling to get right in any way was Vd/F. Therefore I really need to ask what method bear uses to calculate it. […]


❝ or is it simply Vd/F=D/(AUC0inf*Ke)?


Exactly. The wonderful thing with open source is that you can perform a “whitebox” validation. With commercial software you can only perform a “blackbox” validation (rubbish in, rubbish out).

bear’s source code is available at
https://sourceforge.net/projects/yjlee-r-packages/files/bear/src/bear_2.8.3.tar.gz/download

Inspecting the script NCA.r you find in line 284
  VdFRef[j]<-Dose/(aucINF*ke)
and in line 520
  VdFTest[j]<-Dose/(aucINF*ke1)

Hope that helps.


PS: The same formula is used in Phoenix/WinNonlin.
PPS: IMHO, the “volume of distribution” after an extravascular dose (if no i.v. dose was administered in the same study) bears [pun!] limited – if any – meaning. Since F is unknown, what’s the purpose of reporting V/F?

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
yjlee168
★★★
avatar
Homepage
Kaohsiung, Taiwan,
2017-11-03 00:59
(2337 d 19:08 ago)

@ CECIF
Posting: # 17966
Views: 4,462
 

 Vd/F calculation in bear

Hi CECIF,

❝ I have been using bear for bioequivalence calculation since version 6.2.


uhh... since version 6.2? bear? the current version is 2.8.3 (as shown in my signature panel).

❝ ... on exception of some of the AUC0inf and AUC0t mostly on the last sampling times which came out a little higher in my calculations (I can only suposse it is because of the decimal numbers).


What was your selection for trapezoidal method? linear-up/log-down (default method) or all linear? That could make a little differences.

❝ ... I have been trying different equations to get it, but I failed in every one.

❝ is it Vd/F= D*Ka/(I(Ka-Ke))? with ka from the residuals method or the wagner-nelson method?

❝ or is it simply Vd/F=D/(AUC0inf*Ke)?


Simply Vd/F=D/(AUC0inf*Ke). Again please check your selected trapezoidal method. It also affects Vd/F too as you can see. What about Cl/F?

Hope this message can help.

All the best,
-- Yung-jin Lee
bear v2.9.1:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here
CECIF
☆    

Colombia,
2017-11-04 00:54
(2336 d 19:12 ago)

@ yjlee168
Posting: # 17968
Views: 4,343
 

 Vd/F calculation in bear

❝ uhh... since version 6.2? bear? the current version is 2.8.3 (as shown in my signature panel).


You are right, I meant 2.6.2. Sorry for the brain fart. :-(


❝ What was your selection for trapezoidal method? linear-up/log-down (default method) or all linear? That could make a little differences.


My selection was default


❝ [...] What about Cl/F?



Now That I looked for teh code as Helmut told me is all very clear. It seems however that both Vd/F and Cl/F are affected by the same misscalculation from AUC0inf. I will check with the other methods and tell you later.

Anyway, your answers were very helpful. Many thanks
UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,636 registered users;
86 visitors (0 registered, 86 guests [including 7 identified bots]).
Forum time: 20:07 CET (Europe/Vienna)

Nothing shows a lack of mathematical education more
than an overly precise calculation.    Carl Friedrich Gauß

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