## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## -----------------------------------------------------------------------------
library(serodynamics)
library(rjags)
library(runjags)
runjags::findJAGS()

## ----"run-dobson-model", warn = FALSE, message = FALSE------------------------

set.seed(1)
data1 <- rbinom(n = 91, size = 1, prob = .6)
jags_post0 <- run.jags(
  n.chains = 2,
  inits = initsfunction,
  model = serodynamics_example("model.dobson.jags"),
  data = list(r = data1, N = length(data1)),
  monitor = "p"
)

## -----------------------------------------------------------------------------
jags_post0$mcmc |> as.array() |> head()