check_model: Diagnostic panel plot

In tulpa: Templated Unified Library for Posterior Approximation in Bayesian Hierarchical Models

Diagnostic panel plot

Description

Produces a 2x2 or 1x3 panel of diagnostic plots:

1. QQ plot of PIT residuals vs Uniform (with KS p-value)

2. Residuals vs fitted (with lowess smoother)

3. Dispersion histogram (observed variance vs simulated)

4. Spatial correlogram via Moran's I (if coords provided)

Usage

check_model(object, ...)

## Default S3 method:
check_model(object, coords = NULL, nsim = 250L, seed = 123L, ...)

Arguments

object	A fitted model with simulate(), fitted(), residuals()
...	Passed to methods.
coords	Optional N x 2 coordinate matrix for spatial panel
nsim	Number of simulations (default 250)
seed	Random seed (default 123)

Value

Invisible list with ks_p, disp_ratio, moran (if spatial)

Examples

set.seed(1)
d <- data.frame(x = rnorm(80))
d$y <- rpois(80, exp(0.3 + 0.5 * d$x))
fit <- tulpa(y ~ x, data = d, family = "poisson", mode = "laplace")
check_model(fit, nsim = 50)

tulpa documentation built on Oct. 19, 2026, 5:09 p.m.