[NMusers] help with indirect reponse model

  Dear all:

  With the following code and data, the predicted fit is poor and no
standard error output for parameter. I highly appreciate your kind
help. Thanks a lot in your busy time.

  $PROB PLK INDIRECT RESPONSE MODEL
$INPUT ID TIME DV AMT CMT
$DATA PLKDOSE10SIMUL.CSV IGNORE=C
$SUBROUTINE ADVAN6 TOL=5

  $MODEL
COMP=(CENTRAL,DEFDOSE)
COMP=PERIPH
COMP=EFFECT

  $PK
K10=THETA(1)*EXP(ETA(1))
V1=THETA(2)*EXP(ETA(2))
K12=THETA(3)*EXP(ETA(3))
K21=THETA(4)*EXP(ETA(4))
SYNTH=THETA(5)*EXP(ETA(5))
LOSS=SYNTH
IC50=THETA(6)*EXP(ETA(6))
S1=V1

  $DES
DADT(1)=K21*A(2)-(K12+K10)*A(1)
DADT(2)=K12*A(1)-K21*A(2)
INH=IC50/(IC50+A(1)/S1)
DADT(3)=SYNTH-LOSS*INH*A(3)

  $ERROR
CP=A(1)/S1
Q=1
IF (CMT .EQ. 2) Q=0
Y=Q*F*EXP(ERR(1))+(1-Q)*F*EXP(ERR(2))
IPRE=Y

  $THETA
(0,1.3)   ;K10
(0,3)     ;S1
(0,0.3)   ;K12
(0,0.37)  ;K21
(0,1)     ;SYNTH
(0,10)    ;IC50

$OMEGA
0 FIX
0 FIX
0 FIX
0 FIX
0 FIX
0 FIX

  $SIGMA
0.04
0.5

  $ESTIMATION MAXEVAL=9999 POSTHOC PRINT=5
$COV
$TABLE NOPRINT FILE=PLKPKPD.txt ONEHEADER ID TIME AMT CMT CP IPRE

        CID TIME DV AMT CMT
        1 0 0 16807 1
        1 0 0 1 3
        1 0 5 0 3
        1 0.0833 4403 0 1
        1 0.0833 5.82 0 3
        1 0.5 2007 0 1
        1 0.5 5.23 0 3
        1 1 948 0 1
        1 1 8.1 0 3
        1 2 323 0 1
        1 2 8.99 0 3
        1 8 27 0 1
        1 8 15.86 0 3
        1 16 3 0 1
        1 16 33.39 0 3
        1 24 33.75 0 3
        1 48 24.19 0 3
        1 72 10.43 0 3
        1 96 8.03 0 3

lu

Received on Tue Nov 27 2007 - 23:43:03 EST