prior {pompp}R Documentation

Build a joint prior for pompp model parameters

Description

Constructor for pompp_prior objects, which is used in the pompp_fit function. The generated prior is so that Beta, Delta and LambdaStar are indepdendent a priori.

Usage

prior(beta, delta, lambdaStar, marksMean, marksPrecision)

Arguments

beta

An S4 object whose class inherits from BetaDeltaPrior.

delta

An S4 object whose class inherits from BetaDeltaPrior.

lambdaStar

An S4 object whose class inherits from LambdaStarPrior.

marksMean

An S4 object of class NormalPrior.

marksPrecision

An S4 object of class GammaPrior.

Value

A pompp_prior object with the adequate slots. It is ready to be included in a model via the pompp_model function.

See Also

fit_pompp, NormalPrior, GammaPrior and pompp_model.

Examples

# Let us say there are 3 intensity covariates and 4 observability covariates.
# One more element is included in both sets due to the intercepts.
new_prior <- prior(
  NormalPrior(rep(0, 4), 10 * diag(4)),
  NormalPrior(rep(0, 5), 10 * diag(5)),
  GammaPrior(0.0001, 0.0001),
  NormalPrior(0, 100), GammaPrior(0.001, 0.001)
)

[Package pompp version 0.1.3 Index]