set.seed(123)
data("data", package = "ameras")

test_that("vcov & plot & summary_table binomial", {
  skip_on_cran()

  fit <- ameras(
    Y.binomial ~ dose(V1:V10) + X1 * X2,
    data = data,
    family = "binomial",
    methods = c("RC", "FMA")
  )

  expect_no_error({
    vcov(fit)
    vcov(fit, methods = "FMA")
    summary(fit)
    summary_table(fit)
    plot(fit, ask = FALSE)
    plot(fit, methods = "RC", ask = FALSE)
  })

  expect_no_error(
    plot(fit, which = "qq", ask = FALSE)
  )

  expect_no_error(
    plot(fit, type = "deviance", ask = FALSE)
  )
})

test_that("plot works with supplied data when keep.data is FALSE", {
  dat <- data
  fit <- ameras(
    Y.binomial ~ dose(V1:V10),
    data = dat,
    family = "binomial",
    methods = "RC",
    keep.data = FALSE
  )

  # plot() resolves the supplied data once and reuses that resolved version
  # while computing residuals; the caller's data frame should stay untouched.
  expect_no_error(plot(fit, data = dat, ask = FALSE))
  expect_false("rcdose_ameras" %in% names(dat))
})


test_that("vcov & plot & summary_table multinomial", {
  fit <- ameras(
    Y.multinomial ~ dose(V1:V10) + X1 * X2,
    data = data,
    family = "multinomial",
    methods = "RC"
  )

  expect_no_error({
    vcov(fit)
    summary_table(fit)
    plot(fit, ask = FALSE)
  })

  expect_no_error(
    plot(fit, which = "qq", ask = FALSE)
  )

  expect_no_error(
    plot(fit, type = "deviance", ask = FALSE)
  )
})


test_that("vcov & plot & summary_table prophaz", {
  fit <- ameras(
    Surv(time, status) ~ dose(V1:V10) + X1 * X2,
    data = data,
    family = "prophaz",
    methods = "RC"
  )

  expect_no_error({
    vcov(fit)
    summary_table(fit)
    plot(fit, ask = FALSE)
  })
})