test_that("plot_diagnostics runs without error", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  expect_no_error(plot_diagnostics(fit))
})

test_that("plot_diagnostics works with param subset", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  expect_no_error(plot_diagnostics(fit, params = c("random_error", "twofold_error")))
})

test_that("plot_diagnostics errors on invalid param names", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  expect_error(plot_diagnostics(fit, params = c("nonexistent")), "Unknown")
})

test_that("plot_diagnostics errors for multi", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1,
                          group_by = ~age_group)

  expect_error(plot_diagnostics(fit), "not supported")
})

# ---- plot_trajectory tests ----

test_that("plot_trajectory runs without error", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  expect_no_error(plot_trajectory(fit, id = 1))
})

test_that("plot_trajectory works with custom arguments", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  expect_no_error(plot_trajectory(fit, id = 50, n_samples = 20,
                                  main = "Custom title", show_legend = FALSE))
})

test_that("plot_trajectory errors on invalid id", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  expect_error(plot_trajectory(fit, id = 0), "must be an integer")
  expect_error(plot_trajectory(fit, id = 999999), "must be an integer")
  expect_error(plot_trajectory(fit, id = 1.5), "must be a single integer")
})

test_that("plot_trajectory errors for multi", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1,
                          group_by = ~age_group)

  expect_error(plot_trajectory(fit), "not supported")
})

# ---- plot_trajectory subject_ids tests ----

test_that("plot_trajectory works with character subject_ids", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  n <- nrow(inputdata)
  char_ids <- paste0("S", seq_len(n))

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1,
                          subject_ids = char_ids)

  # Character lookup should work
  expect_no_error(plot_trajectory(fit, id = "S1"))
  # Row index should still work
  expect_no_error(plot_trajectory(fit, id = 1))
})

test_that("plot_trajectory works with numeric subject_ids out of range", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  n <- nrow(inputdata)
  num_ids <- 10000L + seq_len(n)  # e.g. 10001, 10002, ...

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1,
                          subject_ids = num_ids)

  # Numeric ID out of row range should trigger lookup
  expect_no_error(plot_trajectory(fit, id = 10001))
  # Row index should still work
  expect_no_error(plot_trajectory(fit, id = 1))
})

test_that("plot_trajectory works with subjects argument", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  n <- nrow(inputdata)
  ext_ids <- paste0("EXT", seq_len(n))

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  # subjects argument should enable lookup
  expect_no_error(plot_trajectory(fit, id = "EXT1", subjects = ext_ids))
})

test_that("plot_trajectory errors on nonexistent subject_id", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  n <- nrow(inputdata)
  char_ids <- paste0("S", seq_len(n))

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1,
                          subject_ids = char_ids)

  expect_error(plot_trajectory(fit, id = "NONEXISTENT"), "not found")
})

test_that("plot_trajectory errors on out-of-range numeric without subject_ids", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  # No subject_ids → out of range numeric should error
  expect_error(plot_trajectory(fit, id = 99999), "must be an integer")
})

test_that("plot_trajectory works in multi-panel layout", {
  data("inputdata", package = "seroreconstruct")
  data("flu_activity", package = "seroreconstruct")

  fit <- sero_reconstruct(inputdata, flu_activity,
                          n_iteration = 200, burnin = 100, thinning = 1)

  par(mfrow = c(1, 2))
  expect_no_error({
    plot_trajectory(fit, id = 1, show_legend = FALSE)
    plot_trajectory(fit, id = 50, show_legend = FALSE)
  })
  # Verify mfrow was preserved (not reset)
  expect_equal(par("mfrow"), c(1, 2))
  par(mfrow = c(1, 1))
})