stephenslab / ebnm

R package to fit Empirical Bayes Normal Means model.
https://stephenslab.github.io/ebnm/
12 stars 9 forks source link

Likelihood in summary? #102

Closed pcarbo closed 1 year ago

pcarbo commented 1 year ago
> summary(fit_normal)
Call:
ebnm_normal(x = x, s = s, mode = "estimate")

EBNM model was fitted to 688 observations with _heteroskedastic_ standard errors.

The fitted prior is of class _normalmix_.

2 degrees of freedom were used to estimate the model.
The likelihood is 989.59.

Available posterior summaries: _mean_, _sd_.
Use method fitted() to access available summaries.

A posterior sampler is _not_ available.
One can be obtained using function ebnm_add_sampler().

The summary here says the likelihood is 989.59, but isn't it the log-likelhood?

willwerscheid commented 1 year ago

whoops, yes you are correct

willwerscheid commented 1 year ago

fixed