logmeanexp | R Documentation |
logmeanexp
computes
log mean exp(x_i),
avoiding over- and under-flow in doing so. It can optionally return an estimate of the standard error in this quantity.
logmeanexp(x, se = FALSE)
x |
numeric |
se |
logical; give approximate standard error? |
When se = TRUE
, logmeanexp
uses a jackknife estimate of the
variance in log(x).
log(mean(exp(x)))
computed so as to avoid over- or
underflow. If se = FALSE
, the approximate standard error is
returned as well.
Aaron A. King
## Not run: ## an estimate of the log likelihood: po <- ricker() ll <- replicate(n=5,logLik(pfilter(po,Np=1000))) logmeanexp(ll) ## with standard error: logmeanexp(ll,se=TRUE) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.