test_that("Test fitAbn.bayes()", {

  load(file="testdata/fitabn_ex0.Rdata")
  # load(file='tests/testthat/testdata/fitabn_ex0.Rdata')

  myres.c.test <- fitAbn(dag=mydag, data.df=mydat, data.dists=mydists, method = "bayes")

  expect_equal(unclass(myres.c.test)[["mliknode"]][[1]], myres.c[[1]])
  expect_equal(unclass(myres.c.test)[["mliknode"]][[2]], myres.c[[2]])
  expect_equal(unclass(myres.c.test)[["mliknode"]][[3]], myres.c[[3]])
  expect_equal(unclass(myres.c.test)[["mliknode"]][[4]], myres.c[[4]])
  expect_equal(unclass(myres.c.test)[["mliknode"]][[5]], myres.c[[5]])
  expect_equal(unclass(myres.c.test)[["mliknode"]][[6]], myres.c[[6]])
  expect_equal(unclass(myres.c.test)[["mliknode"]][[7]], myres.c[[7]])
  expect_equal(unclass(myres.c.test)[["modes"]], myres.c[[8]])
  expect_equal(unclass(myres.c.test)[["error.code"]], myres.c[[9]])
  expect_equal(unclass(myres.c.test)[["hessian.accuracy"]], myres.c[[10]])
  expect_equal(unclass(myres.c.test)[["error.code.desc"]], myres.c[[11]])
  expect_equal(unclass(myres.c.test)[["mlik"]], myres.c[[12]])
  expect_equal(unname(unclass(myres.c.test)[["used.INLA"]]), myres.c[[13]]) # historical reasons. Can be updated in the future.
})