It calculates the marginal density (Lewis and Raftery, 1997) from the DIFM sample using C++.
marginal_d_cpp(data, attributes, hyp_parm, Gibbs, burnin = -1L, verbose = TRUE)
A list of 4 items: Laplace-Metropolis predictive density of the given DIFM, integrated likelihood, the maximum of the predictive densities and determinant of the covariance matrix of the parameters.
The dataset
Model attributes generated from difm.model.attributes
.
Hyperparameters generated from difm.hyp.parm
.
Result of Gibbs sampler from DIFM function.
Burn-in period. If not specified, one tenths of the iterations will be the burn-in period.
Print out the process.