Perform a detailed trace of the comparison process for debugging.
Usage
debug_comparison(pair, y, par, direction = c("A_to_B", "B_to_A"))Examples
exp_spec <- model_spec(
loglik_fn = function(y, lambda) sum(dexp(y, rate = lambda, log = TRUE)),
par_names = "lambda",
par_bounds = list(lambda = c(1e-6, Inf)),
name = "Exponential"
)
gamma_spec <- model_spec(
loglik_fn = function(y, rate) sum(dgamma(y, shape = 1, rate = rate, log = TRUE)),
par_names = "rate",
par_bounds = list(rate = c(1e-6, Inf)),
name = "Gamma(1)"
)
pair <- equivalence_pair(exp_spec, gamma_spec)
set.seed(123)
y <- rexp(50, rate = 2)
trace <- debug_comparison(pair, y, par = c(lambda = 2), direction = "A_to_B")
print(trace)