test_that("FERS runs and returns expected structure", {
  data("data.sim.fixed")
  shape <- FERS(y = "ySim",
                xMain = "hormone",
                xConf = c("age", "invwt", "race2", "race3", "race4", "race5",
                          "season2", "season3", "season4", "smoking1", "ovum1", "diabetes1"),
                dataset = data.sim.fixed,
                nBasis = 5,
                nIter = 50)
  expect_true(!any(is.na(shape$PValueIncr)))
  expect_true(!any(is.na(shape$PValueDecr)))
  expect_true(!any(is.na(shape$PValueConv)))
  expect_true(!any(is.na(shape$PValueConc)))
})