Run all stationarity-related tests for INAD

Usage

run_stationarity_tests_inad(
  y,
  order = 1,
  thinning = "binom",
  innovation = "pois",
  blocks = NULL,
  verbose = FALSE,
  ...
)

Arguments

y: Integer matrix.
order: Model order (1 or 2).
thinning: Thinning operator.
innovation: Innovation distribution.
blocks: Optional block assignments.
verbose: Logical.
...: Additional arguments.

Value

A list with class "stationarity_tests_inad".

Examples

set.seed(1)
y <- simulate_inad(
  n_subjects = 25,
  n_time = 5,
  order = 1,
  thinning = "binom",
  innovation = "pois",
  alpha = 0.3,
  theta = 2
)
out <- run_stationarity_tests_inad(
  y,
  order = 1,
  thinning = "binom",
  innovation = "pois",
  verbose = FALSE,
  max_iter = 15
)
out$summary
#>       constraint method statistic df      LRT   p_value BIC_selected
#> alpha      alpha    lrt  3.099178  3 3.099178 0.3765853  constrained
#> theta      theta    lrt  1.468989  4 1.468989 0.8321201  constrained
#> both        both    lrt  4.245761  7 4.245761 0.7510753  constrained