Skip to contents

Compile a DPMirt model specification

dpmirt_compile.Rd

Takes a dpmirt_spec object and compiles the NIMBLE model and MCMC engine. This is the expensive step (~30-120 seconds) that only needs to be done once per model specification.

Usage

dpmirt_compile(
  spec,
  sampler_config = NULL,
  use_centered_sampler = "auto",
  enable_waic = TRUE,
  enable_logprob_monitor = TRUE,
  verbose = TRUE,
  ...
)

# S3 method for class 'dpmirt_compiled'
print(x, ...)

Arguments

spec

A dpmirt_spec object from dpmirt_spec.

sampler_config

Optional advanced hook for NIMBLE sampler customization. Must be NULL or a function with signature function(conf, model, spec). The function receives the configured NIMBLE MCMCconf, the uncompiled NIMBLE model, and the DPMirt spec, and should return the modified MCMCconf or NULL after mutating conf in place. List-based sampler configuration is reserved but not implemented.

use_centered_sampler

Character or logical. Whether to use the centered sampler for SI parameterization. "auto" enables it when appropriate (SI parameterization with 2PL/3PL models).

enable_waic

Logical. Whether to enable WAIC computation.

enable_logprob_monitor

Logical. Whether to add log-probability monitoring samplers.

verbose

Logical. Print progress messages.

...

Additional arguments (currently unused).

x

A dpmirt_compiled object.

Value

A dpmirt_compiled S3 object containing the compiled model, compiled MCMC, specification reference, and session signature.

Details

The compiled NIMBLE objects contain external C++ pointers that cannot be serialized across R sessions. The compile-once pattern is therefore a within-session optimization. For cross-session workflows, save the dpmirt_spec object and recompile in the new session.

See also

dpmirt_spec, dpmirt_sample, dpmirt

Other model fitting: dpmirt(), dpmirt_fit-methods, dpmirt_resume(), dpmirt_sample(), dpmirt_spec()

Examples

if (FALSE) { # \dontrun{
sim <- dpmirt_simulate(200, 20, model = "rasch", seed = 42)
spec <- dpmirt_spec(sim$response, model = "rasch", prior = "normal")

# Compile (takes 30-120 seconds)
compiled <- dpmirt_compile(spec)
print(compiled)

# Compile-once, sample-many pattern
samples1 <- dpmirt_sample(compiled, niter = 5000, nburnin = 1000, seed = 1)
samples2 <- dpmirt_sample(compiled, niter = 5000, nburnin = 1000, seed = 2)
} # }