Re:

From: Shankar Lanke Date: March 05, 2011 technical Source: mail-archive.com
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
Mar 05, 2011 Susan Al Taher (no subject)
Mar 05, 2011 Shankar Lanke Re:
Mar 05, 2011 Joseph . Kahn Re:
Mar 05, 2011 Isingh716 Re:
Mar 05, 2011 Nathan Teuscher Re: