diff --git a/tests/testthat/test-DAISIE_ML2.R b/tests/testthat/test-DAISIE_ML2.R index 38635aa0..c310a0cc 100644 --- a/tests/testthat/test-DAISIE_ML2.R +++ b/tests/testthat/test-DAISIE_ML2.R @@ -19,38 +19,27 @@ test_that("use", { idparsopt = c(2, 4, 5, 6, 7), parsfix = c(0, Inf), idparsfix = c(1, 3), + methode = 'odeint::runge_kutta_cash_karp54', tol = c(0.01, 0.1, 0.001), res = 15, - tolint = c(0.1, 0.01), + tolint = c(1E-16, 1E-10), verbose = 0 ) expected_MLE <- data.frame( lambda_c = c(0.0, - 0.315015693879739, + 0.1198600880941345, 0.0, 0.0), - mu = c(2.67928476648249, - 2.67928476648249, - 2.67928476648249, - 2.67928476648249), - K = c(Inf, - Inf, - Inf, - Inf), - gamma = c(0.141748213992604, - 0.338044376412124, - 0.141748213992604, - 0.141748213992604), - lambda_a = c(1.10809210249609, - 1.10809210249609, - 1.10809210249609 , - 1.10809210249609), - loglik = c(-409.60117813893, - -409.60117813893, - -409.60117813893, - -409.60117813893), - df = c(5L, 5L, 5L, 5L), - conv = c(0L, 0L, 0L, 0L) + mu = rep(1.0604828698048443,4), + K = rep(Inf,4), + gamma = c(0.0544106971300572, + 0.1706160417441286, + 0.0544106971300572, + 0.0544106971300572), + lambda_a = rep(0.4749396123469573,4), + loglik = rep(-449.6525311168169310,4), + df = rep(5L,4), + conv = rep(0L,4) ) testthat::expect_equal(tested_MLE, expected_MLE) })