-
Notifications
You must be signed in to change notification settings - Fork 0
/
run155.mod
80 lines (68 loc) · 2.03 KB
/
run155.mod
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
;; 1. Based on: run154
;; 2. Description: Good IV+IG, Three dose, CL, BIO, 7 CMTs
;; x1. Author: user
$PROBLEM METABOLISM OF ARCTIGENIN IN RATS AFTER IV AND IG - ADVAN6
$DATA Arctigenin_IG+IV_for_nonmem.csv IGNORE=@
$INPUT ID TIME DV DOSE=AMT MDV CMT EVID DL RTE
$SUBROUTINE ADVAN6 TRANS1 TOL=9
$MODEL
COMP = (LUMEN DEFDOSE) ; 1. intestine lumen
COMP = (AR_CENTR) ; 2. AR in central
COMP = (AG_CENTR) ; 3. AG in central
COMP = (AA_CENTR) ; 4. AA in Central
COMP = (AR_INTES) ; 5. AR in Intestine
COMP = (AG_INTES) ; 5. AG in Intestine
COMP = (AA_INTES) ; 5. AA in Intestine
$PK
KA = THETA(1)*EXP(ETA(1)) ; absorption rate constant of AR
K20= THETA(2)*EXP(ETA(2)) ; elimination rate constant of AG from Central
K30= THETA(3)*EXP(ETA(3)) ; elimination rate constant of AA from central comp
V2 = THETA(4)*EXP(ETA(4)) ; volume of central
S2 = V2
S3 = V2
S4 = V2
Ki= THETA(5)*EXP(ETA(5)) ; Metabolism rate constant in intestine
Kb= THETA(6)*EXP(ETA(6)) ; Metabolism rate constant in Blood
Kl= THETA(7)*EXP(ETA(7)) ; Metabolism rate constant in Liver
F1 = THETA(8)*EXP(ETA(8)) ; Bioavailability
K10= THETA(9)*EXP(ETA(9)) ; elimination rate constant of AR from Central
K52= THETA(10)*EXP(ETA(10))
$DES
DADT(1)=-KA*A(1)
DADT(2) = K52*A(5)-Kb*A(2)-Kl*A(2)-K10*A(2)
DADT(3) = Ki*A(6)+Kl*A(2)-K20*A(3)
DADT(4) = Kb*A(7)+Kb*A(2)-K30*A(4)
DADT(5) = KA*A(1)-Ki*A(5)-Kb*A(5)-K52*A(5)
DADT(6) = Ki*A(5)-Ki*A(6)
DADT(7) = Kb*A(5)-Kb*A(7)
$ERROR
Y = F+ERR(1)
$THETA
(0, 0.03) ; KA(1/min)
(0, 0.1) ; K20(1/min)
(0, 0.1) ; K30(1/min)
(0, 0.3) ; V2(L)
(0, 0.15) ; Ki(1/min)
(0, 0.15) ; Kb(1/min)
(0, 0.05) ; Kl(1/min)
(0, 0.01,1) ; F1
(0 FIX) ; K10
(0 FIX) ; K52
$OMEGA
0 FIX ; IIV KA
0 FIX ; IIV K20
0 FIX ; IIV K30
0 FIX ; IIV V2
0 FIX ; IIV Ki
0 FIX ; IIV Kb
0 FIX ; IIV Kl
0 FIX ; IIV BIO
0 FIX ; IIV K10
0 FIX ; IIV K52
$SIGMA
0.01 ; Residual error
$EST METHOD=0 MAXEVAL=999 NOABORT SIG=5 PRINT=1
$COV
; Xpose
$TABLE ID TIME DV MDV CMT DL RTE ONEHEADER NOPRINT FILE=sdtab155.tab
$TABLE ID Ka K20 K30 V2 ONEHEADER NOPRINT FILE=patab155.tab