Learn R Programming

celda (version 0.0.0.9000)

calculate_marginal_likelihood: Calculate the marginal likelihood from a single celda model

Description

Marginal likelihood is estimated as the harmonic mean of the (non-log) likelihood over all iterations of Gibbs sampling.

Usage

calculate_marginal_likelihood(completeLogLik)

Arguments

completeLogLik

The complete Gibbs sampling history of log-likelihoods for a single celda chain

Value

The estimated marginal likelihood as an mpfr number