Re: Simulation with uncertainty + estimation
Matts,
You should be cautious about implementing this example as you've parameterized
the model. There is an inconsistency between your model parameterization and
the multivariate-normal assumption made when you use the var-cov matrix of the
estimates to incorporate parameter uncertainty. See earlier discussions on this
topic:
http://www.cognigencorp.com/nonmem/current/2009-July/1816.html
Marc
Quoted reply history
On Apr 8, 2010, at 4:20 AM, Kågedal, Matts wrote:
>
> Hi all,
> A model for simulation with uncertainty was previously posted by Mats
> Karlsson in this forum.
> After a few bug fixes (with help from Mats) it now actually works. See below:
> Matts Kågedal
> AstraZeneca
>
> $PROB
> $DATA data1.csv IGNORE=@
> $INPUT ID DV TIME XVCL XVV AMT
>
> ;data1=
> ; 1 0 0 0 0 100
> ; 1 0 1 0 0 0
> ; 1 0 21 0 0 0
> ; 2 0 0 0 0 100
> ; etc
> $SIM (20090726 NEW) SUBPROBLEMS=30 ;
> $THETA ; estimates of THETA and OMEGA from previous run
> (0,1) ; POP_CL theta1
> (0,10) ; POP_V theta2
> $OMEGA
> 0.01 ; PPV_CL eta1
> 0.01 ; PPV_V eta2
> ;variance-covariance matrix of the THETA estimates from previous run
> $OMEGA BLOCK(2)
> 0.02 ; UNC_POP_CL eta3
> 0.0001 0.2 FIX; UNC_POP_V eta 4
>
> $SIGMA .001
> $SUB ADVAN1
> $PK
> CL=THETA(1)*EXP(ETA(1)) ;
> V=THETA(2)*EXP(ETA(2)) ;
> ;S1=V
> TCLE=THETA(1); estimated typical value
> TVE=THETA(2); estimated typical value
>
> IF (ICALL.EQ.4.AND.NEWIND.EQ.0) THEN ; do this just once per
> subproblem
> DO WHILE (ETA(3).LT.-.99.OR.ETA(4).LT.-9.9) CALL SIMETA(ETA) ;resample if
> "bad" values
> END DO
> TVCL=THETA(1)+ETA(3)
> TVV =THETA(2)+ETA(4)
>
> WRITE (53,*) TVCL, TVV ; Output simulated THETAs to file
> ENDIF
>
> IF (ICALL.EQ.4) THEN ; do this for all simulated values
> CL=TVCL*EXP(ETA(1)) ; with uncertainty for CL
> V =TVV*EXP(ETA(2)) ; with uncertainty for V
> ENDIF
>
> S1=V
> K=CL/V
>
> $ERROR
> Y =F*EXP(EPS(1))
> $INFN
> IF (NEWIND.EQ.0) THEN
> WRITE (52,*) THETA ; Output estimated THETAs to file
> ENDIF
> $EST MAX=9999 METH=1 INTER
>
>
>
>
> Confidentiality Notice: This message is private and may contain confidential
> and proprietary information. If you have received this message in error,
> please notify us and remove it from your system and note that you must not
> copy, distribute or take any action in reliance on it. Any unauthorized use
> or disclosure of the contents of this message is not permitted and may be
> unlawful.
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>