test_that("Simple model with cvll runs", {

  set.seed(123456)
  b0 <- .2 # True value for the intercept
  b1 <- .5 # True value for the slope
  n <- 500 # Sample size
  X <- runif(n, -1, 1)

  Y <- b0 + b1 * X + rnorm(n, 0, 1) # N(0, 1 error)

  obj_cvll_ols <- cvll(Y ~ X, data.frame(cbind(Y, X)), method = "OLS")
  obj_cvll_mr <- cvll(Y ~ X, data.frame(cbind(Y, X)), method = "MR")

  obj_cvlldiff <- cvlldiff(obj_cvll_ols$cvll, obj_cvll_mr$cvll,
                           obj_cvll_ols$df)

  test_stat <- obj_cvlldiff$test_stat
  check_against <- c(3.45354)
  expect_equal(round(as.numeric(test_stat), 5),
               check_against)

  test_pval <- obj_cvlldiff$p_value
  check_against <- c(0.00030)
  expect_equal(round(as.numeric(test_pval), 6),
               check_against)

})