test_that("sd_interpret_estimates() returns the expected data frame", {

  estimates <- c(par_beta = 0,
                 par_rho  = 0.8472979,
                 I0       = 0,
                 inv_phi  = -2.3025855)

  par_list <- list(list(par_name  = "par_beta",
                        par_trans = "exp"),
                   list(par_name  = "par_rho",
                        par_trans = "expit"),
                   list(par_name  = "I0",
                          par_trans = "exp"),
                   list(par_name  = "phi",
                        par_trans = c("exp", "inv")))

 actual   <- sd_interpret_estimates(estimates, par_list)

 expected <- data.frame(par_beta = 1, par_rho = 0.7, I0 = 1, phi = 10)

 expect_equal(actual, expected, tol = 1e-4)


 estimates <- data.frame(c(0, 0.6931472),
                         c(0.8472979, 0.8472979),
                         c(0, 0.4054651),
                         c(-2.3025855, -2.3025855))

 actual   <- sd_interpret_estimates(estimates, par_list)

 expected <- data.frame(par_beta = c(1, 2),
                        par_rho = 0.7,
                        I0 = c(1, 1.5),
                        phi = 10)

 expect_equal(actual, expected, tol = 1e-4)
})