From 4f62f39b8979cb3a5c52d65d425e530dc0770665 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Sun, 16 Feb 2025 17:27:18 +0100 Subject: Fix solution type "deSolve" for a special case The special case occurs if not observed data for time zero are available. Fixes #16. Write NEWS and update docs. A lot of html was regenerated, because I had checked out the main branch in between, which made pkgdown believe that I have changed the help files. --- tests/testthat/test_deSolve.R | 19 +++++++++++++++++++ 1 file changed, 19 insertions(+) create mode 100644 tests/testthat/test_deSolve.R (limited to 'tests/testthat') diff --git a/tests/testthat/test_deSolve.R b/tests/testthat/test_deSolve.R new file mode 100644 index 00000000..3d15de35 --- /dev/null +++ b/tests/testthat/test_deSolve.R @@ -0,0 +1,19 @@ +context("Solutions with deSolve") + +test_that("Solutions with deSolve work if we have no observations at time zero", { + skip_on_cran() + + # For testing purposes, we replace 0 values in the time column with 0.1 + # This confused mkinfit with solution_type "deSolve" up to version 0.1.3.0 + FOCUS_D_nozero <- FOCUS_D + FOCUS_D_nozero[FOCUS_D$time == 0, "time"] <- c(0.1, 0.1) + + f_sfo_sfo_nozero <- mkinfit(SFO_SFO, FOCUS_D_nozero, quiet = TRUE) + f_sfo_sfo_nozero_deSolve <- mkinfit(SFO_SFO, FOCUS_D_nozero, + solution_type = "deSolve", quiet = TRUE) + expect_equal( + parms(f_sfo_sfo_nozero), + parms(f_sfo_sfo_nozero_deSolve) + ) +}) + -- cgit v1.2.1