Re: Bayesian fits
From: Pierre Maitre <maitre@cdg.ch>
Subject: Re: [NMusers] Bayesian fits
Date: Mon, 05 Nov 2001 12:25:28 +0100
Paul Hutson a écrit :
>
> However, it is not clear to me how to code the control stream when
previous
> data regarding the estimates of THETA and ETA are known,
Long time ago, I used to do this with a PRED routine that I wrote
myself, and I must say that I have never tried with NMTRAN.
However, I believe that the following NMTRAN code should work:
$PROB xxxxxxxxxn 2cpt sc and iv
$DATA ../data/xxxxxx_PK.nm NOREWIND
$INPUT ID TIME SESS CMT TTMT GEND AGE HEIG WEIG DOSE=AMT RATE CP=DV
AS NEOP B2MI MDV EVID
$SUBROUTINES ADVAN4 TRANS4
$PK
CALLFL=-2
;
CL = THETA(1)*EXP(ETA(1))
V2 = THETA(2)
Q = THETA(3)*EXP(ETA(2))
V3 = THETA(4)
KA = THETA(5)*EXP(ETA(3))
S2 = V2
$ERROR
IPRED=F
IRES=CP-IPRED
IWRES=IRES/IPRED
Y = F*EXP(EPS(1))
;" WRITE(6,*) Y
$THETA 4.65 FIXED
10.4 FIXED
17.8 FIXED
287 FIXED
0.359 FIXED
$OMEGA 0.301 FIXED 0.0673 FIXED 1.68 FIXED
$SIGMA 0.199
$ESTIM PRINT=3, MAXEVAL=9999, POSTHOC
$COVARIANCE OMITTED
$TABLE UNCONDITIONAL NOPRINT ONEHEADER FILE= sdtab
ID TIME SESS CMT TTMT DOSE RATE CP IPRED IRES IWRES
$TABLE UNCONDITIONAL NOPRINT ONEHEADER FILE= indiv_param
ID CL V2 Q V3 KA
$SCATTER RES VS PRED
$SCATTER PRED VS CP UNIT
$SCATTER WRES VS PRED
$SCATTER RES VS TIME
$SCATTER RES VS ID
$SCATTER WRES VS ID
Best regards
Pierre Maitre
--
Dr Pierre-O. Maitre
Privat Docent
FMH Anesthésiologie
FMH Pharmacologie Clinique
Cabinet médical / A la Joy
CH-1272 Genolier Switzerland