panel_logmeanexp: Log-mean-exp for panels

In cbreto/panelPomp: Inference for Panel Partially Observed Markov Processes

Log-mean-exp for panels

Description

This function computes the logmeanexp for each column or row of a numeric matrix and sums the result. Because the loglikelihood of a panelPomp object is the sum of the loglikelihoods of its units, this function can be used to summarize replicated estimates of the panelPomp model likelihood. If se = TRUE, the jackknife SE estimates from logmeanexp are squared, summed and the squared root is taken.

Usage

panel_logmeanexp(x, MARGIN, se = FALSE)

Arguments

x	Matrix with the same number of replicated estimates for each panel unit loglikelihood.
MARGIN	The dimension of the matrix that corresponds to a panel unit and over which averaging occurs (1 indicates rows, 2 indicates columns).
se	logical; whether to give standard errors.

Value

A numeric value with the average panel log likelihood or, when se = TRUE, a numeric vector adding the corresponding standard error.

Author(s)

Carles \Breto

See Also

panel_loglik

Examples

ulls <- matrix(c(1,1,10,10),nr=2)
panel_logmeanexp(ulls,MARGIN=2,se=TRUE)

cbreto/panelPomp documentation built on April 13, 2024, 12:23 a.m.