View source: R/likelihoods_origdata.R
mash_compute_loglik | R Documentation |
Compute loglikelihood for fitted mash object on new data.
mash_compute_loglik(g, data, algorithm.version = c("Rcpp", "R"))
g |
A mash object or the fitted_g from a mash object. |
data |
A set of data on which to compute the loglikelihood. |
algorithm.version |
Indicate R or Rcpp version |
The log-likelihood for each element is p(Bhat_j |
Shat_j,g,\alpha)
where Bhat_j | B_j, Shat_j \sim N(B_j,
Shat_j)
and B_j/Shat_j^\alpha | Shat_j \sim g
.
The log-likelihood for data computed using g.
simdata = simple_sims(50,5,1)
data = mash_set_data(simdata$Bhat, simdata$Shat)
m = mash(data, cov_canonical(data))
mash_compute_loglik(m,data)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.