Monte Carlo integration against draws from the alternative

Nodes are n draws from the alternative with equal weights, fixed for the fit (reproducible under set.seed()). Certification re-resolves the spec to get an independent sample.

Usage

mc_engine(n)

Arguments

n

Number of draws.

Value

An engine spec, to pass as the engine argument of ripr_init().

Examples

set.seed(1)
fam <- multinomial_family(n_trials = 3L, k = 2L)
null <- null_model(fam, simplex_region(vertices = rbind(c(0.5, 0.5), c(0, 1))))
ripr_init(fam(c(0.6, 0.4)), null, engine = mc_engine(200L))
#> <ripr_state>
#>   null     multinomial_family in 1 part
#>   iterate  1 atom, KL 0.0543246
#>   gap      none recorded
#>   steps    fw 0, lb 0, em 0, weight 0, sfw 0, ct 0 (0.004s)
#>   trace    1 row, 0 snapshots