RE: Model estimates based on 1 sample per subject
From: tgordi@buffalo.edu
Subject: RE: [NMusers] Model estimates based on 1 sample per subject
Date: 8/19/2003 6:55 PM
I apologize for not having provided the control stream. It is now given
at the end of this mail. As you might guess, the V1 and CL values are
arbitrary. Beside this model, I also tried one where DADT(1) was set
equal to -A(1)*CL/V1. That model gives almost identical result as the
one presented here.
Thank you!
Toufigh
$PROB CELL
$INPUT ID AMT TIME DV REG CMT
$DATA ind.csv IGNORE=#
$SUBROUTINE ADVAN6 TRANS1 TOL=5
$MODEL NCOMP=1
COMP=CENT
$PK
V1=1
CL=0.00001
S1=V1
IMAX=THETA(1)*EXP(ETA(1))
KM=THETA(2)
KNS =THETA(3)
$DES
DADT(1)=0
$ERROR
CONC=F
LAM=(IMAX*CONC/(KM+CONC))+KNS*CONC
Y=LAM*(1+EPS(1))
$THETA (0, 500) ;1 IMAX
$THETA (0, 10000) ;2 KM
$THETA (0, 0.001) ;3 KNS
$OMEGA 0.1 ;1 IIV IMAX
$SIGMA 0.1 ;1 PROP ERROR
$ESTIMATION NOABORT POSTHOC MAXEVAL=9999 PRINT=3 METHOD=1 MSFO=MSF1
$COVARIANCE
$TABLE ID TIME Y NOPRINT ONEHEADER FILE=SDTAB1
$TABLE ID TIME IMAX KM KNS ETA1 NOPRINT ONEHEADER FILE=PATAB1
$TABLE ID REG NOPRINT ONEHEADER FILE=CATAB1