Cmax, Tmax, and Time and AUC above threshold such as MIC
Title: Paul R
I would like to identify maximal concentrations in tissue,
corresponding time of Cmax tissue, and time and AUC of these
concentrations above a target concentration (eg, MIC for antibiotics or
1 uM for taxol). I identified the Cmax and associated Tmax in
simulated data for 'tissue' in the following code adapted from earlier
posts 99jul292004.html
. It works, but in the output table I am then required to
extract the Cmax and Tmax as the maximum values in the $TABLE (eg MAX()
in EXCEL, accomplishing little with the NONMEM code), or manually as
the last value in columns 'Cmax' and 'Tmax' before the null values
begin.
So to my questions:
1) How can the singular Cmax value identified by the last non-null
value in my IF statement be captured (apparently the "read-out" that
Prof Sheiner referred to in his comments on Cmax in earlier threads)
with its corresponding Tmax value and saved, for example, in a file
such as 'patab1'?
2) In a somewhat related fashion, has anyone code they can suggest to:
a) Start and stop the elapsed time for which concentrations are
above an MIC or other specified concentration (or effect)
b) Store previous time and concentration values so that the AUC
during this time above MIC can be calculated using trapezoid or
log-trapezoid methods?
Many thanks, as always.
Paul
$SUBROUTINES ADVAN6 TOL=7
$MODEL
COMP=(CENTRAL DEFDOSE DEFOBS)
COMP=(TISSUE)
$PK
CL=THETA(1)*WGT*EXP(ETA(1)); CL
V1=THETA(2)*WGT*EXP(ETA(2)); VPARENT
Q= THETA(3)*WGT; Q
V2=THETA(4)*WGT; V2
K10=CL/V1
K12=Q/V1
K21=Q/V2
S1=V1
$DES
D1=-A(1)*(K10+K12)+A(2)*K21
DADT(1)=D1
D2=A(1)*K12 - A(2)*K21
DADT(2)=D2
IF (D2.GT.0) THEN
CMax= A(2)/V2
TMax= TIME
ELSE
CMax=0
TMax=0
ENDIF
$TABLE ID TIME AMT RATE BOLT WGT CNC1 CNC2 Y1 Y2 CMaX TMaX NOPRINT
ONEHEADER
--
Paul R.
Hutson, Pharm.D.
Associate
Professor
UW School
of Pharmacy
777
Highland Avenue
Madison
WI 53705-2222
Tel 608.263.2496
Fax
608.265.5421
Pager
608.265.7000, p7856