RE: Bayesian fits

From: Ruediger Port <rep@hsc.usc.edu> Subject: RE: [NMusers] Bayesian fits Date: Wed, 7 Nov 2001 10:38:18 -0800 (PST) Hi Nicolas, here as a $PK control file for estimating individual Bayes parameters when the population parameters are supposedly known. The principle is just to set MAXEVALS = 0 in $ESTIMATION (to prevent the population parameters from being re-estimated) and to request the POSTHOC step. The individual Bayes parameters can then be written to the output table by specifying them in $TABLE. Good luck! Ruedi (r.port@dkfz.de) ------------------------------------------------------------------------------ ; ka1.PK.control: linear one-compartment model, central elimination, ; peroral administration $PROBLEM first-order absorption $INPUT ID TIME EVID CMT PCMT CP=DV MG=AMT RATE $DATA data IGNORE=# $SUBROUTINES ADVAN2 TRANS2 $PK KA = THETA(1)*EXP(ETA(1)) CL = THETA(2)*EXP(ETA(2)) ; CL/F V = THETA(3)*EXP(ETA(3)) ; V/F S2 = V $ERROR ;" if (F.LE.0.) print *,"F",F ; verbatim print statement - sometimes ; useful when things are going wrong ... ; use it without $ESTIMATION otherwise ; report file will clog your hard disk mep = F ; "mixed-effects prediction" Y = F + F*EPS(1) $THETA 1.52 ; mean ka (1/h) 5.46 ; mean CL (L/h) 23.2 ; mean V (L) $OMEGA .66 .246 .32 ; omega^2: ka, CL, V $SIGMA .026 ; sigma^2 $ESTIMATION PRINT=1 MAXEVALS=0 POSTHOC $TABLE ID TIME PCMT EVID AMT KA CL V mep NOPRINT ONEHEADER FILE=anmtable