mcmc_temporal_catmark {stpphawkes}

Bayesian Estimation of Temporal Hawkes Model Parameters with Categorical Marks

Description

This function computes the posterior of the parameters of a temporal exponential decay Hawkes model using Metropolis-with-in-Gibbs sampling.

Usage

mcmc_temporal_catmark(
  times,
  marks,
  t_max = max(times),
  t_mis = NULL,
  param_init = NULL,
  mcmc_param = NULL,
  branching = TRUE,
  print = TRUE
)

Arguments

times

- vector of arrival times

marks

- vector of marks

t_max

- maximum time value (default = max(times))

t_mis

- mx2 matrix, mth row contains two elements describing the mth missing time range (default = NULL)

param_init

- list of parameters of initial guess (default = NULL, will start with MLE)

mcmc_param

- list of mcmc parameters

branching

- using branching structure in estimation (default = TRUE)

print

- print progress (default = TRUE)

Details

The default is to estimate the branching structure which is much more computationally efficient. The model will also account to missing data if t_mis is provided.

Value

A DataFrame containing the mcmc samples