+--------------------------------------------------------------------+
 |                                                                    |
 |                          BAYES EXAMPLE 5                           |
 |                                                                    |
 +--------------------------------------------------------------------+

 This is example5.ctl from the NONMEM 7 distribution medium.  It, along
 with the data file, can be found in the examples directory.

 ;Model Desc: Population Mixture Problem in 1 Compartment model,
 ; with rate constant parameter mean modeled for two
 ; sub-populations, but its inter-subject variance is the same in
 ; both sub-populations Project Name: nm7examples
 ;Project ID: NO PROJECT DESCRIPTION

 $PROB RUN# example5 (from ad1tr1m4t)
 $INPUT C SET ID JID TIME CONC=DV DOSE=AMT RATE EVID MDV CMT VC1
        K101 VC2 K102 SIGZ PROB
 $DATA example5.csv IGNORE=C

 $SUBROUTINES ADVAN1 TRANS1

 $MIX
 P(1)=THETA(4)
 P(2)=1.0-THETA(4)
 NSPOP=2

 $PK
 Q=1
 IF(MIXNUM.EQ.2) Q=0
 MU_1=THETA(1)

 ; Note that MU_2 can be modeled as THETA(2) or THETA(3),
 ; depending on the MIXNUM value.
 ; Also, we are avoiding IF/THEN blocks.

 MU_2=Q*THETA(2)+(1.0-Q)*THETA(3)
 V=DEXP(MU_1+ETA(1))
 K=DEXP(MU_2+ETA(2))
 S1=V

 $ERROR
 Y = F + F*EPS(1)

 $THETA
 (-1000.0  4.3 1000.0)  ;[MU_1]
 (-1000.0 -2.9 1000.0)  ;[MU_2-1]
 (-1000.0 -0.67 1000.0) ;[MU_2-2]
 (0.0001 0.667 0.9999)  ;[P(1)]

 $OMEGA BLOCK(2)
 0.04 ;[p]
 0.01  ;[f]
 0.04 ;[p]

 $SIGMA
 0.01 ;[p]

 $EST METHOD=ITS INTERACTION NITER=100 PRINT=1 NOABORT SIGL=8
      FILE=example5.ext CTYPE=3

 $EST METHOD=IMPMAP INTERACTION NITER=20 ISAMPLE=300 PRINT=1
      NOABORT SIGL=8

 $EST METHOD=IMP INTERACTION NITER=20 ISAMPLE=1000 PRINT=1
      NOABORT SIGL=6

 $EST NBURN=500 NITER=500 METHOD=SAEM INTERACTION PRINT=10
      SIGL=6 ISAMPLE=2

 $EST METHOD=IMP INTERACTION NITER=5 ISAMPLE=1000 PRINT=1
      NOABORT SIGL=6 EONLY=1 MAPITER=0

 $EST METHOD=BAYES INTERACTION NBURN=2000 NITER=5000 PRINT=10
      FILE=example5.txt SIGL=8

 $EST MAXEVAL=9999 NSIG=2 SIGL=8 PRINT=10 FILE=example5.ext
      METHOD=CONDITIONAL INTERACTION NOABORT

 $COV MATRIX=R

REFERENCES: Guide Introduction_7

  
Go to main index.
  

    Created by nmhelp2html v. 1.0 written by Niclas Jonsson (Modified by AJB 5/2006,11/2007,10/2012)