RE: Query on model development

From: Atul Bhattaram Venkatesh Date: July 09, 2007 technical Source: mail-archive.com
Hello Tausif Please check the scaling factor that you included in the model. It should be S2=V/1000. I think this should correct the differences that you are observing. Atul Venkatesh Atul Bhattaram Pharmacometrics US Food and Drug Administration "The contents of this message are mine personally and do not necessarily reflect any position of the Government or the Food and Drug Administration." ________________________________
Quoted reply history
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Tausif Ahmed Sent: Monday, July 09, 2007 7:27 AM To: [email protected] Subject: [NMusers] Query on model development Dear All, I am new to the field of POP-PK. We are developing a POP-PK model for our in-house drug. The conc. are in ng/mL and dose is in mg. The drug is administered orally and follows 1-compartment model. It has a short half-life of 2-3 h. We have developed the basic model using rich data from serial sampling of healthy subjects. -We have determined the Bayesian estimates of CL and V. But I think, we are overestimating these values. If we determine CL , using non-compartmental analysis (rich data )we are getting lower estimates (800 L/h), while using NONMEM (using same rich data), we are getting quite high estimates (mean for 25 mg- 4000 L/hr). Kindly give your inputs on further modifying this model to get an irreducible model and also validating the model. Thanking you in anticipation. Thanks and regards Tausif Ahmed, Ph.D. Senior Research Scientist Metabolism and Pharmacokinetics Department, Ranbaxy Research Lab., India Model: PROB Oral RSd-RMD-PK01-healthy study $DATA RSD_PK01ctdata.csv IGNORE=C $INPUT ID TIME CONC=DV DOSE AMT AGE WT MDV $SUBROUTINES ADVAN2 TRANS2; One Compartment Linear Model $PK TVCL=THETA(1) ; mean clearance CL=TVCL*EXP(ETA(1)) TVV=THETA(2) ; mean central volume V=TVV*EXP(ETA(2)) K=CL/V ;reparameterization required S1=V/1000 KA=THETA(3) $THETA (0, 100) ; clearance estimate (L\hr) (0, 100) ; volume estimate (L) (0, 1.5) ; Absorption rate constant $OMEGA .16 .16 ; forty percent cv $SIGMA .16 ; forty percent cv $ERROR IPRED=F Y=F+ERR(1) $ESTIMATION METHOD=0 MAXEVAL=9999 PRINT=5 POSTHOC SIGDIGITS=4 $COVARIANCE $TABLE ID TIME DV DOSE CL V ETA1 ETA2 AMT AGE WT MDV NOPRINT ONEHEADER FILE=CS1_oralRSDRMDPK01.PAR $SCAT CONC VS TIME OUTPUT MINIMUM VALUE OF OBJECTIVE FUNCTION ************************************************** 6971.242 ************************************************** ************************************************************************ ************************************************ ******************** ******************** FINAL PARAMETER ESTIMATE ************************************************************************ ************************************************ THETA - VECTOR OF FIXED EFFECTS PARAMETERS ********* TH 1 TH 2 TH 3 4.97E+03 8.02E+03 1.70E-01 OMEGA - COV MATRIX FOR RANDOM EFFECTS - ETAS ******** ETA1 ETA2 ETA1 + 2.57E+00 ETA2 + 0.00E+00 1.14E+00 SIGMA - COV MATRIX FOR RANDOM EFFECTS - EPSILONS **** EPS1 EPS1 + 4.72E+02 (i) The information contained in this e-mail message is intended only for the confidential use of the recipient(s) named above. This message is privileged and confidential. If the reader of this message is not the intended recipient or an agent responsible for delivering it to the intended recipient, you are hereby notified that you have received this document in error and that any review, dissemination, distribution, or copying of this message is strictly prohibited. If you have received this communication in error, please notify us immediately by e-mail, and delete the original message. (ii) The sender confirms that Ranbaxy shall not be responsible if this email message is used for any indecent, unsolicited or illegal purposes, which are in violation of any existing laws and the same shall solely be the responsibility of the sender and that Ranbaxy shall at all times be indemnified of any civil and/ or criminal liabilities or consequences there. <<Public.gif>>
Jul 09, 2007 Tausif Ahmed Query on model development
Jul 09, 2007 Atul Bhattaram Venkatesh RE: Query on model development
Jul 09, 2007 Leonid Gibiansky Re: Query on model development
Jul 09, 2007 Krishna Devarakonda Re: Query on model development
Jul 09, 2007 Nitin Mehrotra Re: Query on model development