This function generates Laplace samples of posterior distribution for a vector of reported incidence.
sample_laplace_log_incidence_poisson(
beta_hat,
beta_cov,
reported,
Q,
num_samps_per_ar = 10
)
Maximum likelihood solution for beta parameter.
Covariance of objective solution (either Fisher information or Hessian inverse).
An integer vector of reported cases.
Spline basis matrix.
Number of Laplace samples to return for each AR path.
A matrix of `num_samps_per_ar` log incidence curve samples from laplace approximation of distribution.