1
2
3
4
5
6
7
8
9
10
11
12
13
|
library(pfm)
context("Actual and time weighted average concentrations for SFO kinetics")
test_that("SFO_actual_twa calculates correctly", {
test_times <- c(0, 1, 7, 21, 42)
# This was calculated with the CRD spreadsheet for multiple applications
reference <- data.frame(
actual = c(10, 9.330, 6.156, 2.333, 0.544),
twa = c(NaN, 9.661, 7.923, 5.267, 3.248),
row.names = test_times)
result <- round(10 * SFO_actual_twa(10, times = test_times), 3)
expect_equal(result, reference)
})
|
