Re: multiple problems in a control file
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
[email protected] tel:+64(9)923-6730 fax:+64(9)373-7090
mobile: +64 21 46 23 53
http://www.fmhs.auckland.ac.nz/sms/pharmacology/holford