tests/testthat/helper-equivalence.R

In jlme: Regression Modelling with 'GLM.jl' and 'MixedModels.jl' in 'Julia'

loadNamespace("broom")
loadNamespace("broom.mixed")

mod_deframe <- function(df, ignore_names = FALSE) {
  df <- tidy(df) # resolves to appropriate method for {base}, {lme4}, {jlme}
  df$estimate <- signif(df$estimate, 1)
  if (ignore_names) {
    df$estimate
  } else {
    vec <- stats::setNames(df$estimate, df$term)
    vec[order(names(vec))]
  }
}

expect_similar_models <- function(x, y, ignore_names = FALSE) {
  if (is.list(x) && length(x) == 2 && missing(y)) {
    y <- x[[2]]
    x <- x[[1]]
  }
  x_res <- mod_deframe(x, ignore_names)
  y_res <- mod_deframe(y, ignore_names)
  expect_identical(x_res, y_res)
}