The CMT data item works differently with DOSE and OTHER event records (EVID=1 and 2) than it does with Observation event records (EVID=0). You need to use the PCMT data item with EVID=2. Look for PCMT in the help or html directories.
On Fri, 22 Apr 2011 11:29 +0200, "Németh István (FGYFK)" <[email protected]> wrote: Dear Users, After a succesful estimation of PK data by using ADVAN6 I try to get some more detailed predictions for the central and peripherial compartment as well. Therefore some "blank" records were inserted to the initial data file, like this: ID AMT EVID TIME DV UVOL CMT MDV 1 4e+06 1 0.0 . 0 1 1 1 0e+00 0 0.0 0.00000 0 2 0 1 0e+00 0 0.0 0.00000 0 4 0 1 0e+00 2 0.0 . 0 2 1 1 0e+00 2 0.0 0.00000 0 3 0 1 0e+00 2 0.0 . 0 4 1 1 0e+00 0 0.5 13.93836 0 2 0 1 0e+00 0 0.5 44.12349 0 4 0 1 0e+00 0 1.0 15.06895 0 2 0 1 0e+00 0 1.0 101.76711 0 4 0 1 0e+00 2 1.0 . 0 2 1 1 0e+00 2 1.0 . 0 4 1 1 0e+00 0 2.0 13.14943 0 2 0 1 0e+00 0 2.0 108.56616 0 4 0 1 0e+00 2 2.0 . 0 2 1 1 0e+00 2 2.0 . 0 4 1 1 0e+00 2 3.0 . 0 2 1 1 0e+00 2 3.0 . 0 4 1 1 0e+00 0 4.0 10.92107 0 2 0 1 0e+00 0 4.0 117.83321 0 4 0 CMT=2 is the central CMT=4 is a peripherial comparment. The predictions for the central compartment were OK, but the predictions for peripherial compartment were identical with the central. On the top of this, there were no errors with those records where the DV was a measured value. Is there any idea? Thanks, István The model: $PROBLEM $INPUT C ID AMT EVID TIME DV UVOL CMT $DATA DATANONMEM.CSV IGNORE=C $SUBROUTINE ADVAN6 TRANS1 TOL=4 $MODEL NCOMP=4 COMP=(DEPOT,DEFDOS) COMP=(CENTRAL,DEFOBS) COMP=(URINE) COMP=(PERIPH) $PK TVKA=THETA(1) TVCL=THETA(2) TVCLR=THETA(3) TVV2=THETA(4) TVKCS=THETA(5) TVKSC=THETA(6) TVV4=THETA(7) KA=TVKA*EXP(ETA(1)) CL=TVCL*EXP(ETA(2)) CLR=TVCLR*EXP(ETA(3)) V2=TVV2*EXP(ETA(4)) KCS=TVKCS*EXP(ETA(5)) KSC=TVKSC*EXP(ETA(6)) V4=TVV4*EXP(ETA(7)) S2=V2 S3=UVOL S4=V4 $DES DADT(1)=-KA*A(1) DADT(2)= KSC*A(4)+KA*A(1)-A(2)*(CLR+CL+KCS) DADT(3)= CLR*A(2) DADT(4)= KCS*A(2)-KSC*A(4) $THETA (0,0.1) (0,0.1) (0,0.1) (0,160) (0,0.1) (0,0.1) (0,10) $OMEGA .1 .1 .1 .1 .1 .1 .1 $SIGMA 0.1 10 $ERROR IPRE=F Y=IPRE*(1+ERR(1))+ERR(2) $ESTIMATION PRINT=5 METHOD=1 POSTHOC $TABLE C ID TIME CMT EVID DV IPRE NOPRINT ONEHEADER FILE=benzopred.tab István Németh PhD Biostatistician Developmental Drug Metabolism & Pharmacokinetics Pharmacological and Drug Safety Research Gedeon Richter Plc Tel: 00 36 1 505 7082 Fax: 00 36 1 889 8526 -- Alison Boeckmann [email protected]
