RE: Problem (bug ?) with RATE and F1 in $PK

From: "Piotrovskij, Vladimir [JanBe]" <VPIOTROV@janbe.jnj.com> Subject: RE: [NMusers] Problem (bug ?) with RATE and F1 in $PK Date: Mon, 5 Nov 2001 11:03:17 +0100 Pierre, What you described is indeed a NONMEM feature, but not a bug. In your case of SC infusion you know just duration, not the rate (it depends on the bioavailability). Below is the way how to code your case. I reduced the model and dropped ETAs for simplicity; duration is 5; MODE=3D1 corresponds to SC and MODE=3D2 to IV: $PROB sc and iv $DATA nmd.ssc $INPUT ID TIME CMT DOSE=3DAMT RATE CP=3DDV EVID MODE $SUBROUTINES ADVAN2 TRANS2 $PK CL =3D THETA(1) V =3D THETA(2) KA =3D THETA(3) D1 =3D 5 F1 =3D THETA(4) S2 =3D V $ERROR =20 Y =3D F * (1 + ERR(1)) $THETA (0 20) (0 100) (0 5) (0 .5) $OMEGA .01 $SIM (34906)=20 $EST PRINT=3D10 MAX=3D9990=20 $TABLE ID TIME CMT DOSE RATE EVID MODE FILE=3Dtab.ssc ONEHEADER NOPRINT An example (single-subject) data set is as follows: 1 0 1 100 -2 0 4 1 1 1 2 0 0 1 0 1 1 2 2 0 0 1 0 1 1 3 2 0 0 1 0 1 1 4 2 0 0 1 0 1 1 5 2 0 0 1 0 1 1 6 2 0 0 1 0 1 1 7 2 0 0 1 0 1 1 8 2 0 0 1 0 1 1 9 2 0 0 1 0 1 1 10 2 0 0 1 0 1 1 0 2 100 0 0 4 2 1 1 2 0 0 1 0 2 1 2 2 0 0 1 0 2 1 3 2 0 0 1 0 2 1 4 2 0 0 1 0 2 1 5 2 0 0 1 0 2 1 6 2 0 0 1 0 2 1 7 2 0 0 1 0 2 1 8 2 0 0 1 0 2 1 9 2 0 0 1 0 2 1 10 2 0 0 1 0 2 Best regards, Vladimir