tests/testthat/test-enzyme-rate.R

In BradyAJohnston/biochemr: Tools For User-Friendly & Reproducible Analysis of Biochemistry Experiments

testthat::test_that("bio_enyme_rate is not equal to drc::drm()", {
  values1 <- drc::drm(rate ~ conc,
    data = Puromycin,
    curveid = state,
    fct = drc::MM.2(names = c("Vmax", "Km"))
  ) %>%
    coefficients() %>%
    unname() %>%
    sort() %>%
    round(4)

  values2 <- Puromycin %>%
    dplyr::group_by(state) %>%
    bio_enzyme_rate(conc, rate) %>%
    bio_coefs() %>%
    dplyr::pull(estimate) %>%
    sort() %>%
    round(4)

  testthat::expect_equal(values1, values2)
})