test_that("FERSInt runs and returns expected structure", {
  data("data.sim.fixed.int")
  shape <- FERSInt(y = "ySim",
                   xExp = "pesticide1",
                   xMed = "hormone",
                   xConf = c("pesticide1", "age", "invwt", "race2", "race3", "race4", "race5",
                             "season2", "season3", "season4", "smoking1", "ovum1", "diabetes1"),
                   dataset = data.sim.fixed.int,
                   nBasis = 5,
                   nIter = 1000)
  expect_true(!any(is.na(shape[[1]]$PValueIncr)))
  expect_true(!any(is.na(shape[[1]]$PValueDecr)))
  expect_true(!any(is.na(shape[[1]]$PValueConv)))
  expect_true(!any(is.na(shape[[1]]$PValueConc)))
  expect_true(!any(is.na(shape[[2]]$PValueIncr)))
  expect_true(!any(is.na(shape[[2]]$PValueDecr)))
  expect_true(!any(is.na(shape[[2]]$PValueConv)))
  expect_true(!any(is.na(shape[[2]]$PValueConc)))
})