Hello This is a code + data example for an IV infusion, PARENT = 1 = cpt1 METAB = 2 = cpt2
Regards Saik Saïk URIEN, MD, PhD Directeur de Recherche, INSERM Service de Pharmacologie Centre René Huguenin 35 rue Dailly F92210 St Cloud Tel. 01 4711 1684 Fax 01 4711 1617 email : [EMAIL PROTECTED] $PROB R IV INFUS $INPUT ID TIME AMT RATE DV CMT MDV $DATA DATA125.csv $SUBROUTINE ADVAN6 TOL=3 $MODEL COMP=(PARENT) COMP=(METAB) $PK V1=THETA(1)*EXP(ETA(1)) CL=THETA(2)*EXP(ETA(2)) P12=THETA(3)*EXP(ETA(3)) K20=THETA(4)*EXP(ETA(4)) S1=V KEL =CL/V K12 =P12*KEL K10 =KEL - K12 V2 = 1 $DES DADT(1) = -(K10+K12)*A(1) DADT(2) = K12*A(1) - K20*A(2) $ERROR C1=A(1)/V1 C2=A(2)/V2 IPRED = C1 IF (CMT.EQ.2) IPRED=C2 Y=C1*(1+EPS(1))+EPS(2) IF (CMT.EQ.2) Y=C2*(1+EPS(3))+EPS(4) $THETA (0,5);V $THETA (0,10);CL $THETA (0,.005,1);f12 $THETA (0,.01);K20 $OMEGA 0.4;oV $OMEGA 0.4;oCL $OMEGA 0.1;of12 $OMEGA 0.1;oK20 $SIGMA .4;sPp $SIGMA 1;sPa $SIGMA .4;sMp $SIGMA .1;sMa $EST METHOD=1 INTERACTION MAXEVAL=4500 PRINT=10 POSTHOC $COV $TABLE ID TIME CMT NOPRINT FILE=r DATA CORRESPONDING TO THE CODE ABOVE #ID TIME AMT RATE DV CMT MDV 1101 0 433.36 216.68 0 1 1 1101 1 0 0 11.84 1 0 1101 1 0 0 0.44 2 0 1101 2 0 0 41.17 1 0 1101 2 0 0 1.05 2 0 1101 2.08 0 0 21.08 1 0 1101 2.08 0 0 0.91 2 0 1101 2.25 0 0 11.78 1 0 1101 2.25 0 0 0.82 2 0 1101 2.5 0 0 6.38 1 0 1101 2.5 0 0 0.6 2 0 1101 3 0 0 1.48 1 0 1101 3 0 0 0.26 2 0 1101 4 0 0 0.17 1 0 1101 4 0 0 0.06 2 0 1101 6 0 0 0.02 1 0 1102 0 1132.65 566.32 0 1 1 1102 1 0 0 33.49 1 0 1102 1 0 0 1.17 2 0 1102 2 0 0 38.02 1 0 1102 2 0 0 1.65 2 0 1102 2.08 0 0 45.7 1 0 1102 2.08 0 0 1.77 2 0 1102 2.25 0 0 46.49 1 0 1102 2.25 0 0 3.22 2 0 1102 2.5 0 0 40.18 1 0 1102 2.5 0 0 1.82 2 0 1102 3 0 0 13.41 1 0 1102 3 0 0 1.22 2 0 1102 4 0 0 2.7 1 0 1102 4 0 0 0.46 2 0 1102 6 0 0 0.12 1 0 1102 6 0 0 0.09 2 0 -----Message d'origine----- De : [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] la part de Vijay V. Upreti Envoyé : lundi 12 février 2007 19:20 À : [email protected] Objet : [NMusers] Simultaneous drug and metabolite POPPK Dear NMusers, I have plasma concentration versus time data for a drug and its active metabolite. I would like to simulataneously fit this drug and its metabolite using NONMEM to obtain estimates of prameters and associated variability for both drug and its active metabolite. I am sure guys have done this nicely in past. I was wondering if someone can share a control stream and sample data for my guidance. Thanks _______________ Vijay V Upreti PhD Candidate Pharmacokinetics-Biopharmaceutics Laboratory Dept. of Pharmaceutical Sciences University of Maryland, School of Pharmacy Rm S509, 20 Penn St., Baltimore, MD 21201 Voice: 410.706.7388; Fax: 410.706.5017 ********************************************************************************************** IMPORTANT: The contents of this email and any attachments are confidential. They are intended for the named recipient(s) only. If you have received this email in error, please notify the system manager or the sender immediately and do not disclose the contents to anyone or make copies thereof. *** eSafe scanned this email for viruses, vandals, and malicious content. *** **********************************************************************************************
