-
Notifications
You must be signed in to change notification settings - Fork 0
/
study_data.R
56 lines (51 loc) · 2.08 KB
/
study_data.R
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
# Sum of exponential data for sourcing
# -----------------------------------------------------------------------------
# The datasets in order are:
# d1b - one compartment kinetics given iv
# d2b - two compartment kinetics given iv
# d3b - three compartment kinetics given iv
# d1a - one compartment kinetics given orally
# d2a - two compartment kinetics given orally
# d3a - three compartment kinetics given orally
# -----------------------------------------------------------------------------
time.samp <- seq(0, 24, by = 0.05)
# Two Compartment Kinetics
pred.d2b <- function(x, p) {
exp(p[1]*x + p[3]) + exp(p[2]*x + p[4])
}
d2b.p <- c(-0.06, -0.5, 5, 6)
d2b <- pred.d2b(time.samp, d2b.p)
#with(d2b, plot(time, log(conc)))
d2b.t <- c(0, 0.5, 1, 3, 5, 8, 12, 16, 24)
# Three Compartment Kinetics
pred.d3b <- function(x, p) {
exp(p[1]*x + p[4]) + exp(p[2]*x + p[5]) + exp(p[3]*x + p[6])
}
d3b.p <- c(-0.03, -0.15, -0.5, 3.8, 4.7, 6)
d3b <- pred.d3b(time.samp, d3b.p)
#with(d3b, plot(time, log(conc)))
d3b.t <- c(0, 0.5, 1, 2, 4, 8, 12, 16, 24)
# One Compartment Kinetics w/ Absorption
pred.d1a <- function(x, p) {
exp(p[1]*x + p[3]) - exp(p[2]*x + p[3])
}
d1a.p <- c(-0.2, -0.4, 4)
d1a <- pred.d1a(time.samp, d1a.p)
#with(d1a plot(time, log(conc)))
d1a.t <- c(0, 0.5, 1, 3, 5, 7, 10, 16, 24)
# Two Compartment Kinetics w/ Absorption
pred.d2a <- function(x, p) {
exp(p[1]*x + p[4]) + exp(p[2]*x + p[5]) - exp(p[3]*x + log(sum(exp(p[4]), exp(p[5]))))
}
d2a.p <- c(-0.06, -0.4, -0.6, log(exp(4)*0.15), log(exp(4)*0.85))
d2a <- pred.d2a(time.samp, d2a.p)
# with(d2a, plot(time, log(conc)))
d2a.t <- c(0, 0.5, 1, 3, 5, 8, 12, 16, 24)
# Three Compartment Kinetics w/ Absorption
pred.d3a <- function(x, p) {
exp(p[1]*x + p[5]) + exp(p[2]*x + p[6]) + exp(p[3]*x + p[7]) - exp(p[4]*x + log(sum(exp(p[5]), exp(p[6]), exp(p[7]))))
}
d3a.p <- c(-0.05, -0.25, -0.6, -1, log(exp(4)*0.10), log(exp(4)*0.35), log(exp(4)*0.55))
d3a <- pred.d3a(time.samp, d3a.p)
#with(d3a, plot(time[-1], log(conc)[-1]))
d3a.t <- c(0, 0.5, 1.5, 3, 5, 8, 12, 16, 24)