test_that("LRMed2 runs and returns expected structure", {
  data("data.sim.med")
  medModel <- LRMed2(data = data.sim.med,
                     exposure = "pesticide1",
                     mediator = "hormone",
                     outcome = "ySim",
                     confounderVec = c("age", "invwt", "race2", "race3", "race4", "race5",
                                       "season2", "season3", "season4", "smoking1", "ovum1",
                                       "diabetes1"),
                     mValue = 0.15, eValueLow = 0.1, eValueHigh = 1.1)
  expect_true(!any(is.na(medModel$CDE)))
  expect_true(!any(is.na(medModel$NDE)))
  expect_true(!any(is.na(medModel$NIE)))
})