ubiquity/inst/ubinc/systems/system-one_cmt_cl.txt at master · john-harrold/ubiquity · GitHub

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
80
81
82
83
84
85
86
87
# Author: John Harrold <john.m.harrold@gmail.com>
#
# One compartment model with absorption compartment (At). The central
# compartment (Cp) had a volume (Vp). The system is parameterized in terms of
# macro cosntances (CL).
#
# While the the default dosing is a 1 mg IV dose into the central compartment,
# the model is written to accept dosing into an absorption compartment and
# through continuous IV infusion.
#
#       _________
#      |         |
#      |   At    |
#      |         |
#      |_________|
#           |
#           | ka, fb
#           |
#           V
#       _________
#      |         |
#      |   Cp    |
#      |   Vp    |
#      |_________|
#           |
#           | CL
#           |
#           V
#
# System Units:
#
#   mass          [=] mg
#   volume        [=] ml
#   concentration [=] mg/ml
#   time          [=] hr


# #-------------#
# | Parameters  |
# #-------------#
#
# System parameters
#    name              value     lower      upper   units editable    grouping
#                                bound      bound

<P>  Vp                  1.0       eps        Inf   ml         yes     System
<P>  CL                  1.0       eps        Inf   ml/hr      yes     System
<P>  ka                  1.0       eps        Inf    1/hr      yes     System
<P>  fb                  1.0       eps        Inf   --         yes     System


# #-------------------#
# |Input Information |
# #-------------------#
#
# Bolus Events
# ------------
# times/events state   values        scale      units
<B:times>;             [  0  ];      1;         hours
<B:events>;      Cp;   [1.0  ];      1/Vp;       mg
<B:events>;      At;   [0.0  ];      1;          mg

# Infusion Rate
# ------------
# #  name     time/levels  values  scale    units
<R:Dinf>;     times;       [0];     1;    hours
<R:Dinf>;     levels;      [0];     1;    mg/hour

# #-----------------------------#
# | ODEs, and State Information |
# #-----------------------------#

<ODE:At> -ka*At
<ODE:Cp>  ka*At*fb/Vp - CL/Vp*Cp + Dinf/Vp

# #---------#
# | Outputs |
# #---------#

<O> Cp_mg_ml    = Cp

# #---------#
# | Options #
# #---------#
# specify different time scales
<TS:hours> 1.0
<TS:days>  1.0/24