(no subject)

5 messages 5 people Latest: Mar 05, 2011

(no subject)

From: Susan Al Taher Date: March 05, 2011 technical
Hi I am trying to build a structural model for an oral drug. I am getting extremely large and unreasonable values for vd, eg a billion L, putting an upper limit on theta will give a near boundary warning. Does anyone have an idea why this might occur and how to solve it. This is part of the input file: $SUBROUTINE ADVAN2 TRANS2 $PK TVCL = THETA(1) CL = TVCL*EXP(ETA(1)) V = THETA(2) KA = THETA(3) S2 = V $THETA (0,10) $THETA (0,20) $THETA 2 FIXED Best regards Susan

Re:

From: Shankar Lanke Date: March 05, 2011 technical
Dear Susan I am a beginner in NONMEM here is my input from one of my online resources. I appreciate if you others comment on this. Go through individual profiles one by one, perform NCA analysis in WINONLIN for initial estimates (for further information go through Chapter 10 of ENE ETTE BOOK) 1.Try different starting estimates. Sometimes NONMEM is very sensitive to the starting estimates of the parameters. For variability parameters, NONMEM tends to perform better if the initial estimates are too high, rather than too low. 2. Reduce the number of THETAs in the model. If a model is over-specified, then NONMEM can thrash trying to find the correct value for THETA. Reduce a 3 compartment model to a 2 compartment model. Eliminate covariate effects. The fewer parameters NONMEM tries to estimate, the more likely it is to successfully estimate the value of the parameters. 3. Reduce the number of variance parameters. Variance parameters can be particularly troublesome. Try having only 1 ETA and 1 EPS in the model. 4. Fix the values of ETA at 0. This is the naïve pooled data (NPD) approach, and it works surprisingly well. NONMEM sometimes has real problems estimating the structural model and the variance model at the same time. When this is the case, you may be able to get NONMEM to run by fixing every ETA value at 0 (e.g., $OMEGA 0 FIX 0 FIX 0 FIX ) and so on. With the ETA estimates fixed at 0, estimate the values of THETA. Then, fix the THETA values at those estimated, and now estimate the OMEGA matrix (i.e. ETA) for the fixed values of THETA. This way NONMEM isn’t trying to estimate both THETA and OMEGA at the same time. For further detailed information you should go through NONMEM manual 5 and this chapter helps you well in NONMEM 5 (Chapter 9 - Additional NM-TRAN Records) All the best.
Quoted reply history
On Sat, Mar 5, 2011 at 12:02 PM, Susan Al taher <[email protected]>wrote: > Hi > > I am trying to build a structural model for an oral drug. I am getting > extremely large and unreasonable values for vd, eg a billion L, putting an > upper limit on theta will give a near boundary warning. > > Does anyone have an idea why this might occur and how to solve it. > > This is part of the input file: > $SUBROUTINE ADVAN2 TRANS2 > $PK > TVCL = THETA(1) > CL = TVCL*EXP(ETA(1)) > V = THETA(2) > KA = THETA(3) > S2 = V > > > $THETA (0,10) > $THETA (0,20) > $THETA 2 FIXED > > > Best regards > Susan > -- Regards, Shankar Lanke Ph.D. University at Buffalo Office # 716-645-4853 Fax # 716-645-2886 Cell # 678-232-3567

Re:

From: Joseph . Kahn Date: March 05, 2011 technical
1) Likely your data units are off by a few orders of magnitude. (maybe concentrations are shown in ng/mL or maybe in ug/mL or pg/mL) 2) perhaps your bioavailability is also small 3) perhaps your drug units are in ug instead of mg? 4) Try scaling S2=V/1000000 to compensate 5)Don't put upperbounds on your volume till you know the right scale and order of magnitude Good luck! - Joseph Joseph Kahn, Ph.D. Novartis Pharmaceuticals Corporation Novartis Pharmaceuticals Corporation One Health Plaza East Hanover, NJ 07936-1080 USA Phone: +1 862 7783939 Email : [email protected] Susan Al taher <[email protected]> Sent by: [email protected] 03/05/2011 12:34 PM To <[email protected]> cc Subject Hi I am trying to build a structural model for an oral drug. I am getting extremely large and unreasonable values for vd, eg a billion L, putting an upper limit on theta will give a near boundary warning. Does anyone have an idea why this might occur and how to solve it. This is part of the input file: $SUBROUTINE ADVAN2 TRANS2 $PK TVCL = THETA(1) CL = TVCL*EXP(ETA(1)) V = THETA(2) KA = THETA(3) S2 = V $THETA (0,10) $THETA (0,20) $THETA 2 FIXED Best regards Susan

Re:

From: Isingh716 Date: March 05, 2011 technical
It might be helpful if you could show some of the data records. I am guessing it has something to do with the units of conc and dose. You can also simulate PK with your initial guesses and overlay on your data to see if these guesses are even reasonable to begin with or not ? Best Jeet Sent from my iPhone
Quoted reply history
On Mar 5, 2011, at 12:02 PM, Susan Al taher <[email protected]> wrote: > Hi > > I am trying to build a structural model for an oral drug. I am getting > extremely large and unreasonable values for vd, eg a billion L, putting an > upper limit on theta will give a near boundary warning. > > Does anyone have an idea why this might occur and how to solve it. > > This is part of the input file: > $SUBROUTINE ADVAN2 TRANS2 > $PK > TVCL = THETA(1) > CL = TVCL*EXP(ETA(1)) > V = THETA(2) > KA = THETA(3) > S2 = V > > > $THETA (0,10) > $THETA (0,20) > $THETA 2 FIXED > > > Best regards > Susan

Re:

From: Nathan Teuscher Date: March 05, 2011 technical
Susan, In my experience very large estimates for volume of distribution are almost always a results of inaccuracies with units between concentration observations (DV), dosing amounts (AMT), and the scale parameter (S). I believe you can correct your issue by adjusting the S2=V code in your control stream. Please see my blog (blog.learnpkpd.com) for a complete description. The specific post on the scale parameter is located here: http://wp.me/p1fZzZ-4O Best of luck! Nathan _____________________________ Nathan Teuscher Founder PK/PD Associates [email protected] www.learnpkpd.com blog.learnpkpd.com www.youtube.com/learnpkpd www.twitter.com/learnpkpd
Quoted reply history
On Mar 5, 2011, at 11:02 AM, Susan Al taher wrote: > Hi > > I am trying to build a structural model for an oral drug. I am getting > extremely large and unreasonable values for vd, eg a billion L, putting an > upper limit on theta will give a near boundary warning. > > Does anyone have an idea why this might occur and how to solve it. > > This is part of the input file: > $SUBROUTINE ADVAN2 TRANS2 > $PK > TVCL = THETA(1) > CL = TVCL*EXP(ETA(1)) > V = THETA(2) > KA = THETA(3) > S2 = V > > > $THETA (0,10) > $THETA (0,20) > $THETA 2 FIXED > > > Best regards > Susan