NONMEM Users Network Archive

Hosted by Cognigen

Re: multiple problems in a control file

From: Nick Holford <n.holford>
Date: Fri, 02 Oct 2009 17:09:10 +1300

Pierre-Olivier,

If you want to fit each subject in a file as a single subject then you
can use something like the following pair of control streams. The method
relies on the use of INCLUDE to create a multiple problem control
stream. Note that the second and subsequent problems cannot have $PK or
$ERROR records. The second and subsequent problems are automatically
included in the primary control stream with the INCLUDE statement. The
number of times this is done is determined by the number on the INCLUDE
line. Typically this will be the number of subjects in the data file
minus one.

Nick

========= Run1 - run this primary control stream ==========
$PROB Single sub
$INPUT ID AMT TIME DV MDV
$DATA data10.csv RECS=ID NOREWIND ; 10 subjects
$EST METHOD=ZERO
MAXEVAL=9999 SIG=3 NOABORT
$COV

$THETA
(0,0.5,);K
(0,3000,);V
(0, 1.3);KA
$OMEGA
0.1 ; eps

$SUBR ADVAN2 TRANS1
$PK
K=THETA(1)
V=THETA(2)
KA=THETA(3)
S2=V
$ERROR
Y = F + F*ERR(1)

INCLUDE run2.ctl 9 ; Replace number with N subjects minus 1


====== The following is a secondary control stream with the name run2.ctl
====== The data file and this control stream must be located in the same
directory
====== that NM-TRAN is running in or include a path from the NM-TRAN run
directory to
====== the control stream in the primary control stream and to the data
file in this
====== second control stream

$PROB Single sub
$INPUT ID AMT TIME DV MDV
$DATA data10.csv RECS=ID NOREWIND 10 subjects
$EST METHOD=ZERO
MAXEVAL=9999 SIG=3 NOABORT
$COV

$THETA
(0,0.5,);K
(0,3000,);V
(0, 1.3);KA
$OMEGA
0.1 ; eps




Tremblay, Pierre-Olivier wrote:
>
> Hi,
>
> I'm trying to run multiple single-subject problems within a single
> control file but there is something I'm not doing right because when I
> include the $SUBROUTINE LIBRARY option, then NONMEM won't estimate the
> model parameters. However if I omit that option, I can't run the
> second problem. I pasted my test control file below .
>
> $PROB RUN# 001
> $INPUT C ID SUBJ AMT RATE EVID DAY DAT1=DROP TAD TIME DV MDV
> $DATA TESTDATA.CSV RECS=ID IGNORE=C LAST20=50
> $SUBROUTINES ADVAN2 TRANS1 SUBROUTINE=LIBRARY
> $PK
> K=THETA(1)
> V=THETA(2)
> KA=THETA(3)
> S2=V
> $ERROR
> Y = F + F*ERR(1)
> IPRED=F
> $EST PRINT=5 MAX=9999 SIG=3
> $THETA
> (0,0.5,);[K]
> (0,3000,);[V]
> (0, 1.3);[KA]
> $OMEGA
> 0.1 ;[P]
> $PROB RUN# 002
> $INPUT C ID SUBJ AMT RATE EVID DAY DAT1=DROP TAD TIME DV MDV
> $DATA TESTDATA.CSV RECS=ID NOREWIND IGNORE=C LAST20=50
> $PK
> K=THETA(1)
> V=THETA(2)
> KA=THETA(3)
> S2=V
> $ERROR
> Y = F + F*ERR(1)
> IPRED=F
> $EST PRINT=5 MAX=9999 SIG=3
> $THETA
> (0,0.5,);[K]
> (0,3000,);[V]
> (0, 1.3);[KA]
> $OMEGA
> 0.1 ;[P]
> I'd really appreciate if somenone could point me in the right
> direction with this.
>
> Tx :o)
>
> Pierre-Olivier
>
>
> -----------------------------------------
> Ce courriel et toutes ses pièces jointes sont confidentiels et ne
> peuvent être utilisés que par le destinataire ou l'organisation à qui
> ils sont adressés. Si vous n'êtes pas le destinataire visé, veuillez
> en informer l'expéditeur et supprimez immédiatement le courriel. Ce
> message confirme également que ce courriel et ses pièces jointes ont
> été inspectés par notre logiciel d'antivirus pour votre sécurité.
> -----------------------------------------
> This e-mail and any files transmitted with it are confidential and
> intended solely for the use of the individual or entity to whom they
> are addressed. If you have received this email in error please notify
> the sender and delete it immediately. This footnote also confirms that
> this e-mail message and its attached file(s) have been swept for your
> security by our antivirus software using the latest virus signatures.
> Anapharm Par Proofpoint

--
Nick Holford, Professor Clinical Pharmacology
Dept Pharmacology & Clinical Pharmacology
University of Auckland, 85 Park Rd, Private Bag 92019, Auckland, New Zealand
n.holford
mobile: +64 21 46 23 53
http://www.fmhs.auckland.ac.nz/sms/pharmacology/holford
Received on Fri Oct 02 2009 - 00:09:10 EDT

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.