test_that("fqardl runs for three quantiles", {
  set.seed(1)
  TT <- 60
  y  <- cumsum(rnorm(TT))
  x  <- cumsum(rnorm(TT))
  res <- fqardl(y = y, x = matrix(x, ncol = 1),
                tau = c(0.25, 0.5, 0.75),
                pmax = 2L, qmax = 2L, maxk = 1)
  expect_s3_class(res, "fqardl")
  expect_equal(ncol(res$beta), 3L)
  expect_length(res$rho, 3L)
})

test_that("fqardl no_fourier option works", {
  set.seed(2)
  TT <- 50
  y  <- cumsum(rnorm(TT))
  x  <- cumsum(rnorm(TT))
  res <- fqardl(y = y, x = matrix(x, ncol = 1),
                tau = 0.5, pmax = 2L, qmax = 2L,
                no_fourier = TRUE)
  expect_equal(res$kstar, 0)
})

test_that("fqardl long-run beta has correct dimensions", {
  set.seed(3)
  TT <- 60
  y  <- cumsum(rnorm(TT))
  x1 <- cumsum(rnorm(TT))
  x2 <- cumsum(rnorm(TT))
  X  <- cbind(x1, x2)
  res <- fqardl(y = y, x = X,
                tau = c(0.25, 0.75),
                pmax = 1L, qmax = 1L, maxk = 1)
  expect_equal(dim(res$beta), c(2L, 2L))
})

test_that("print.fqardl does not error", {
  set.seed(4)
  TT <- 50
  y  <- cumsum(rnorm(TT))
  x  <- cumsum(rnorm(TT))
  res <- fqardl(y = y, x = matrix(x, ncol = 1),
                tau = 0.5, pmax = 1L, qmax = 1L, maxk = 0.5)
  expect_output(print(res))
})

test_that("multipliers data frame is returned", {
  set.seed(5)
  TT <- 60
  y  <- cumsum(rnorm(TT))
  x  <- cumsum(rnorm(TT))
  res <- fqardl(y = y, x = matrix(x, ncol = 1),
                tau = c(0.25, 0.5, 0.75),
                pmax = 2L, qmax = 2L, maxk = 1)
  expect_true(!is.null(res$multipliers))
  expect_true(is.data.frame(res$multipliers))
})