tests/testthat/test-submerge.R

In daniel1noble/orchaRd: Orchard Plots and Prediction Intervals for Meta-analysis and Meta-regression

context("Checking sub_merge function..")

data(fish)
warm_dat <- fish
model <- metafor::rma.mv(yi = lnrr, V = lnrr_vi,random = list(~1 | group_ID, ~1 | es_ID),
                         mods = ~ experimental_design + trait.type + deg_dif + treat_end_days, method = "REML", test = "t", data = warm_dat, control=list(optimizer="optim", optmethod="Nelder-Mead"))

# Intercept only
overall <- orchaRd::mod_results(model, group = "group_ID", data = warm_dat)

# Moderators
#Trait type
across_trait <- orchaRd::mod_results(model, group = "group_ID", mod = "trait.type", data = warm_dat)

merge <- orchaRd::submerge(across_trait, overall)

testthat::test_that("Checking submerge function for fish dataset ..", {

  testthat::expect_equal(
    dim(merge$mod_table)[[1]], 5,
    info = "Check that merging of 'overall' (1) and 'across_trait' (4) for fish have dimensions pf 5...")

})